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Dear Colleagues, 


Aircraft design is, as we know, the first fascinating step in the life of an aircraft, where visions 
are converted into reality. 


In a practical sense, aircraft design supplies the geometrical description of the aircraft. 
Traditionally, the output is a three-view drawing and a list of aircraft parameters. Today, the 
output may also be an electronic 3D model. In the case of civil aircraft, a fuselage cross- 
section and a cabin layout are provided in addition. 


In an abstract sense, aircraft design determines the design parameters to ensure that the 
requirements and constraints are met and design objectives are optimized. The fundamental 
requirements for civil aviation are payload and range. Many constraints come from 
certification rules demanding safety. The objectives are often of a financial nature, like 
lowest operating costs. Aircraft design always strives for the best compromise among 
conflicting issues. 


The design synthesis of an aircraft goes from the conceptual design to the detailed design. 
Frequently, expert knowledge is needed more than computing power. Typical work involves 
Statistics, the application of inverse methods, and use of optimization algorithms. Proposed 
designs are analyzed with respect to aerodynamics (drag), structure (mass), performance, 
Stability and control, and aeroelasticity, to name just a few. A modern aircraft is a complex, 
computer-controlled combination of its structure, engines, and systems. Passengers demand 
high comfort at low fares, society demands environmentally friendly aircraft, and investors 
demand a profitable asset. 


Overall aircraft design (OAD) comprises all aircraft types in civil and military use, considers all 
major aircraft components (wing, fuselage, tail, undercarriage) as well as the integration of 
engines and systems. The aircraft is seen as part of the air transport system and beyond 
contributing to multimodal transport. Aircraft design applies the different aerospace 
sciences and considers the aircraft during its whole life cycle. Authors from all economic 
sectors (private, public, civic, and general public) can submit to this Special Issue (Sl). 


Education and training in aircraft design is considered as important as research in the field. 


The SI can be a home for those active in the European Workshop on Aircraft Design 
Education (EWADE) or the Symposium on Collaboration in Aircraft Design (SCAD), both 
independent activities under the CEAS Technical Committee Aircraft Design (TCAD). Please 


see http://AircraftDesign.org for details. 


Following the successful initial Special Issue on “Aircraft Design (SI-1/2017)” and the 
relaunch with “Aircraft Design (SI-2/2020)”, this is already the third SI in sequence named 
“Aircraft Design (SI-3/2021)”. The editorial “Publishing in ‘Aircraft Design’ with a Continuous 
Open Access Special Issue” describes the history, the set up, and idea behind this SI. It was 
published on 14 Jan 2020 as https://doi.org/10.3390/aerospace7010005. 


Activities in the past showed that aircraft design may be a field too small to justify its own 
(subscription-based) journal. A continuous open access special issue may fill the gap. As 
such, the Special Issue “Aircraft Design” can be a home for all those working in the field who 
regret the absence of an aircraft design journal. 


The Special Issue “Aircraft Design" is open to the full range of article types. It is a place to 
discuss the "hot topics" (zero-emission airplanes, electric flight, urban air mobility — you 
name it). The classic topics in aircraft design remain: 

e Innovative aircraft concepts 

e Methodologies and tools for aircraft design and optimization 


e Reference aircraft designs and case studies with data sets 


4 


It is up to us as authors to shape the Special Issue “Aircraft Design” according to our interests 
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Abstract: During its investigations into a series of ten aircraft crashes from 1979 to 1981, US National 
Transportation Safety Board (NTSB) officials were presented with a hypothesis that “several” of the 
crashes could have been caused by pilot impairment from breathing oil fumes inflight. The NTSB and 
their industry partners ultimately dismissed the hypothesis. The authors reviewed the crash reports, 
the mechanics of the relevant engine oil seals, and some engine bleed air data to consider whether 
the dismissal was justified. Four of the nine aircraft crash reports include details which are consistent 
with pilot impairment caused by breathing oil fumes. None of the tests of ground-based bleed air 
measurements of a subset of oil-based contaminants generated in the engine type on the crashed 
aircraft reproduced the inflight conditions that the accident investigators had flagged as potentially 
unsafe. The NTSB’s conclusion that the hypothesis of pilot incapacitation was “completely without 
validity” was inconsistent with the evidence. Parties with a commercial conflict of interest should 
not have played a role in the investigation of their products. There is enough evidence that pilots 
can be impaired by inhaling oil fumes to motivate more stringent design, operation, and reporting 
regulations to protect safety of flight. 


Keywords: aircraft; accident; flight safety; fumes; engine oil; hydraulic fluid 


1. Introduction 


With the exception of the Boeing 787, commercial and military aircraft are designed to 
“bleed” (or extract) ventilation air, either off the main aircraft engine compressors or from 
an auxiliary compressor when the aircraft is on the ground. This hot compressed “bleed 
air” is then cooled, dehumidified, typically mixed with some fraction of recirculated air, 
and routed to the cabin and flight deck for ventilation and pressurization [1]. To greater or 
lesser degrees—engine oil can contaminate the compressed air, whether it migrates across 
engine seals, spills from an overserviced reservoir, or is vented improperly, for example [2]. 
The problem with this design is that a fraction of that compressed air is then bled off the 
engine and routed to the air conditioning system for cabin ventilation and pressurization. 
The “bleed air” stream is not filtered, so oil fumes from the engine can be delivered directly 
to the occupied zones of the aircraft. 

Concerns about aircrew exposure to either suspected or confirmed bleed-sourced 
fumes inflight have been raised globally since the 1930s on all types of aircraft [3-13]. One 
of the early references to pilots breathing “hot oil fumes” inflight acknowledged that “the 
symptoms in these cases have been similar to those of carbon monoxide poisoning,” but 
noted that various aldehyde breakdown products are “probably the causative agents,” 
which the author describes as present at concentrations that are “obviously sufficient to be 
dangerous to safety inflight” [14] (pp. 178-180). 

Within this history, there are two documents which explicitly refer to fatal accidents 
associated with pilot impairment caused by oil-contaminated bleed air. One of those 
documents refers to “several unexplained fatal crashes involving single pack, carrier-based, 
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turbine-powered aircraft ... attributed (rightly or wrongly) to contaminated engine bleed 
air” [15] (p. 1). In response to those crashes, the United States (U.S.) Air Force, Navy, and 
manufacturers apparently “initiated efforts to resolve the contamination problem” [15] 
(p. 1). The other document describes an investigation by the U.S. National Transportation 
Safety Board (NTSB) conducted in response to ten crashes on privately-operated turbo- 
prop planes, all equipped with the Garrett TPE331 engine [16]. NTSB investigators had 
hypothesized that “several” of the crashes could have been the result of the pilots being 
impaired by breathing oil fumes inflight. The suggested source of the oil fumes was a 
cracked front main shaft compressor carbon seal in the engine. In response, the NTSB 
initiated ground-based engine testing to assess what the pilots may have inhaled [16]. In 
conjunction, the U.S. Federal Aviation Administration (FAA) exposed test animals to oil 
fumes and observed the physiological effects [17]. Ultimately, the NTSB and their industry 
partners dismissed the hypothesis that inhaling oil fumes could have impaired any of the 
pilots [16]. The authors reviewed the evidence to determine if the dismissal was supported 
by the evidence that was available at the time. The authors also reviewed the relevant 
aviation regulations in the United States to consider whether they are sufficient to prevent 
inflight exposure to oil fumes. 


2. Materials and Methods 


Between August 1979 and April 1981, ten private aircraft, each outfitted with Garrett 
TPE331 engines lubricated with Exxon 2380 oil, crashed. The dates, location, aircraft types, and 
registration numbers are listed in Table 1, reproduced from the NTSB investigative reports [16]. 


Table 1. Accidents investigated for pilot incapacitation. 


Crash No. ae dd) Location Aircraft Type oo 
1 1979-08-03 Hays, Kansas Mitsubishi MU-2B N208MA 
2 1979-11-01 Nashville, Tennessee Mitsubishi MU-2F N8730 
3 1979-12-21 Provo, Utah Mitsubishi MU-2B-20 N2-OBR 
4 1980-01-11 Atlantic Ocean Cessna 441 N441NC 
5 1980-02-14 Near Houston, Texas Mitsubishi MU_2B-35 N346MA 
6 1980-02-23 A Retisale Mitsubishi MU-2-40 N962MA 

Louisiana 
7 1980-12-06 Ramsey, Minnesota Mitsubishi MU-2-40 N969MA 
8 1980-12-15 Richmond, Indiana Mitsubishi MU-2B-30 N93UM 
9 1981-01-07 Burns, Oregon Aero Commander 690B N81521 
10 1981-04-22 Alpena, Michigan Mitsubishi MU-2B-20 NYS 


The authors searched online NTSB databases for the crash reports, nine of which 
were fatal accidents and one of which was a serious incident. The crash reports were only 
available on microfilm and, because of the COVID-19 pandemic, the NTSB library which 
houses microfilm reports was closed. Instead, the authors were able to purchase nine of 
the ten crash reports through a private company (General Microfilm, West Virginia). The 
remaining report was unavailable. Additional details of the 10 crashes are listed in Table 2. 

The authors researched the design of the Garrett TPE331 engine and the associated 
potential for bearing and seal failures which could cause oil fumes to contaminate the bleed 
air supply. 

The authors reviewed the NTSB investigative report which described the results of 
bleed air quality engine testing for selected oil-based contaminants in the bleed air supplied 
by a Garrett TPE331 engine with and without oil contamination [16]. The test protocol 
had been designed to “investigate the hypothesis that toxic or anesthetic gases could be 
generated from engine oil that leaked into the engine airflow through a broken seal and 
that these gases might adversely affect the crew’s capacity to control the aircraft during 
critical phases of flight” [16] (p. 8). Key features of the bleed air tests are listed in Table 3. 
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Table 2. Additional details of ten accidents investigated for pilot incapacitation. 
NTSB Report Flight Phase, Time of Pilot Age (Years), Probable Causes/Factors ee 
CHSNG: No. Day, Weather Flight Time (h) (per NTSB) No: Deets 
Forward main shaft bearing 
me ‘e , 34; failed; improper inflight 
: BOSC (EEO cTuise, night, “clear 2168. decisions; failed to obtain or ss 
maintain flying speed. 
PIC age redacted; Pilot-in-command (PIC) 
descent, night, 1498. misjudged the distance and 
2 TAD80FA007 three-mile visibility, dry, Copilot age altitude on final approach; 5 
light fog redacted; inadequate supervision of 
VSPA flight; failed to use checklist. 
Inadequate pre-flight 
preparation or planning; 
improper Instrument Flight 
; tek wee «| Rules (FR) operations 
y ae 12,833 misjudged altitude; Pilot 
one mile 
could not find runway; 
aircraft came to rest in 
water. 
Miscellaneous, 
undetermined; unable to 
- 47; obtain a response from crew; 
4 IAD80A A018 descent, night, check 8000 fmeantre lied descent 3 
aircraft came to rest in 
water. 
approach, evening, ss fe 
5 FTW80FA042 —- two-mile visibility, rain oe pO ee 4 
12,500 altimeter setting incorrect. 
showers 
ee approach, morning, ol; Improper IFR operation; 
: Saas 0.25-mile visibility, fog hrs. not reported crashed into water. . 
Failed to obtain/maintain 
J approach, afternoon, 54; speed; icing conditions, 
ae “good” weather 4949 including sleet, frozen rain, 
etc. 
; Improper IFR operation; 
8 CHigireces: ““PPTosthy ences ee icing conditions, including 0 
fog 10,587 
sleet, frozen rain, etc. 
approach, afternoon, 39: Miscellaneous /undetermined; 
OFF SEA81FA015 overcast, 15-mile 01977 uncontrolled collision with 2. 
visibility, wind calm eround / water. 
' Improper IFR operation; 
10 CHI81FA051 rec pen ae ht a crash on fina approach, 3 
eink oe 18 16,766 1.6 miles short of runway; 


cause unknown. 


* Microfilm report was unavailable for crash #4. ** Details in the reports for crash #1, 6, 7, and 9 are consistent with engine oil fumes as a 
contributory factor. Additional details are provided in Table 5. *** In all but crash #8, everyone onboard was killed. In crash #8, the pilot 
and one of two passengers were injured. 


Aerospace 2021, 8, 389 4 of 20 


Table 3. Key features of engine bleed testing, per the NTSB investigative report [16]. 


Reflects Potential 





Carbon Seal Intact Oil Added to Glass Filter in Bleed Air 
Test No. a ; Onboard 
or Missing? Bleed Air? Line? i 
Conditions? 
1 intact No yes no 
Z intact No yes no 
3a—3f intact yes; 2-12 lb/hr yes no 
3g intact No no no 
4 intact No yes no 
5a missing No yes no 
5b missing yes; “dirty start” yes no 


The authors reviewed the relevant aviation regulations for U.S.-registered aircraft 
operated on routine commercial flights (Table 4). Globally, aviation regulations are harmo- 
nized, so most of these regulations are the same in other countries. This is especially true 
for design and construction-related regulations which are classified under Part 25 of the 
Federal Aviation Requirements in the United States and as CS-25 (Certification Standards) 
in Europe, published by the European Union Aviation Safety Agency (EASA). 


Table 4. U.S. aviation regulations relevant to air quality on routine commercial flights. 


Regulation Design/Operation Description 





“The aircraft air supply system must be 


14 CFR § 25.831(b) * designed to ensure that crew and passenger 








ae Desi — 
“Ventilation” ane compartment air [is] free from harmful or 
hazardous concentrations of gases or vapors.” 
The carbon monoxide concentration in the cabin 
es si Design and flight deck supply air must not exceed 
25.831(b)(1)” Ventilation” 
50 ppm. 
“Warning information must be provided to alert 
the crew to unsafe system operating conditions, 
14 CFR § 25.1309(c) and to enable them to take appropriate 
“Equipment, systems, Design corrective action. Systems, controls, and 
and installation” associated monitoring and warning means must 
be designed to minimize crew errors which 
could create additional hazards.” 
Airlines “shall report the occurrence or detection 
14 CER § 121.703(a)(5) of each failure, malfunction, or defect concerning 
- nes ... [a] aircraft component that causes 
Service difficulty Operation ; é 
ance! accumulation or circulation of smoke, vapor, or 
P toxic or noxious fumes in the crew compartment 
or passenger cabin during flight.” 
Airlines shall report each “interruption to a 
14 CFR § scheduled flight,” such as a diversion, 
121.705“Mechanical cancellation, or tail swap, caused by known or 
Operation eve 
interruption summary suspected mechanical difficulties or 
report” malfunctions that are not required to be reported 


under the 14 CFR § 121.703. 


Airlines shall report accidents and occurrences 
which are associated with the operation of an 
aircraft and affect (or could affect) the safety of 
operation, including smoke/fumes. 


FAA Order 8020.11D 
(Chapter 6 and FAA Operation 
Form 8020-23) 


* Note that “14 CFR” refers to the Code of Federal Regulations Title 14 which are aviation and aerospace regulations 
published by the FAA for U.S.-registered aircraft. The symbol “§”is shorthand for “Part” when referring to U.S. 
regulations. The relevant regulations include “Part 25” (Airworthiness Standards: Transport Category Airplanes) 
and “Part 121” (Operating Requirements: Domestic, Flag, and Supplemental Operations). 
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Finally, the authors reviewed some proposals to implement engineering and opera- 
tional control measures intended to prevent exposure to oil fumes on aircraft. 


3. Results 
3.1. Accident Reports 


For nine of the ten aircraft crashes, the authors were able to obtain the NTSB investi- 
gation reports (Tables 1 and 2). The microfilm report for crash number 4 was unavailable. 
Of these nine crashes, eight were fatal accidents (35 people were killed) and one (crash 
number 8) was a serious incident (the pilot and one of the two passengers were injured). 
Details in the accident reports for crash numbers 1, 6, 7, and 9 (Tables 1 and 2) suggest that 
oil fumes were either a contributory or casual factor (Table 5). 


Table 5. NTSB accident report details consistent with exposure to oil fumes. 


Crash No. Description of Crash Conditions 


The NTSB report describes how, when the pilot was en route, at what may have 
been top of climb, the pilot “reported a loss of oil pressure in the right engine.” 
At the same time, “the pilot reported smoke and fumes in the cabin.” A few 

1 minutes later, the pilot reported that he had shut down the right engine. He then 
stopped communicating with air traffic control and crashed into a field without 
putting the landing gear down (per witness reports). The subsequent teardown 
of the right engine “revealed that the forward main shaft bearing had failed.” 





The NTSB report describes how, on approach, the pilot stopped responding to 
air traffic control. The aircraft descended into a lake and “all aboard perished.” 
The report goes on to say that “[d]uring the course of the investigation, it was 
reported by various persons that the aircraft had a history of smoke, fumes, and 
carbon monoxide collecting in the cockpit and cabin area...Examination of the 
interior of the subject aircraft environmental control system revealed an oily 
residue in the portion which is supplied by the left engine [with a history of a 
cracked carbon seal] while that supplied by the right engine was clean.” 


The NTSB report describes how, on approach, the pilot did not respond to air 
traffic control instructions. Witness reports described how the aircraft came out 
of the clouds spinning with the nose down and crashed into a field. In the 

7 engine tear down report, there were references to “black dirt deposits found 
throughout the entire gas path of the engine ... [including] the bleed air ports 
... ” The NTSB report stated that the pilot had failed to obtain/maintain speed 
and the cause of the crash was unknown. 


The NTSB report describes how ground witnesses observed the aircraft flying 
away from the airport maneuvering erratically. The aircraft crossed a highway, 
pitched up steeply, fell, and crashed. The engine tear down was “to determine 
the condition of the carbon seal and bearing located between the compressor 

9 section and gearbox section.” The report concluded that “one carbon seal, due to 
the discoloring of the oil slinger, was leaking some oil past the face.” The report 
also noted that “Bonneville Power maintenance personnel felt the amount of 
leakage was insignificant,” but the basis for their claim and their qualification to 
make it were not provided. 


In addition to the four crashes described in Table 5, crash numbers 5 and 10 
(Tables 1 and 2) include details which could be related to oil fumes, although the sig- 
nificance of those details is less clear. Specifically, the NTSB report for crash number 
5 describes how, without warning, the aircraft crashed into the trees near the runway on 
final approach. Another pilot who observed the crash stated: “We feel that the pilot of 
the MU-2 got disoriented and just flew the airplane into the ground.” The autopsy report 
for the pilot reported 5% carboxyhemoglobin in his blood. The NTSB report for crash 
number 10 describes how the aircraft crashed on final approach, after which the engine 
manufacturer performed the engine tear down to determine if oil fumes played a role in the 
crash. The manufacturer concluded that there was “no indication” of anything other than 
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normal engine operation and they attributed the engine compressor carbon seal damage to 
impact. The Civil Aeromedical Institute (CAMI) toxicology lab (the research arm of the FAA) 
tested a component of the air supply system from the aircraft and found a “trace quantity of 
petroleum base constituent,” although the significance of that finding is not clear. 


3.2. Description of the TPE331 Engine and Seal Assembly 


The Garrett TPE331 is a fixed-shaft constant-speed turboprop engine (Figure 1). Its gas 
turbine consists of a compressor, combustion chamber, and turbine. Ambient air is directed 
to the compressor section through the engine inlet. A two-stage centrifugal compressor 
increases air pressure and directs it to the combustion chamber. In the combustion chamber, 
fuel is added to the air through the fuel nozzles. On engine start-up, the gas mixture 
is ignited by igniter plugs. In a normal operation the igniter is not in use because the 
combustion is self-sustained. The hot and high-velocity combustion gases flow through 
the turbine rotors, where the energy of the gases is converted to torque exerted on the main 
shaft (A in Figure 1). The reduction gear is designed as a planetary gear on the propeller 
shaft (B in Figure 1). It converts the low torque (at high rpm) of the main shaft to high 
torque (at low rpm) on the propeller shaft and drives the propeller. The combustion gases 
leave the turbine to the atmosphere via the exhaust. 


Exhaust 
outlet 


Figure 1. Garrett TPE331 turboprop engine (based on [18], p. 15-3). A. main shaft (engine shaft) with the gas turbine. B. 


propeller shaft with the reduction gear. 1,2. bearings that support the main shaft. 3,4. bearings that support the propeller 


shaft. 1. compressor bearing. 


The TPE331 is called a fixed-shaft engine because the propeller is firmly connected 
to the gas turbine. The constant-speed engine maintains its speed by a governor on the 
propeller. The propeller shaft rotates at a constant speed of 1591 rpm in cruise flight. The 
main shaft of the engine rotates at a constant speed of 41,730 rpm. Power changes are made 
by increasing the fuel flow (which increases the torque) rather than the engine speed. 

Most of the air (70%) passing through the engine provides internal cooling. Only 
about 10% of the air that passes through the engine is actually used in the combustion 
process. Up to approximately 20% of the compressed air may be bled off for the purpose 
of heating, cooling, cabin pressurization, and pneumatic systems [18], which appears to 
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be within the average range for other engine types [2]. If the bleed air is contaminated 
with oil fumes, then the air in the cabin is also contaminated. Normal bleed air pressure is 
approximately 157 psi (10.82 bar) and the temperature is 360 °C in cruise flight [16]. 

The engine shafts are supported by rolling bearings. The compressor bearing (1 in 
Figure 1) is a ball-bearing (Figure 2). The same is true for bearings 2 and 3 in Figure 1. 





Figure 2. A typical ball-bearing [19]. 


The engine cutaway drawing is given in Figure 3 and a detail of that picture enlarged 
is given in Figure 4. 


A 


Figure 3. Garrett TPE331 turboprop engine cutaway drawing [20]. 


In the TPE331 front main shaft compressor seal assembly, two seals are used, acting in 
series. They prevent the escape of the engine lubricant present in the reduction gearbox via 
the main shaft and through the compressor bearing into the compressor (Figure 4). The 
following detailed explanations are reproduced with text elements from [16]: 

“First, a mechanical (carbon) seal is provided to prevent leakage if operating oil pressures 
exceed normal pressures and to prevent potential leakage that might occur during engine 
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shutdown. The mechanical seal potentially is subject to wear and damage” [16] (p. 2). “When 
the engine is shutdown, oil drains from the engine walls which increases the oil level above 
the shaft installation level of the compressor front main shaft; without a mechanical seal, 
the oil would leak from the reduction gearbox, out the engine inlet and onto the ground. 
Therefore, a mechanical (carbon-element) seal is necessary to prevent oil leakage when the 
engine is not operating” [16] (p. 33). This is illustrated in Figure 5. 


~ 
A \ 


Figure 4. Garrett TPE331 turboprop engine cutaway drawing (detail based on [20]). A. main shaft. 





B. first-stage centrifugal compressor. 1. compressor bearing (a ball-bearing). 2. carbon seal. 3. labyrinth seal. 


Accessory 
Gearbox 
Housing 













Washer 
Compressor 
Bearing O-Ring 
Spring 
Retainer 


Shaft 


a % om ee ee ee 


Figure 5. Mechanical (carbon) seal assembly next to the main shaft compressor bearing on the 
TPE331 engine (based on Figure 3 in [16]). 


“Secondly, the TPE331 has a pressurized knife-edge labyrinth air seal that is specifically 
designed to prevent passage of air/oil mist from the reduction gearbox into the compressor 
chamber during engine operation” [16] (p. 2). “This pneumatic-type seal is pressurized to 
approximately 26 psi (1.79 bar) at the inner knife edges. When the engine is operating, the 
pneumatic seal is independently capable of preventing the air/oil mist from passing out of 
the reduction gearbox because of a flow of pressurizing air from the seal into the reduction 
gearbox. However, the labyrinth seal has no sealing capability to prevent engine oil from 
exiting the reduction gearbox when the engine is not operating and, therefore, requires a 
mechanical seal to prevent oil leakage when the engine is shut down. The symptom of 
a failed mechanical seal is oil leaking past the labyrinth seal and running out the engine 
inlet onto the ground when the engine is not operating” [16] (p. 35). The bleed air passage 
directs sealing air from the compressor to the labyrinth seal. All of this is illustrated in 
Figure 6 and additional details are shown in Figure 7. 
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Figure 6. Labyrinth seal close to the main shaft compressor bearing (based on Figure 11 in [16]). Note: The carbon seal is located in 


between the compressor bearing (on the left) and the labyrinth seal (on the right) in the empty space but is not shown in this figure. 


AIR/OIL MIST 
TURNED AWAY 





SUPPLY 


BLEED AIR 


28 psi 


Figure 7. Flow through the labyrinth seal in the TPE331 engine next to the compressor bearing and resulting pressure. The 


flow to the left pushes back the oil that still comes through the carbon seal (based on Figure 13 in [16]). 


“The pressure balances that are maintained may be clearly understood by considering 
that the chamber containing the air/oil mist operates at approximately 12 psia (—2 psig) 
while the center of the labyrinth is pressurized to 26 psia (+12 psig). Since air will always 
flow from a higher pressure area to a lower pressure area, air will flow outward from the 
center of the labyrinth. Each of the knife-edges operates with only slight clearance from 
the outer wall, resulting in a high local velocity as the air crosses into the next chamber. 
With a high local velocity, the air/oil mist cannot flow past the knife edges and enter the 
compressor air flow” [16] (p. 36). This is illustrated in Figure 7. 

“In addition to the labyrinth seal, there are two other pneumatic sealing actions 
incorporated in the oil containment design of the TPE331 engine ... [First, there is] a 
negative pressure of approximately 2 psi inside the reduction gearbox with respect to 
atmospheric conditions. This is accomplished by pumping the oil out of the reduction 
gearbox at a rate twice that at which it is being pumped in. (Nine gallons per minute in, 
versus 18 gallons per minute out.) The extra nine gallons pumped per minute is air and 
this creates a negative pressure (slight vacuum) in the reduction gearbox. If the mechanical 
seal should fail, the overscavenging will draw air into the reduction gearbox and will 
prevent the air/oil mist from flowing out ... Secondly the negative pressure effect is 
supplemented by raising the pressure in the chamber on the air side of the mechanical seal 
to a positive six psi above ambient” [16] (p. 36). This higher pressure is produced by the 
centrifugal compressor. 
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3.3. Bleed Air Testing 


In July 1981, the NTSB and a team of government and industry partners formed an 
“Ad Hoc Investigative Committee” to respond to the hypothesis that oil fumes could con- 
taminate pilots’ breathing air through the main shaft compressor carbon seal in the Garrett 
TPE331 turboprop engine. All the crashed aircraft flew with this engine type. The team 
proposed to measure selected gaseous contaminants (carbon monoxide, carbon dioxide, 
nitrous oxides, and total hydrocarbons) in the bleed air produced by that engine type. 
The goal was to attempt to reproduce the conditions on the crashed aircraft to determine 
whether the pilots could have been subjected to toxic or anesthetic gases sufficient to impair 
their ability to operate the aircraft. 

Because oil leakage through the carbon seal was the suspected source of oil fumes in 
the bleed air, selected gaseous contaminants were measured in the bleed air downstream of 
the compressor, first through an intact seal and second with the seal removed (considered 
worst-case), during conditions that ranged from no oil to 12 pounds of oil per hour (lb/hr). 

The sampling methods were developed by the Exxon Research and Engineering 
Company and the test conditions, summarized in Table 3, were performed on an engine 
test stand at the Garrett Turbine Engine Company. During two of the tests, oil mist 
samples were also collected and analyzed. Largely though, gases were sampled because 
the investigators “expected that particulate matter would not form” [16] (p. 9). 


3.3.1. Primary Bleed Air Testing Conditions 


The key features of the 12 primary engine bleed testing conditions, as described in the 
NTSB report [16] are summarized in Table 3. Additional salient details of the engine bleed 
testing are summarized in Table 6. 


Table 6. Description of TPE331 engine bleed testing, per NTSB report [16]. 


Trial No. Description of Sampling Conditions 





These two trials were intended to characterize background levels of contaminants in 
the bleed air stream. The engine compressor carbon seal was intact, and a glass wool 
filter was installed in the bleed line. Moreover, oil was not purposefully injected into 
the bleed air stream. 





These six trials were intended to characterize levels of oil-based contaminants in the 
bleed air stream through an intact engine compressor carbon seal and a glass wool 

3a—3f filter. Oil was injected into the bleed air stream for 90 min at a rate of 2-12 pounds 
per hour while the bleed air was sampled for carbon monoxide, carbon dioxide, 
nitrous oxides, and hydrocarbons. 





This trial was intended to characterize levels of oil-based contaminants in the bleed 
air stream through an intact engine compressor carbon seal but without a glass wool 
filter in the bleed stream. The tester injected oil for 15 min, during which time the 
same gaseous bleed air measurements were made. 





This trial was intended to simulate the potential for worst-case bleed air 
4 contamination by removing the engine compressor carbon seal. Of note, though, a 
glass wool filter was installed in the sampling line and no oil was injected. 





During this trial, the tester did not collect measurements but ran the engine for 





me 10 min and then shut it down to prepare for a “dirty start” in the following test. 
During this trial, the engine underwent a simulated “dirty start” which involved 
5b internal oil ingestion. The tester measured the concentrations of carbon monoxide, 


carbon dioxide, nitrous oxides, and hydrocarbons for the first 10 min after the 
engine was started. 


3.3.2. Additional Engine Bleed Air Testing 


During two of the oil injection tests (3c and 3f), the tester collected oil mist samples 
on “membrane filters” which were subsequently analyzed for tricresyl phosphates. The 
authors of the report noted that no “para or ortho isomers” of tricresyl phosphate were 
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present above the detection limit [16] (p. 24), although the authors did not define the 
detection limit or mention the concentrations of meta and mixed meta/para isomers which 
comprise more than 99% of commercial blends [21]. 

During tests 1-3, carbon sorbent tubes and impingers were used to collect samples 
that were subsequently analyzed by GC-MS [16]. The authors stated that these “analyses 
performed by the Environmental Protection Agency (EPA) [lab] ... did not lend themselves 
to meaningful interpretation because of apparent contamination of the samples and lack 
of parallel quantitation of known compounds... ” [16] (p. 27). Still, the authors affirmed 
that “based on the analyses, there was no significant toxicological gaseous content of the 
TPE331 engine bleed air... ” [16] (p. 27). The authors also referred to the presence of “a 
number of compounds in the bleed air under various conditions of the test protocol,” but 
described occupational exposure limits as protective for workers. The report noted that, 
by comparison, “the concentration of those substances [measured by the EPA] appears to 
be at such a low level as to cause no acute degradation of pilot performance” [16] (p. 66). 
Specifically, the authors stated that “it is apparent that toxic effects to aircraft crewmembers 
would only result from breathing air contamination of sufficient concentrations to cause 
acute effects” [16] (p. 66). 


3.4. Regulations 


Aircraft regulations stipulate that the cabin and flight deck ventilation supply systems 
must be designed to provide air that is “free from harmful or hazardous concentrations of 
gases and vapors” (14 CFR § 25.831(b)), which includes no more than 50 ppm of carbon 
monoxide and 5000 ppm of carbon dioxide (Table 4). Moreover, crews must be provided 
with “warning information” to alert them to unsafe conditions and enable them to take 
“appropriate corrective action” (14 CFR § 25.1309(c)) (Table 4). To date, though, these 
regulations have not been applied to monitor contaminants in the cabin air and alert pilots 
to the presence of oil fumes. 

In addition to the design regulations, U.S. airlines are required to comply with three 
FAA reporting regulations for fume events (Table 4), but all are underutilized [22]. 


4. Discussion 
4.1. Accident Reports 


In the investigated crashes (Table 1), eight of the ten aircraft were Mitsubishi MU-2 air- 
craft. The Mitsubishi MU-2 is a twin-engine turboprop aircraft with a pressurized cabin 
(Figure 8). The aircraft is known for its difficult handling and high rate of accidents, includ- 
ing fatal accidents [23]. Crash numbers 7 and 8 with this aircraft were in icing conditions 
which may be relevant, given an Australian Bureau of Air Safety Investigation report 
on two fatal crashes on MU-2 aircraft, both of which were also in icing conditions [24]. 
However, even for the 1988 and 1990 fatal crashes in Australia, “icing on the airframe” 
was only listed as the “probable” cause and, for one of those two crashes, pilot fatigue 
was also referenced extensively. Thus, the reference to icing conditions in crash 7 in this 
investigation should be acknowledged but does not rule out cabin air contamination, either 
as a causal or contributory factor. 

Regarding the remaining two crashed aircraft (Table 1), one was a Cessna 441 Conquest 
II and the other was an Aero Commander 690B. As with the MU-2, each of these is a 
turboprop aircraft with a pressurized cabin. 

In at least four of these ten crashes, there are details which are consistent with impairment 
caused by oil fumes (Table 5). For each of the crashes, the NTSB identified “probable causes” 
(including “undetermined”), all of which are reproduced in Table 2. In many cases, though, 
the “probable cause” is more of a description than an explanation. The real question is, why? 
Why did the pilot in command “misjudge distance and altitude” or “fail to obtain/maintain 
speed”? Why was air traffic control “unable to get response from crew”? Was it because these 
pilots were distracted or fatigued, or were they impaired by fumes? 
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Figure 8. Mitsubishi MU-2B (Photograph by Alan Lebeda [25], trimmed, GFDL 1.2). 


It is tempting to blame crashes on “pilot error,” typically attributed to either to inexperi- 
ence or fatigue. Regarding inexperience, the flight time for the pilots on the crashed planes 
ranged from 1498 to 16,766 h (average 7999) (Table 2). One analysis estimated that, during the 
time of these crashes, turboprop pilots logged an average flight time of 528 h per year [26], so 
if the pilots on the crashed planes reflect this average, they are not novices. Regarding fatigue, 
the accident reports did not consistently cite either each pilot’s total flight time specific to 
the aircraft type they crashed or the number of hours they flew during the 24 h prior to the 
accident flights. As a result, the impact of pilot fatigue cannot be assessed. 

The authors of this paper are unaware of other challenges to the basis for the claims in 
the NTSB report [16] that it is not possible for oil fumes to either have contaminated the 
bleed air on these aircraft or to have played any role in any of these accidents. 

Theirs is not the first report, however, to either overlook or downplay the flight 
safety implications of breathing oil or hydraulic fluid fumes. For example, one published 
case study investigated the findings of six documented fume event investigations and 
identified specific oversights and misconceptions about the potential for oil to contaminate 
the bleed air, as well as a tendency to dismiss the crew-reported symptoms as stress 
reactions [27]. Another investigation described an event during which pilots reported a 
blue haze shortly after reaching cruise altitude, prompting a precautionary landing [28]. 
The initial ground-based engine runs failed to identify any fault. Only additional (and 
non-standard) high-power engine runs identified the source of fumes—a fractured bearing 
seal in the compressor. This is important because high-power runs were not part of the 
fault-finding procedure endorsed by the aircraft manufacturer. 

Another investigation described how both pilots reported fumes during approach and 
felt dizzy and nauseous, with irritated eyes and throat [29]. The pilots donned oxygen and 
requested priority landing clearance. The auxiliary power unit (APU) was not operating 
during the approach phase, so it was not investigated as a potential source of oil contami- 
nation. However, the aircraft manufacturer has published bulletins alerting airlines that, 
when the APU is in use, oil fumes can contaminate the downstream environmental control 
system and ducting, and fumes from that initial contamination can continue to manifest 
inflight, even when the APU is not in use [30]. Still, the APU was not checked, and the 
investigation was inconclusive; “no explanation” for either the fumes or the symptoms 
was found. 

In 2020, the French Bureau of Enquiry and Analysis for Civil Aviation Safety (BEA) 
investigated a serious incident on a commercial flight during which the pilots reported a 
strong, acrid odor accompanied by irritant symptoms during taxi, shortly after the APU 
had been turned on [31]. During the climb out, both pilots reported dizziness, among other 
symptoms. They eventually donned oxygen masks and diverted to the nearest airport. 
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Although the details are consistent with breathing bleed-sourced oil fumes, maintenance 
inspections did not identify an obvious mechanical failure. As a result, the accident 
investigator concluded that the incident—including fumes, persistent symptoms, and 
a diversion—was most likely caused by the pilots inhaling “an excessive quantity of 
carbon monoxide.” Failing to find a more plausible explanation, the accident investigators 
hypothesized that the source of the fumes may have been engine exhaust from a small 
business jet as it taxied about 80 m away from the incident aircraft in the presence of a 
“calm wind.” However, the report added that it “cannot be excluded that the crews were 
intoxicated by another substance.” This incident highlights the value of installing and 
operating sensors onboard, both to provide real-time information to the pilots and assist 
maintenance in subsequent troubleshooting. 

Other than its investigation into these ten turboprop crashes in 1979-1981, the NTSB 
has not weighed in on the flight safety implications of pilots inhaling oil-contaminated 
bleed air inflight. This is contrary to its counterparts in Australia [32,33], Germany [4], 
Iceland [34], Ireland [35], New Zealand [36], Spain [5], Sweden [37], Switzerland [7], and 
the UK [3,6], for example. Each of those safety boards has investigated one or more 
commercial flights involving pilots whose inflight impairment compromised flight safety 
and has concluded that breathing contaminated bleed air either likely or definitively played 
a role. The NTSB does require airlines to report onboard fume/smoke events, but only if 
either the aircraft sustained structural damage or if specific health impacts were reported 
by one or more occupant (49 CFR § 830.2). However, oil fume events that meet the NTSB 
reporting criteria have not been investigated. 


4.2. Bearing and Seal Failures in the TPE331 Engine 


The description in [16] of the two seals in series contains some contradictions. Accord- 
ing to the explanation, the carbon seal is necessary because “when the engine is shutdown, 
oil drains from the engine walls which increases the oil level above the shaft installation 
level of the compressor” [16] (p. 33). The labyrinth seal is necessary because the carbon 
seal is potentially “subject to wear and damage” [16] (p. 2). So, if the carbon seal is worn or 
damaged, then the oil will flow out of the reduction gearbox into the compressor and con- 
taminate the cabin air when the engine is next started up. This situation was investigated 
in test 5 [16]. 

It is known that all seals leak in small quantities [38]. Pressure differentials assumed by 
the manufacturer may not be present in failure cases. Moreover, pressure differentials may 
not be sufficient during rapid thrust reduction or at idle thrust. In test 4, it was assumed 
only that one of the two seals in series (the carbon seal) would have failed. The labyrinth 
seal was left intact. Certainly, a mechanical failure that causes the carbon seal to fail may 
also cause the labyrinth seal to fail at the same time. For example, a failure of the main 
shaft compressor bearing could cause the main shaft to vibrate or rotate eccentrically which 
could severely damage (or even destroy) both seals. 

In addition to crash number 1 (Table 5), an TPE331 compressor bearing failure which 
caused the engine to fail was documented by the Australian Transport Safety Bureau 
(ATSB) [39], although without either fatalities or reference to oil fumes. That report 
described the history of compressor bearing failures on the TPE331 engine which peaked 
in the early to mid-1980s. Garrett attributed the compressor bearing failures to propeller 
strikes and subsequently revised the engine maintenance manual, which—by the 1990s— 
reduced the number of failures. In the accident investigated by the ATSB, when the 
compressor bearing failed, the bearing cage fractured which caused the balls to come loose. 
Such loose components can collide with and cause damage to numerous internal rotating 
engine components, including the compressor seals. 


4.3, Bleed Air Testing 


Regarding the engine bleed air testing reported by the NTSB [16], the conditions in 
trials 1-4 did not represent the conditions of concern on the crashed aircraft because the 
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engine compressor carbon seal was intact and, in all but one case, the air was filtered for 
vapor and particulate. No measurements were collected during trial 5a. Regarding the 
relevance of conditions tested during trial 5b, on the one hand, the bleed air would have 
contained oil from the dirty engine start, but it is not clear if this is the same volume and 
pattern of oil as with a faulty engine compressor carbon seal inflight. Moreover, the bleed 
air had been passed through a glass wool filter for vapors and particulate, unlike on the 
crashed aircraft. Finally, it is not clear if the temperatures in the test stand compressor 
reflected inflight conditions. 

It is important to carefully consider the question of the temperature in the engine 
compressor on these crashed aircraft. In its report, the NTSB referenced some certification 
data collected by a major engine manufacturer after injecting oil into a large turbofan engine. 
The tests showed that that carbon monoxide formation was detectable starting around 
316 °C [16] (p. 7). The report also notes that TPE331 engine has a standard compressor 
discharge temperature of 360 °C [16] (p. 6) with a maximum temperature of 386 °C [16] 
(p. 25). The compressor temperature was high enough for carbon monoxide to be generated 
in the presence of oil, but the authors claim that the concentration of carbon monoxide 
would have been too low to cause impairment [16] (p. 2). The authors also rule out the 
presence of other oil decomposition products (such as acrolein) sufficient to irritate the 
mucous membranes by noting that “a number of participants sniffed the bleed air lines ... 
[and] no one described an acrid or irritating quality”; [rather] “the consensus was that the 
odor was that of a warm oil, not a decomposed oil” [16] (p. 27). 

In a companion study, CAMI researchers exposed test animals to Exxon 2380 oil 
for seven hours without, what the NTSB called, “any immediate or delayed behavioral 
change” [16] (p. 28). The CAMI report does describe an experiment with rats and chickens 
exposed to aerosolized (but unheated oil) in which the animals seemed unaffected [17] 
(p. 14). However, in the same study the CAMI researchers reported that, when the engine 
oil was heated, carbon monoxide started to form at 306 °C and when rats in one trial were 
exposed to oil fumes heated to 350 °C for 30 min, “it was obvious that the animals were 
approaching incapacitation” [17] (p. 8). Moreover, when the temperature was further 
increased, the animals “expired.” Unlike the rats, the pilots on the crashed planes would 
have inhaled oil fumes in a reduced pressure environment with a corresponding reduction 
in the partial pressure of oxygen in the ambient air. 

More recent sampling data has highlighted the presence and potential health impact 
of elevated levels of ultrafine particles in the bleed air supplied to the cabin and flight 
deck on commercial aircraft [40-44]. As noted above, the authors of the NTSB report 
assumed that “particulate matter would not form” [16] (p. 9), and so did not assess any 
associated toxicity. 


4.4, Regulations 


In the 1960s, a report from an aircraft manufacturer (Douglas Aircraft Corporation) 
describes the need to “show that the level of contamination required for olfactory warning 
is well below the generally accepted tolerance limits for the toxic materials produced” [15] 
(p. 3). In 1960, U.S. aircraft manufacturers would have been required to demonstrate to the 
U.S. Civil Aeronautics Board (which predated the FAA) that aircraft systems were designed 
to provide ventilation air “free from harmful or hazardous concentrations of gases or 
vapors”, which included (but was not limited to) the carbon monoxide concentration not 
exceeding 50 ppm [45]. Presumably, this is what motivated Douglas Aircraft Corporation 
to propose testing to demonstrate that, “under all possible oil leakage rates,” olfactory 
warnings would give pilots sufficient time to “permit effective crew countermeasures.” [15] 
(p. 3). Essentially, tests needed to show that pilots would smell oil fumes before they would 
be impaired, proving that the presence of oil-based contaminants in bleed air was safe. 

Since then, many studies and reports discuss the presence of airborne contaminants 
in cabin air and comment on whether the levels are “safe” or “acceptable” [40,46—48]. 
“Acceptability” has more commonly been framed in the context of comparing aircraft data 
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to published exposure limits for individual chemical constituents [16,40,49,50], although 
this approach is not without criticism [51,52]. 

The quest to define “safe” or “acceptable” concentrations of oil-based compounds 
is likely borne from the need for industry to comply with the regulation that the aircraft 
supply system be designed to supply air that is “free from harmful or hazardous concen- 
trations of gases or vapors” (14 CFR § 25.831(b)). If something is acceptable—if it is not 
harmful—then the design meets the regulation. Interestingly, though, the FAA does not de- 
fine “harmful or hazardous” beyond its limits for carbon monoxide (14 CFR § 25.831(b)(1)) 
and carbon dioxide (14 CFR § 25.831(b)(2)) and it had not done so in the early 1980s either. 
Applying published exposure limits for individual chemicals to complex mixtures [53] 
(Appendix 10) is problematic [54], especially in a reduced pressure and safety-sensitive en- 
vironment [52]. This is in part because exposure limits do not exist for every constituent in 
these mixtures and the various constituents have a variety of toxicity endpoints. Moreover, 
the diversity of exposure limits published by different authorities for a single chemical 
illustrates the fallacy that a single number can define the boundary between what is and is 
not acceptable. In the case of the FAA, their limits for carbon monoxide and carbon monox- 
ide are as high—or higher—than published chemical exposure limits for ground-based 
application. So, given this gap in the regulatory framework, in 1981, the NTSB and their 
industry colleagues defined “safety” for themselves—not in the context of certification 
regulations, but in the context of whether pilots could have been impaired. 

The authors of the 1984 NTSB report acknowledged that carbon monoxide would 
be generated when oil was purposefully introduced into the compressor section of the 
TPE331 engine, but they noted that it should be “well below permissible limits in the 
bleed air” [16] (p. 2). However, the permissible exposure limit for carbon monoxide is 
not a useful benchmark for safety. An applied research study into the effects of exposure 
to carbon monoxide inflight concluded that “the maximum permissible concentration of 
carbon monoxide in pressurized passenger airplane cabins should be 0.01 mg/liter,” which 
is equivalent to 8.7 ppm [55]. By comparison, the permissible aircraft design limit for 
carbon monoxide is 50 ppm (14 CFR § 25.831(b)(1)) which is, inexplicably, as high, or higher 
than ground-based limits. This, even though it is primarily applied in a reduced pressure 
environment (14 CFR § 25.841(a)) with a corresponding reduction in the partial pressure of 
oxygen available to occupants. Safety of flight depends on pilots’ alertness and reaction 
time which can be compromised by exposure to asphyxiants, such as carbon monoxide [56]. 
Moreover, aircraft occupants have no means of egress once the doors are closed, such that 
an additional safety factor regarding exposure to carbon monoxide is warranted. Finally, 
carbon monoxide is only one element of a complex chemical mixture. 


4.5. Exposure Control Measures 


Just as the flight safety implications of breathing oil and hydraulic fluid fumes have 
been well documented, the need for control measures has also long been recognized. In 
1966, an aircraft manufacturer reported that, in response to the evidence that contaminated 
bleed air could have caused some fatal crashes, the Navy required that crewmembers 
breathe 100% oxygen from takeoff to landing [15]. In 1955, an engineer with North Ameri- 
can Aviation described the outcome of a two-year investigation into engineering solutions 
to prevent bleed air contamination, recommending either a separate compressor or a bleed 
air filter [57]. In response to one oil fume event during the descent phase of a commercial 
flight, the UK Air Accidents Investigation Board recommended that the FAA and EASA 
“consider requiring” a system to warn pilots about the presence of “smoke or oil mist in the 
air delivered from each air conditioning unit” [538]. Engineering and operational measures 
continue to be called for, including bleed air filtration, sensors to provide early warning of 
airborne contaminants, improved maintenance procedures, more targeted troubleshooting 
procedures, relocating the air inlet for the auxiliary power unit, and airline worker training 
and education [59,60]. 
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5. Conclusions 


Between 1979 and 1981, NTSB investigators suggested that inhaling oil fumes inflight 
may have been a causal factor in “several” fatal crashes of turboprop aircraft. In response, 
the agency partnered with the very companies that had a commercial interest in the 
outcome of the investigation. Aside from noting the potential for some people to develop 
“extreme chemical sensitivity” [16] (p. 28) to chemicals in oil fumes, the NTSB and their 
industry partners soundly dismissed the hypothesis that oil fumes may have impaired some 
of the pilots on the crashed planes, affirming that it was “completely without validity” [16] 
(p. 3). However, their conclusions are not supported by the evidence for three key reasons. 

First, it was not—and is not—possible to draw any definitive conclusions regarding 
the cause of these crashes because real-time bleed air testing and suitable post-mortem 
blood analyses were not available. However, the accident reports for four of the fatal 
crashes on these turboprop aircraft include details consistent with (and suggestive of) the 
pilots being impaired by oil fumes (Table 5). Moreover, this is consistent with similar 
concerns and reports of pilot impairment documented since the 1930s, such that “pilot 
error” and “undetermined” are insufficient explanations. Without air supply monitoring 
equipment to provide real-time warning, and without options for a blood test to investigate 
inhalation of oil fumes, aircraft crashes that are attributable to bleed air contaminants will 
not be recognized as such. 

Second, none of the tests of ground-based bleed air measurements of a subset of 
oil-based contaminants generated in the engine type that had been on the crashed aircraft 
reproduced the types of inflight conditions that the accident investigators had flagged as 
potentially unsafe (Tables 3 and 5). Specifically, the engine test stand conditions did not 
assess the impact of oil seeping through a cracked or otherwise damaged compressor seal 
on the quality of unfiltered bleed air downstream. As a result, the bleed air testing results 
are not relevant to the question of whether oil fumes could have impaired the pilots inflight. 
As such, the argument that they somehow discount the potential for impairment is invalid. 

Finally, to interpret the bleed air data they collected, the authors of the report used pub- 
lished exposure limits as a benchmark for whether the concentration of gaseous compounds 
in oil fumes would have been sufficient to cause pilot impairment. However, published ex- 
posure limits are not appropriate for assessing the risks associated with inhaling a complex 
mixture of compounds in an enclosed and reduced pressure environment. 

The significant concerns raised by these crashes and the history of pilot impairment 
associated with breathing oil fumes all support more stringent design, operation, and 
reporting regulations to protect safety of flight. 


6. Recommendations 


e For current aircraft that are equipped with bleed air systems, engineering control 
measures such as sensors and filters should be mandated to prevent inflight exposure 
to fumes and, thus, improve flight safety. For new aircraft types, non-bleed air supply 
systems should be standard. 

e Given the flight safety implications, all crewmembers should be trained to recognize 
and respond to the presence of bleed air contaminants [61]. 

e The NTSB should reopen the 1984 report and update the conclusions based on what 
was known then about oil fumes and flight safety, and considering the data and reports 
published since then. The NTSB should also issue recommendations to the FAA on 
necessary actions to minimize the flight safety impacts of exposure to oil fumes. 

e Until suitably protective measures are implemented fleet-wide, the NTSB should consider 
pilot impairment from bleed-sourced fumes as a potential casual factor in future aircraft 
accident investigations. This is particularly important when an accident includes either 
relevant maintenance history or a pilot’s failure to communicate, for example. 

e If only to avoid the appearance of bias, future investigations should be independent 
of all commercial interests. Even though it is not unusual for the NTSB to include the 
FAA and manufacturers in accident investigations, both the poorly conceived design 
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of the air sampling trials and the sweeping conclusions, suggest that commercial 
conflicts influenced this project. 
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Abbreviations 


The following abbreviations are used in this manuscript. 


CAMI United States Civil Aeromedical Institute 

CFR United States Code of Federal Regulations 

Ge) Certification Specification (issued by EASA) 

EASA — European Union Aviation Safety Agency 

FAA United States Federal Aviation Administration (issued by the FAA) 
FAR United States Federal Aviation Regulation 

lb/hr pounds per hour 

NTSB United States National Transportation Safety Board 


psia pound-force per square inch absolute (relative to zero or absolute vacuum) 
psig pound-force per square inch gauge (relative to atmospheric pressure; at sea level 14.7psi) 
rpm revolution per minute 


TPE turboprop engine 
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Abstract: For aircraft employing the fly-by-wire technique, the closed-loop dynamic characteristics 
are determined by both the configuration design and the flight control system. As the capacity 
of the control system has certain limitations, the configuration parameters are also constrained by 
the requirements of the closed-loop flying qualities. This paper presents an aircraft configuration 
parameter boundaries determination method based on closed-loop flying qualities requirements 
independent of the actual flight control law design, mainly aiming at the parameters that affect 
the stability and control characteristics. First, a nonlinear dynamic inversion-based flight control 
law is adopted to decouple the control law gains from the configuration parameters and to study 
the relationship between the configuration parameters and closed-loop flying qualities. Second, a 
flying qualities evaluation scheme is established by selecting the most severe flight conditions and 
the evaluation criteria that are most sensitive to changes in the parameters. Finally, the parameter 
boundaries according to the requirements of Level 1 flying qualities are determined by searching for 
the critical values that lead to degradation of the flying qualities. The proposed method is verified by 
an application example of the design ranges of a sample aircraft’s wing position, horizontal tail area, 
center of gravity, vertical tail area and vertical tail position. 


Keywords: configuration design; flying qualities; flight control law 


1. Introduction 


Aircraft configuration design depends on the design requirements of different types of 
parameters [1,2]. For example, the determination of wing area, fuselage geometry and cross 
section mainly considers the payload requirement; the airfoil, aspect ratio and sweep angle 
of the wing are designed according to the aerodynamic characteristics requirements; wing 
position, horizontal tail area, center of gravity (CG), vertical tail area and position have 
considerable influences on the dynamic characteristics of aircraft, and the design of these 
parameters is conducted mainly in consideration of stability and control requirements [3]. 

The flying qualities evaluation can comprehensively assess the stability and control 
characteristics, the performance of the flight control system, and the adaptability to the 
flight missions [4—8]. Thus, the configuration design should take into account the flying 
qualities requirements. As the developments of digital flight control systems and the 
fly-by-wire technique allow more advanced flight control laws (FCL) to be implemented 
on aircraft [9-15], the configuration design has become more diversified. However, the 
capacity of FCL has certain limitations, and flying qualities cannot be guaranteed only rely- 
ing on the control system. Considering that the flying qualities of the closed-loop aircraft 
are determined by both the configuration design and the flight control system, the configu- 
ration parameters are also constrained by the closed-loop flying qualities requirements [16]. 
Therefore, it is necessary to establish the relationship between the configuration parameters 
and the closed-loop flying qualities requirements in the conceptual aircraft design. 


https://www.mdpi.com/journal/aerospace 
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Recently, there has been renewed interest in the closed-loop flying qualities prediction 
and assessment in conceptual aircraft design [17-22]. The Institute of Aeronautics and 
Astronautics of RWTH Aachen University developed a methodology for the prediction 
and assessment of flying qualities and implemented it into a preliminary aircraft design 
tool [20]. Arthur Rizzi et al. achieved the multidisciplinary optimization of a canard con- 
figured TransCruiser with flight control by using CEASIOM, a multidisciplinary software 
environment for aircraft design [21]. Stefano Cacciola et al. realized the optimization of the 
longitudinal control law, as well as the three-surface configuration for better flying quali- 
ties [22]; however, the FCLs involved in these studies were mainly stability augmentation 
systems (SAS) or stability and control augmentation systems (SCAS), which are based on 
linear control theory and well-established gain-scheduling methods [16]. As a result, the 
controller gains are coupled with the configuration parameters of the aircraft, and thus, 
it is impossible to identify the influence of the configuration parameters separately from 
the influence of the controller gains on the closed-loop flying qualities. Due to this lack 
of flexibility, the traditional FCLs are not suitable to study the relationship between the 
configuration parameters and the closed-loop flying qualities requirements. 

Inverse model-based control techniques, e.g., nonlinear dynamic inversion (NDI) 
as well as incremental nonlinear dynamic inversion (INDI), offer a great advantage to 
continuously changing plants [23,24]. The studies of DLR (German Aerospace Center) 
and others [25,26] show that the flight control law based on the NDI method can effec- 
tively adapt to changing aircraft configurations. Therefore, this paper proposes an aircraft 
configuration parameter boundaries design method based on closed-loop flying qualities 
requirements. This method realizes the configuration parameter boundaries design inde- 
pendent of actual flight control law design by using the NDI method, and the closed-loop 
flying qualities requirements are transformed into the design ranges of target configu- 
ration parameters. The parameter boundaries simplify the stability and control design, 
so that the configuration design can be more focused on aerodynamic, structural and 
other requirements. 

The remainder of this paper is structured as follows. Section 2 introduces the de- 
termination method of configuration parameter boundaries based on flying qualities 
requirements. Section 3 defines the wing position, horizontal tail area, CG position, vertical 
tail area and vertical tail position as the target configuration parameters. Section 4 describes 
the modelling of the closed-loop aircraft and proposes the selection of flight conditions and 
evaluation criteria. An application example comprising the target parameter design of a 
medium-range transport aircraft is then shown in Section 5. Finally, a concluding summary 
is given in Section 6. 


2. Configuration Parameter Boundaries Based on Flying Qualities Requirements 


Conceptual design is the first phase of the aircraft design process. It involves sketching 
a variety of possible configurations that meet the design requirements, such as aerodynam- 
ics, propulsion, flight performance, structural and control systems. 

The steps of conceptual design are described as [1,2,27]: 


e Preliminary Sizing: definition of requirements, selection of an aircraft configuration, 
selection of a propulsion system, etc. 
Layout sizing: fuselage sizing, wing sizing, empennage sizing, design for high lift, etc. 
Specific design: stability and control design, flight control system design, landing gear 
design, structures design, etc. 

e Aircraft performance analysis: aerodynamic analysis, polar estimation, mass estima- 
tion, flying qualities evaluation, mission analysis, etc. 


The whole design process has an iterative nature to realize optimization design [27]. 

Stability and control design and flying qualities evaluation are often in later stages 
of the design phase. Configuration parameters adjustment due to stability and control 
requirements or flying qualities requirements may lead to much iterative work. This paper 
proposes a method to determine the configuration parameter boundaries based on closed- 
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loop flying qualities requirements to simplify the stability and control design. The proposed 
parameter boundaries aim to transform the closed-loop flying qualities requirements into 
the design ranges of target configuration parameters. In this way, the flying qualities 
requirements can be considered at the layout sizing stage during optimization design. 

The determination process of the parameter boundaries is as follows. The first step is 
to construct a high-order closed-loop aircraft model according to the basic configuration 
design; the model can be divided into an aircraft motion model and the FCL based on the 
model reference NDI method. The second step is to establish the flying qualities evaluation 
scheme by selecting the most stringent flight conditions and the evaluation criteria that 
are most sensitive to the parameter changes. The third step is to establish a reasonable 
variation series of the target parameters and to search for the critical values that lead to 
degradation of the flying qualities. Finally, the parameter boundaries according to the 
requirements of the Level 1 flying qualities (Level 1 FQ) are determined. The method for 
determining parameter boundaries is shown in Figure 1. 


evaluation criteria 


Flying qualities 
level down ? 


Yes 
Boundaries of parameters 


Figure 1. Configuration parameter boundaries based on closed-loop flying qualities requirements. 








It should be emphasized that the parameter boundaries proposed in this paper cor- 
respond to specific flying qualities criteria, so they only reflect a part of stability and 
control requirements. Meanwhile, these configuration parameters should also meet other 
requirements. For example, the CG position design should meet the trim requirements, the 
vertical tail design should meet the controllability requirements in the single-engine out 
condition at take-off, and so forth. Designers should seek to reach the design configuration 
that satisfactorily meets all requirements in the optimization design process. 


3. Configuration Parameter Analysis 
3.1. Wing Position and Horizontal Tail Area 


The longitudinal stability and control characteristics of aircraft are determined mainly 
by the stability and damping derivatives of the pitch moment, and these derivatives are 
affected by the wing position, CG position, horizontal tail position, horizontal tail area, and 
other parameters. The horizontal tail is usually installed at the rear end of the fuselage to 
ensure the control power of the elevator, and the variation range of the tail position is small. 
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Thus, the longitudinal characteristics are governed predominantly by the wing position 
and horizontal tail area [1]. These two parameters are usually designed in collaboration 
according to the requirements of the longitudinal flying qualities. 

The stability derivative of the pitch moment Cy, can be expressed as [28]: 


Cra = Cra.wB (Xeg. — Xac.wB) — Cra.tVr + Cina.p (1) 


where Xac.wp is the mean aerodynamic center of the wing-body configuration and is 
affected mainly by the wing position; Cr, we is the lift coefficient derivative of the wing- 
body configuration; X¢.¢. is the CG position of the aircraft; Cy,.7 is the lift coefficient 
derivative of the horizontal tail; Vr is the horizontal tail volume ratio, affected mainly by 
the horizontal tail position and its area; and Cinq» is the pitch moment derivative of the 
propulsion system. Vr is defined as: 


IpSt 

are (2) 
where /r is the distance between the aircraft CG and the aerodynamic center of the hori- 
zontal tail, Sy is the horizontal tail area, c is the mean aerodynamic chord of the aircraft 
wing, and S is the wing area. Consequently, forward movement of the wing position and a 
decrease in the horizontal tail area lead to a decrease in |Cinq| and a reduction in the static 
longitudinal stability. Similarly, rearward movement of the wing position and an increase 
in the horizontal tail area lead to an increase in |C,,| and an increase in the demand on the 
pitch control moment. 

The damping derivatives of the pitch moment, namely, Cing and C,,,,, are generated 
mainly by the horizontal tail, whereas the influence of the wing-body configuration ac- 
counts for only about 10% of the overall influence [28]. The contributions of the tail to Cig 
and C_,, are [28]: 


Cing.T = —2Cia.t Vr C 3) 


_ Ir de 
Cini.T _ —2Cie.t V0 3g 


where de/0d« is the derivative of the downwash angle of the horizontal tail with respect to 
the angle of attack. A decrease in horizontal tail area leads to a decrease in both |Ciq| and 
Ca . 

The stability and damping derivatives of the pitch moment mainly affect the longitudi- 
nal short-period mode characteristics of aircraft. The natural frequency Wpsp and damping 
ratiO Csp of the short-period mode can be expressed as: 


z — My +My — Zo (4) 
is 2Wn sp 


where Mx, Mg, M,, and Z, are the dimensional dynamic derivatives related to Cinw, Cing, 
C., and Cra, respectively. The detailed derivation process from (3) to (4) is described in 
the references [29,30]. 

Combining Equations (3) and (4) demonstrates that a decrease in the horizontal tail 
area leads to a decrease in ¢sp, slows down the short-period convergence, and decreases the 
dynamic stability of the aircraft; in contrast, an increase in the tail area leads to an increase 
in the damping ratio and restrains the pitch response to the pilot control input. 

In a general conceptual design process, the horizontal tail is first of all sized from static 
conditions as done with the “V-diagram” [27,31,32], which corresponds to the basis of the 
longitudinal stability and control requirements. Then, the horizontal tail design should 
consider the dynamic stability and control requirements, which are the primary focus of 
this paper. 
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According to the above analysis, the longitudinal dynamic stability requirements of 
the aircraft should be considered to determine the forward range of the wing position and 
the maximum limit of the horizontal tail area. Furthermore, the rearward range of the wing 
position and the minimum limit of the horizontal tail area should be determined by the 
longitudinal maneuverability requirements. 


3.2. Center of Gravity Position 


The CG is not a fixed point but moves with the initial loading condition and fuel 
consumption of the aircraft. For an aircraft with determined aerodynamic shape, a change 
in the CG position x¢g causes Cj, to vary and thus affects the longitudinal characteristics. 
Therefore, the range of the CG location needs to be strictly designed according to the 
requirements of the longitudinal flying qualities. 

A forward CG represents a large |Ci,| and requires a larger pitch control moment, 
while a rearward CG position represents a small |Cin.| and poor longitudinal stability. 
Therefore, in this paper, the forward CG limit mainly considers the longitudinal maneu- 
verability requirements, and the rearward CG limit is determined by the longitudinal 
stability requirements. 

Additionally, the forward CG limit should allow the aircraft to be trimmed at maximum 
lift coefficient in ground effect. The final forward CG limit will be determined by both the 
longitudinal maneuverability and trim requirements, which leads to future investigations. 


3.3. Vertical Tail Area and Position 


Directional stability, also known as weathercock stability, is provided primarily by the 
vertical tail. Meanwhile, a large vertical tail is accessible spiral stability. However, a large 
vertical tail surface may produce an excessive aerodynamic drag force during high-speed 
flight at high altitudes in cruise, and such a large surface would require a higher structural 
strength and mass of the airframe. Therefore, it is necessary to restrict the vertical tail 
area to optimize both the aerodynamic characteristics in cruise and the airframe structure. 
To ensure the directional stability after the vertical tail area is reduced, it is necessary to 
increase the vertical tail arm by adjusting the position of the vertical tail. Similar to the 
horizontal tail, the vertical tail is also installed at the end of the fuselage to ensure the 
control power of the rudder. Thus, the vertical tail position should be adjusted by changing 
the chord length, sweep angle and other parameters without moving the root position. It 
should be noted that a large sweep angle may lead to the decrease of the vertical tail lift 
coefficient. Thus, the vertical tail design should avoid a large sweep angle [29,30]. 

The contributions of the vertical tail to the stability and damping derivatives of the 
directional moment C,,.y and Cy,,y can be expressed as [28]: 


Cugp.v = —Cyp.vW 
Curvy = ~2Cygy Vv ¥ (5) 


ly S 
YW oo 


where Cy,,v is the lift coefficient derivative of the vertical tail; Vy is the vertical tail volume 
ratio; Sy is the vertical tail area; ly is the force arm of the vertical tail; and b is the wingspan 
of the aircraft. 

The stability and damping derivatives of the directional moment affect mainly the 
Dutch roll mode characteristics of the aircraft. The natural frequency wy 4, and damping 
ratio Cg, of the Dutch roll mode can be expressed as: 


Whar = 4 / Np = Np Yr = N,Yp 
— 2Wn dr 
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where N pr Nr, Y, and Y; are the dimensional dynamic derivatives related to Cup, Cus Cyg 
and Cy,, respectively. The detailed derivation process from (5) to (6) is also described in 
the references [29,30]. Combining Equations (5) and (6) reveals that a reduced vertical tail 
arm and a decrease in the tail area cause a reduction in both iC B| and ¢g, and decreases 
the directional static and dynamic stability. 

In a general conceptual design process, vertical tail sizing for static stability has to 
compensate for the destabilizing effect of the fuselage; vertical tail sizing for static control 
has to compensate the yaw moment from a one-engine out flight case [29,30]; vertical tail 
sizing for roll response control requires the proper balance of the rolling time constant 
and the roll sensitivity; there are also spiral stability and directional dynamic stability 
requirements. Directional dynamic stability mainly refers to the Dutch roll mode stability, 
which is the primary focus of this paper. The eventual vertical tail area will be the minimum 
area that meets all these requirements. 


4. Aircraft Modeling, Selection of Flight Conditions and Evaluation Criteria 
4.1. Aircraft Modelling 


The aircraft motion model is based on the six-degree-of-freedom equations of motion. 
To determine the configuration parameter boundaries based on closed-loop flying qualities 
requirements, the motion model should be constantly updated according to the changes in 
the configuration parameters. 

To easily adapt the control laws to continuously change the aircraft configuration 
during the design process, NDI controllers offer considerable advantages. For instance, the 
model reference NDI control law [33,34] uses transfer function-based reference models to 
describe the desired aircraft dynamics. The reference models can be expressed as [33]: 


qr | KqrWpr(6 - 1/Te21) 
Ac s* 26sp.rWsp.rS le Cie 


Pr _ pr (7) 
Pc S+ Wye 
Ty Wrz 


‘c S+ Wry 


where Pc, Jc, and rc are the rate commands of the pilot; p;, g,, and r; are the corresponding 
reference rate commands generated by the reference models; K,, and K,.; are the control 
gains; Csp.r, Wsp.r, and To2, are the desired damping ratio, natural frequency and time 
constant of the short-period mode, respectively; and wy, and w;.; are the desired roll and 
directional frequencies. The reference commands are achieved by the actual dynamic 
inversion at the core of the control law, which computes the necessary surface positions. 
A proportional-integral (PI) compensator is added to drive down the error between the 
desired dynamics and the actual dynamics. Filters are also needed to attenuate undesirable 
noise from the feedback sensors. The structure of the high-order closed-loop aircraft model 
with the model reference NDI control law is shown in Figure 2. 

The parameters of the model reference NDI control law are designed according to the 
desired dynamic characteristics and do not need to be adjusted according to changes in the 
configuration parameters [25]. As a result, the FCL achieves the desired closed-loop flying 
qualities of the varying aircraft configuration and eliminates the disturbance attributable to 
variation in the control system. Therefore, the model reference NDI method is employed 
to design the control law in this paper, which is convenient for studying the relationship 
between the configuration parameters and the closed-loop flying qualities. 
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Figure 2. Structure of the high-order closed-loop aircraft model. 
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It should be emphasized that the design ranges of the configuration parameters 
derived from the NDI control law represent the capability of control laws that makes the 
closed-loop aircraft meet the flight quality requirements. For aircraft adopting other control 
laws, such as SCAS, the conclusions are still valid. 


4.2. Selection of Flight Conditions for Evaluation 


Any variation in the flight conditions will lead to a change in the flying qualities 
evaluation results, so the requirements for the configuration parameters differ under 
different flight conditions. Consequently, to determine configuration parameter boundaries 
based on closed-loop flying qualities requirements, it is necessary to determine the flight 
conditions with the most severe requirements for the configuration parameters [1,3]. 

The previous discussion shows that the forward range of the wing position, the 
minimum limit of the horizontal tail area, and the rearward CG limit are the main factors 
influencing the longitudinal dynamic stability requirements of the aircraft. In contrast, 
the longitudinal maneuverability requirements are governed primarily by the rearward 
range of the wing position, the maximum limit of the horizontal tail area and the forward 
CG limit. For aircraft whose load and fuel are arranged mainly in the front region of the 
airframe, the forward CG limit appears during the take-off phase [1]. In this case, the 
stability margin of the aircraft is large during the take-off phase, and the longitudinal 
maneuverability is poor. Considering the low flight speed during the take-off phase, the 
aerodynamic efficiency of the control surface is low, which worsens the maneuverability. 
Therefore, the rearward range of the wing position, the maximum limit of the horizontal 
tail area and the forward CG limit should be determined in the take-off phase. The flight 
speed, flight altitude, weight distribution, flaps, landing gear and other settings should 
be consistent with the most forward CG position. With loading and fuel consumption, 
the CG position reaches the rearward limit in landing phase, during which the stability 
margin of the aircraft decreases and the most stringent requirements for the longitudinal 
stability are generated. Therefore, the forward range of the wing position, the minimum 
limit of the horizontal tail area and the rearward CG limit should be determined in landing 
phase, and the exact flight condition should be set to the state with the most rearward 
CG position. The selected flight conditions for evaluating an aircraft whose load and fuel 
are arranged mainly at the rear of the airframe are the opposite. For an aircraft whose 
load and fuel cannot be judged clearly whether they are arranged in the front or rear 
position (for example, the fuel is in the wings and the wingbox), the flight conditions for 
longitudinal flying qualities evaluation can be selected directly according to the change of 
the CG position during flight. 

The vertical tail design for the directional dynamic stability should be determined 
under the flight condition with the most stringent requirements. Due to the decrease in 
damping coefficients during high-speed and high-altitude flight, the Dutch roll mode 
of the aircraft will deteriorate, thereby increasing the dynamic stability requirements. 
Therefore, the cruise phase at high speed and high altitude is selected to determine the 
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design ranges of the vertical tail area and position based on the directional dynamic 
stability requirements. The exact flight condition should be set to the state with the most 
decreased damping coefficients. 


4.3, Selection of the Flying Qualities Evaluation Criteria 


Current flying qualities standards provide multiple evaluation criteria [35-38]. Each 
evaluation criterion has its own inspection focus, and thus, the same configuration design 
may yield different evaluation results according to different criteria. To determine con- 
figuration parameter boundaries based on closed-loop flying qualities requirements, it is 
necessary to determine the evaluation criteria that are most sensitive to changes in the 
configuration parameters to obtain the strictest parameter design ranges. 


4.3.1. CAP Criterion 

The control anticipation parameter (CAP) is defined as the ratio of the initial pitch 
acceleration to the steady-state normal acceleration [38]. The CAP value can reflect the 
short-period mode flying qualities of an aircraft, including the longitudinal stability and 
maneuverability. To apply the CAP criterion to a high-order closed-loop aircraft, the low- 
order equivalent system (LOES) theorem is adopted to match the characteristics of the 
high-order system. The longitudinal LOES of the aircraft can be expressed as: 

Ka 


wo Ket 
F s* + 2spWnsps + Ons 
q K,(s +1/Tp,) 
Fe s* + 2 spWnsp$ + Wisp 


fa?) 


(8) 


— TqS 


where « is the aircraft angle of attack; q is the pitch rate; F. is the pitch control input; 1/To, 
is the zero of the pitch attitude transfer function; Ky and K, are the equivalent system gains 
of the transfer functions; and Ty and T, are the equivalent time delays. 

The CAP of the LOES is defined as: 


ay to - 


where V is the flight speed and g is the acceleration due to gravity. The Level 1 FQ of the 
CAP criterion are defined by the short-period damping ratio ¢sp of the equivalent system 
and the CAP value. 

The upper boundary of the CAP value and the lower boundary of ¢sp reflect the 
stability requirements of an aircraft. A large CAP value and small Csp indicate a lack of 
longitudinal stability, and thus, a pilot’s small control input may lead to a violent pitch 
response. Therefore, the upper boundary of the CAP value and the lower boundary of sp 
can be used to determine the forward range of the wing position, the minimum limit of the 
horizontal tail area and the rearward CG limit. 

In contrast, the lower boundary of the CAP value and the upper boundary of ¢sp 
reflect the maneuverability requirements. A small CAP value and large Csp represent a 
large longitudinal stability margin, which may increase the pilot’s workload to achieve 
pitch motion. Thus, the rearward range of the wing position, the maximum limit of the 
horizontal tail area and the forward CG limit should be designed according to the lower 
boundary of the CAP value and the upper boundary of sp. 


4.3.2. Chalk Criterion 


The Chalk criterion is a time-domain criterion that focuses on the characteristics of 
the pitch rate response and defines three parameters for evaluation: the effective time 
delay t,, the transient peak ratio Aqmax/Agmin and the effective rise time At, as shown in 
Figure 3 [38]. 
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Figure 3. Definitions of the Chalk criterion parameters. 


In Figure 3, gss is the steady-state pitch rate; gmax is the maximum pitch rate; qmin 
is the first minimum pitch rate; the instant ft, is defined by the intersection between the 
tangent at the maximum slope and the time axis; the instant tz is defined by the intersection 
between the tangent at the maximum slope and the horizontal line of the steady-state pitch 
rate; and At is the difference between f; and fp. t; reflects the time delay between the real 
pitch rate response and an ideal pitch rate response with the maximum pitch acceleration, 
Agmax/ Agmin reflects the short-period damping ratio of the pitch response, and At reflects 
the rise time to accelerate to the steady-state pitch rate of the pitch rate ideal response. 

The Chalk criterion is based on the time-domain pitch response of the aircraft without 
identifying the frequency-domain characteristics. Thus, the Chalk criterion can be used 
to avoid the situation in which a frequency-domain criterion is not applicable to aircraft 
with an unconventional response. The design ranges of the configuration parameters 
should be determined according to the more stringent requirements of the CAP and 
Chalk criteria to enhance the universality of the design method based on closed-loop 
flying qualities requirements. 


4.3.3. Dutch Roll Criterion 


The flying qualities with respect to the yaw axis can be evaluated by the characteristics 
of the lateral-directional oscillatory (Dutch roll) response to the yaw controller. The Dutch 
roll frequency and damping ratio can be defined by matching the higher-order sideslip 
response to the yaw control input to the following lower-order form [38]: 


K = 
B = +e T° (10) 
Fr (s° + 2 grWy.drS + W% a.) 


where f is the sideslip angle; F, is the yaw control input; Kg is the equivalent system gain 
of the transfer function, and T,, is the equivalent time delay of the transfer function. 

The flying qualities requirements of the Dutch roll criterion are specified in terms 
of the minimum values of Wy gr, Gar and €g;Wn.a,r, and the evaluation results can reflect 
the lateral-directional stability of the aircraft. This criterion can be used to determine the 
vertical tail area and position boundaries. 


5. Application of the Proposed Parameter Boundaries Design Method 


This section takes a medium-range transport aircraft as an application example of 
the configuration parameter boundaries design method. The definitions of the target 
configuration parameters of the sample aircraft are illustrated in Figure 4. The wing 
position xy and the vertical tail position xy are defined as the distance from the nose to the 
quarter chord point on the mean aerodynamic chord (MAC) of the wing and the vertical 
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tail. The CG position x; is defined as the percentage of the distance from the CG to the 
leading edge of the wing MAC [1]. 





Figure 4. Target configuration parameters of the sample aircraft. 


The configuration parameters of the basic aircraft design are depicted in Table 1. The 
aircraft motion model and the inner loop of the NDI control law are based on the basic 
aircraft design, and the reference model parameters of the NDI control law are designed 
according to the flying qualities requirements, as shown in Table 2. By combining the 
motion model with the FCL, the high-order closed-loop aircraft model is established. 


Table 1. Configuration parameters of the basic aircraft design. 


Parameters Values 

Length overall (m) 32 
Wing position (m) 12.6 
Wing area (m7) 115.3 

Wing mean aerodynamic chord (MAC) (m) 6.2 

Center of gravity 33.4% (TO) 38.2% (LD) 

Horizontal tail position (m) 28.6 
Horizontal tail area (m2) 279 

Horizontal tail MAC (m) 2.9 
Vertical tail position (m) 29.8 
Vertical tail area (m*) 15.9 

Vertical tail MAC (m) o2 


Table 2. Reference model parameters design. 


Wsp.r To2+ Wp.r Wrz 
Kar Csp.x (rad/s) (s) Kps (rad/s) (rad/s) 
75 0.9 3.5 0.5 42 25 3.5 


5.1. Boundaries of the Wing Position and Horizontal Tail Area 


For the sample aircraft, the forward CG limit appears in the take-off phase (sea-level, 
low-speed, forward CG position, maximum take-off weight, maximum thrust, flaps and 
landing gear retracted), while the rearward CG limit appears in the landing phase (sea- 
level, low-speed, rearward CG position, minimum landing weight, landing engine setting, 
flaps and landing gear down). According to the previous discussion, the forward range of 
xw and the minimum limit of Sy; should be designed by evaluating the flying qualities in 
landing phase according to the most stringent requirement among the upper boundary of 
the CAP value, the lower boundary of ¢sp, and the Chalk criterion. Likewise, the rearward 
range of xw and the maximum limit of S; should be designed in take-off phase according 
to the most stringent requirement among the lower boundary of the CAP value, the upper 
boundary of ¢sp and the Chalk criterion. 

Preliminary sizing requires that the CG position is on the MAC of the wing; thus, the 
first estimate of xy ranges from 10 m to 16 m. Then, according to the requirements of the 
longitudinal Level 1 FQ of closed-loop aircraft, the design range of S; corresponding to 
different xw is determined. Since the rearward range of xw and the maximum limit of Sr 
both reflect the maneuverability requirements, the upper boundary of S7 corresponding to 
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different xy can also be regarded as the rearward boundary xw corresponding to different 
Sr. Similarly, the lower boundary of S; corresponding to different xw can also be regarded 
as the forward boundary of xw corresponding to different S;. The design range of xy is 
determined by the design range of Sr. 


e Maximum limit of the horizontal tail area 


The variation series of xy is established as follows: within the preliminary range from 
10 m to 16 m, one point is selected every 1 m, yielding seven points in total. For each Xw, 
a series of Sp with an interval of 2 m? is selected. Then, the CAP criterion is applied to 
evaluate the flying qualities of aircraft with different xy—S7; combinations in the take-off 
phase. Partial evaluation results with xy values of 10 m, 12 m, 14m, and 16 m are shown 
in Figure 5. The CAP criterion characteristics with xw values of 12 m and 14 m are shown 
in Table 3. 
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Figure 5. CAP criterion evaluation results in take-off phase. 


Table 3. CAP criterion evaluation results with xy = 12m and 14 m in the take-off phase. 


xw(m) Sr (Mm?) sp CAP Level xw(m) S7(m’) Csp CAP Level 
12 18 0.96 3.76 2 14 14 1.01 3.07 2 
12 20 0.99 3.56 i 14 16 1.05 3.35 | 
12 PP. 1.04 3.00 1 14 18 1.08 2.96 1 
12 24 1.08 202 1 14 20 1.11 249 il 
12 26 1.13 2.49 1 14 22 1.15 2.57 | 
12 28 1.16 2.06 1 14 24 1.20 2.36 1 
12 30 1,17 1.75 1 14 26 1.24 2.05 | 
12 32 1.21 1.55 1 14 28 1.26 1.79 1 
12 34 1.2) 1.29 iL 14 30 1.28 1.52 | 
12 36 1.29 0.98 1 14 32 1.29 1.27 1 
12 38 132 0.82 2 14 34 1.32 1.09 2 


For a certain xw, the CAP value decreases gradually with increasing S7, and the 
damping ratio ¢sp increases. The maximum limit of Sy defined by the CAP criterion is the 
critical value that leads to the degradation of the flying qualities to Level 2. When xy is set 
to 12 m and 14 m, the corresponding maximum limit of S is 36 m? and 32 m?, respectively. 
Hence, as the wing position moves rearward, the maximum limit of Sy gradually decreases. 

According to all the above evaluation results, the level boundaries of the xw-—Sr 
combination defined by the CAP criterion in take-off phase are presented in Figure 6. 
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According to this figure, to ensure the Level 1 FO of the CAP criterion, when xy is set 
as 10 m, the corresponding design range of Sr is from 24 m* to 38 m*; when xy is set as 
16 m, the S; design range is from 14 m? to 26 m*. The right boundary of Level 1 in Figure 6 
determines the maximum limit of Sy corresponding to different xy, and the rearward 
range of xw corresponding to different Sr. 
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Figure 6. Parameter level boundaries defined by the CAP criterion in the take-off phase. 


The Chalk criterion is also applied to evaluate the flying qualities in the take-off phase. 
Figure 7 shows the influence of changing Sr on the pitch rate response with xw = 12 m. 
The variations in the tail area exert the main influences on the response peak and response 
speed of the pitch rate. As Sr increases, the response peak decreases, and the response 
time delay increases. The influence of the rearward displacement of the wing is similar to 


that of changing S;, which means that decreasing xw will reduce both the response peak 
and the response speed. 
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Figure 7. Pitch rate responses with different Sr. 


Among all Chalk criterion parameters, the effective time delay t; is the most affected 
by the variations in and Sr. The evaluation results of t; are shown in Table 4. The transient 
peak ratio Aqmax/ Agmin and the effective rise time At are also affected, but the FQ Level is 
not degraded. With the rearward movement of the wing and an increase in the tail area, 
the t; of the pitch response increases, which indicates poor maneuverability. The FQ Level 
is degraded to Level 2 (0.12 s < t; < 0.17 s) or even Level 3 (#1; > 0.17 s). In the chosen ranges 


of xw and Sr, an increase in xw and a decrease in Sr will not degrade the flying qualities 
of the Chalk criterion. 
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Table 4. Chalk criterion evaluation results in the take-off phase (t;, seconds). 
ay 18 20 22 24 26 28 30 32 34 36 38 40 
xw (m) 
10 0.049 0.055 0.065 0.071 0.076 0.083 0.091 0.098 0.107 0.112 0.118 0.129 
11 0.055 0.063 0.072 0.082 0.082 0.095 0.103 0.105 0.109 0.115 0.125 On37 
12 0.063 0.071 0.081 0.091 0.089 0.102 0.109 0.109 0.112 0.118 0.132 0.152 
13 0.075 0.079 0.089 0.098 0.095 0.107 0.112 0.113 0.117 0.142 0.156 0.169 
14 0.081 0.087 0.095 0.103 0.106 0.113 0.116 0.118 0.135 0.158 0.171 0.182 
15 0.092 0.096 0.101 0.109 0.112 0.118 0.128 0.145 0.152 0.167 0.185 0.196 
16 0.098 0.102 0.107 0.115 0.119 0.129 0.142 0.155 0.163 0.172 0.192 0.212 
Background Level 1 Level 2 Level 3 


The level boundaries of the xw-S7 combination defined by the Chalk criterion in the 
take-off phase are shown in Figure 8. According to this figure, to ensure the Level 1 FQ of 
the Chalk criterion, the maximum limit of S7 is 38 m* when xw = 10 m, and the maximum 
limit is 26 m? when xw = 16 m. Compared with the parameter boundaries defined by 
the CAP criterion shown in Figure 6, the Chalk criterion defines a more relaxed Level 1 
boundary of the rearward range of xy and the maximum limit of S;. In contrast, the Level 
2 boundary defined by the Chalk criterion is stricter and can be used as a supplement to 
the CAP criterion boundaries. 
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Figure 8. Parameter level boundaries defined by the Chalk criterion in the take-off phase. 


e Minimum limit of the horizontal tail area 


The selection of the xw and Sr series is the same as the design of the maximum 
Sr range. The CAP criterion is applied to evaluate the flying qualities in the landing 
phase. Partial evaluation results with xy values of 10 m, 12 m, 14m, and 16 m are shown 
in Figure 9, where the corresponding minimum limits of Sp are 26 m7, 22 m7, 20 m?, 
and 16 m’, respectively. As the wing position moves forward, the minimum limit of Sr 
gradually increases. 

The level boundaries of the xw-S7 combination defined by the CAP criterion in 
landing phase are shown in Figure 10. According to this figure, to ensure the Level 1 FQ of 
the CAP criterion, the design range of Sr is from 26 m* to 40 m? when xy = 10 m and from 
16 m* to 32 m* when xy = 16 m. The left boundary of Level 1 in Figure 10 can determine 
the minimum limit of Sy corresponding to different xy, and the forward range of xw 
corresponding to different Sr. 
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Figure 9. CAP criterion evaluation results in the landing phase. 
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Figure 10. Parameter level boundaries defined by the CAP criterion in the landing phase. 


According to the evaluation results of the Chalk criterion in the landing phase, the 
effective rise time and the response time delay decrease as the wing moves forward and the 
tail area decreases, which represents a better response speed of the pitch rate. Meanwhile, 
the pitch damping characteristics become worse as the transient peak ratio Agmax / AGmin 
increases, but the FQ Level is not degraded. Therefore, the forward boundary of xw and 
the minimum limit of Sy should be determined only by the CAP criterion. 


e Level boundaries of the wing position and horizontal tail area 


To determine the design ranges of xw and Sy according to the Level 1 FQ requirements 
in the entire flight envelope, the most stringent boundaries in Figures 6, 8 and 10 are 
combined to form the final parameter level boundary diagram, as shown in Figure 11. 

According to Figure 11, the Level 1 boundary of the rearward xw and the maximum 
Sr (right side in Figure 11) is determined by the evaluation results in the take-off phase; 
the upper half is defined by the Chalk criterion and the lower half is defined by the CAP 
criterion. The Level 2 boundary (top right in Figure 11) is determined in take-off phase by 
the Chalk criterion. The Level 1 boundary of the forward xw and the minimum Sr (left 
side in Figure 11) is defined in the landing phase by the CAP criterion. The design range of 
Sr determined by the Level 1 FQ requirements is from 26 m* to 38 m* when xw = 10 mand 
from 16 m2 to 26 m* when xw = 16 m. In the configuration design process, the design range 
of Sr can be defined as the horizontal interval within the Level 1 zone (shaded part in 
Figure 11) corresponding to the designed xy. Similarly, the design range of xy is defined 
as the vertical interval within the Level 1 zone according to the designed Sr. 
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Figure 11. Complete level boundaries of xy and Sr. 


5.2. Boundaries of the Center of Gravity Position 


The forward limit of X-¢ should be designed in the take-off phase according to the most 
stringent requirement among the lower boundary of the CAP value, the upper boundary 
Of Esp, and the Chalk criterion. Likewise, the rearward limit of X-¢ should be designed 
by evaluating the flying qualities in the landing phase according to the most stringent 
requirement among the upper boundary of the CAP value, the lower boundary of ¢sp, and 
the Chalk criterion. 

The series for the forward Xc¢ limit design is set as follows: from 6% to 36%, a point 
is selected every 3% of the wing MAC, yielding ten points in total. Likewise, the series 
for the rearward Xc¢ limit design is set as follows: from 38% to 65%, a point is selected 
every 3% of the wing MAC, also yielding ten points in total. Then, the CAP criterion is 
applied to evaluate the flying qualities of the aircraft with the forward X¢. in the take-oft 
phase, while the rearward X;¢ designs are evaluated in the landing phase. The evaluation 
results are presented in Figure 12. Compared with the basic design of the sample aircraft 
in the take-off phase, as the CG position moves forward (TO phase, X¢.—), the CAP value 
decreases gradually, the damping ratio Csp increases, and the FQ Level is degraded from 
Level 1 to Level 2. In the landing phase, as the CG position moves rearward (LD phase, 
Xce+), the CAP value increases, and the damping ratio ¢sp decreases. 

The forward X-, designs are also evaluated by the Chalk criterion in the take-off 
phase. The effective time delay ft; is the parameter most affected by the variation in X¢¢; 
the evaluation results of f; are shown in Table 5. As CG moves forward, the f; value of the 
pitch response increases, which indicates poor maneuverability. The FQ Level is degraded 
to Level 2 or even Level 3. 


Table 5. Chalk criterion evaluation results of the forward Xc¢g (t1, seconds). 





Xcg (% MAC) t, (s) Level Xcg (% MAC) t, (s) Level 
33 0.05 1 18 0.121 2 
30 0.061 1 15 0.137 2 
27 0.074 1 12 0.155 2 
24 0.088 1 9 0.172 3 
21 0.103 1 6 0.189 3 


According to the evaluation results of the Chalk criterion in the landing phase, the 
pitch damping characteristics of the aircraft become worse as the CG moves rearward, but 
the FO Level is not degraded. Therefore, the Chalk criterion is only applied to design the 
forward Xc¢g limit. 

The level boundaries of X-. are defined as in Figure 13. According to this figure, the 
Level 1 boundary of the forward CG limit is 21% of the wing MAC, determined by the 
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upper boundary of ¢sp in the CAP criterion in the take-off phase. The Level 2 boundary of 
the forward CG limit is 12% of the wing MAC, determined by the Chalk criterion. Finally, 
the Level 1 boundary of the rearward CG limit is 59% of the wing MAC, determined by 
the upper boundary of the CAP value in the landing phase. The X;. range should be kept 
within the Level 1 zone (shaded part in Figure 13). 
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Figure 12. CAP criterion evaluation results of Xcg. 
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Figure 13. Level boundaries of Xcg¢. 


5.3. Boundaries of the Vertical Tail Area and Position 


The vertical tail area Sy and position xy should be designed in the cruise phase (high 
speed, high altitude, cruise engine setting and configuration) according to the requirements 
of the Dutch roll criterion. The variation series of Sy is set as follows: below the basic 
value of Sy, a point is selected every 1 m7, and for each Sy, a series of xy with an interval 
of 1 mis selected. As xy is adjusted by changing the quarter chord point of the vertical 
tail without moving the root position, the design range of xv is limited from 26 m to 31 m. 
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Then, the Dutch roll criterion is applied to evaluate the flying qualities of different Sy—xy 
combinations. Partial evaluation results with Sy values of 10 m2, 12 m2, 14 m’, and 16 m2 
are shown in Figure 14, and the Dutch roll characteristics with Sy values of 12 m? and 


14 m2 are shown in Table 6. 
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Figure 14. Dutch roll criterion evaluation results in the cruise phase. 


Table 6. Dutch roll criterion evaluation results with xy = 12m, 14 m in the cruise phase. 





Sy (m2). xy (m) Cdr Wr.dr (rad/s) Level Sy (m2) xy (m) Ca Wr.dr (rad/s) Level 
DD 31 0.142 TA32 1 14 31 0.166 1.153 1 
12 30 0.138 1.107 1 14 30 0.154 L127 1 
12 29 0.115 1.012 2 14 29 0.150 1.093 1 
jw 28 0.104 0.937 2 14 28 0.145 1.044 il 
12 27 0.097 0.878 2 14 27 0.132 0.981 2 
12 26 0.091 0.835 2 14 26 0.117 0.920 2 


For a certain Sy, both wy.q, and ¢g; of the Dutch roll mode decrease gradually with the 
forward movement of the tail. The forward range of xv is defined as the critical value that 
leads to the degradation of directional flying qualities. When Sy is set as 12 m? and 14 m’, 
the corresponding forward range of xy is 30 m and 28 m, respectively. With decreasing 
Sy, the design range of xy gradually decreases. The minimum limit of Sy is defined as the 
critical value without the corresponding xy that meets the requirements of the Level 1 FQ. 

According to all evaluation results, the level boundaries of the Sy—xy combination are 
defined as in Figure 15. 

Considering that the drag force due to the vertical tail is largely related to the ver- 
tical tail area, it is necessary to restrict the tail area while satisfying the flying qualities 
requirements. According to Figure 15, xy within the design range can ensure Level 1 FQ 
of the closed-loop aircraft when Sy is not less than 15 m*. The minimum limit of Sy is 
11 m? and the corresponding xy is 31 m. In the configuration design process, the design 
range of Sy can be defined as the horizontal interval within the shaded Level 1 zone in 
Figure 15, corresponding to the designed xy. Similarly, the design range of xy is defined 
as the vertical interval within the Level 1 zone according to the designed Sy. 
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Figure 15. Complete level boundaries of Sy and xy. 


6. Conclusions 


In this study, the influences of typical configuration parameters on flying qualities are 
analyzed, and a configuration parameter boundaries design method based on closed-loop 
flying qualities requirements is proposed and verified by an application example. 

(1) Five typical configuration parameters that should be designed based on flying 
qualities requirements are determined: wing position, horizontal tail area, center of gravity 
position, vertical tail area and position. 

(2) The proposed configuration parameter boundaries design method based on closed- 
loop flying qualities requirements is achieved as follows. First, the flight control law is 
designed based on the model reference nonlinear dynamic inversion method, and a high- 
order closed-loop aircraft model is established. Second, the most stringent flight conditions 
and the evaluation criteria that are most sensitive to the target parameters are determined. 
Third, by searching for the critical values that lead to the degradation of flying qualities, 
the closed-loop flying qualities requirements are transformed into the design ranges of 
target configuration parameters. 

(3) An application example involving the typical parameter boundaries design of a 
sample aircraft is provided. The wing position of the sample aircraft ranges from 10 m 
to 16 m, the design range of the horizontal tail area corresponding to the forward limit of 
wing position is from 26 m* to 38 m?, and the horizontal tail area range corresponding 
to the rear limit of wing position is from 16 m? to 26 m*. The design ranges of the wing 
position and horizontal tail area are defined by the level boundary diagram. The forward 
limit of the center of gravity is 21% of the average aerodynamic chord, while the rearward 
limit is 59%. The minimum limit of the vertical tail area is 11 m*, and the design ranges of 
the vertical tail area and position are defined by the level boundary diagram. 
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Abbreviations: 

CG Center of Gravity 

FCL Flight Control Laws 

SAS Stability Augmentation System 

SCAS Stability and Control Augmentation System 
NDI Nonlinear Dynamic Inversion 

INDI Incremental Nonlinear Dynamic Inversion 
FQ Flying Qualities 

PI Proportional-Integral 

CAP Control Anticipation Parameter 


LOES Low-Order Equivalent System 
MAC Mean Aerodynamic Chord 


TO Take-Off 
LD Landing 
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Abstract: The introduction of disruptive innovations in the transport aviation sector is becoming 
increasingly necessary. This is because there are many very demanding challenges that the transport 
aviation system will have to face in the years ahead. In particular, the reduction in pollutant 
emissions from air transport, and its impact on climate change, clearly must be addressed; moreover, 
sustainable solutions must be found to meet the constantly increasing demand for air traffic, and 
to reduce the problem of airport saturation at the same time. These three objectives seem to be in 
strong contrast with each other; in this paper, the introduction of a disruptive airframe configuration, 
called PrandtlPlane and based on a box-wing lifting system, is proposed as a solution to face these 
three challenges. This configuration is a more aerodynamically efficient alternative candidate to 
conventional aircraft, introducing benefits in terms of fuel consumption and providing the possibility 
to increase the payload without enlarging the overall aircraft wingspan. The development and 
analysis of this configuration, applied to a short-to-medium range transport aircraft, is carried out 
through a multi-fidelity physics-based approach. In particular, following an extensive design activity, 
the aerodynamic performance in different operating conditions is investigated in detail, the structural 
behaviour of the lifting system is assessed, and the operating missions of the aircraft are simulated. 
The same analysis methodologies are used to evaluate the performance of a benchmark aircraft with 
conventional architecture, with the aim of making direct comparisons with the box-wing aircraft 
and quantifying the performance differences between the two configurations. Namely, the CeRAS 
CSR-01, an open-access virtual representation of an A320-like aircraft, is selected as the conventional 
benchmark. Following such a comparative approach, the paper provides an assessment of the 
potential benefits of box-wing aircraft in terms of fuel consumption reduction and increase in payload 
capability. In particular, an increase in payload capability of 66% and a reduction in block fuel per 
pax km up to 22% is achieved for the PrandtlPlane with respect to the conventional benchmark, while 


maintaining the same maximum wingspan. 


Keywords: box-wing; PrandtlPlane; disruptive aircraft; innovation; future aviation 


1. Introduction 


In the near future, the aviation industry will have to face several challenging tasks; the 
most relevant issue is related to the environmental impact of transport aviation, which has 
to drastically reduce its pollutant footprint and its impact on climate change [1-6]. At the 
same time, the demand for commercial air transport is growing quickly [7-10], whereas 
airports are gradually addressing saturation problems using aprons [11]. Meeting the 
erowing demand for flights and reducing the environmental impact of air transportation 
seem to be contrasting objectives of aviation research and industry; a possible way to face 
this problem is to explore and study disruptive technologies, both in the field of propulsion 
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and aircraft design. Regarding propulsion, several studies are dedicated to the analysis and 
design of new types of engines, alternative to the conventional internal combustion one, 
such as hybrid and electrical powertrains [12-17]; also, alternative energy sources to fossil 
fuels, such as hydrogen [18-21] or biofuels [22—24], are currently under study. Regarding 
aircraft development, studies have been conducted on the optimisation of the tube-and- 
wing aircraft to reach the maximum potential of the conventional configuration [25-28], but 
ereat effort is currently dedicated to the study and development of new unconventional 
architectures [29-34], such as the blended wing-body configuration [35-37] or the joined 
wings architecture [38,39], which may represent a real breakthrough in aircraft evolution. 

In this paper, the focus is on an unconventional configuration called PrandtlPlane [40-43], 
a concept based on the box-wing lifting architecture, which is represented in Figure 1; for 
the PrandtlPlane (PrP) configuration, the box-wing is designed according to Prandtl’s 
“best wing system” theory [44], which allows one to minimise the induced drag of the 
lifting system [45,46]. This configuration has been extensively studied in the PARSIFAL 
project [47,48], a research project funded by the European Commission in the framework 
of the Horizon2020 program. 





Figure 1. Artistic views of PrandtlPlane aircraft of different categories: ultra-light amphibious 
(left) [43], executive class (centre) [49], hybrid electric regional airliner (right) [15]. 


Some of the main challenges that will be faced by aeronautical research in the near 
future are summarised in Table 1, together with some general possible solutions. 


Table 1. Challenges ahead for the transport aviation and related general solutions. 


Challenge Possible Solution 


To meet the large air traffic demand increase expected in the coming 
years, in particularly for short/medium routes [7] *. 


To design an aircraft with an increased cabin capacity 
compared to the present aircraft operating on 
short/medium routes. 


To avoid airport saturation problems, already relevant today [11]. To limit the size and overall dimensions of the aircraft. 





To reduce the environmental impact of the aircraft [2], thus 
minimising fuel consumption per passenger. 


To increase the aerodynamic efficiency as much as possible 
and/or to adopt new types of propulsion (i.e., electric 
or hydrogen). 


* The forecasts used as reference were provided before the COVID-19 outbreak; the pandemic will have a direct impact on the air traffic 
demand and on the related forecast. The main consequence is related to a shift-forward of the peak increase in the demand with respect the 
previous market outlooks. However, the overall scenario is supposed to be the same, in a delayed time scale. 


To effectively address the issues ahead, aircraft configuration must comply with the 
three indicated solutions at the same time. Concerning the issue of the saturation of airport 
areas, it should be pointed out that using larger existing aircraft to reduce the number of 
flights does not represent a solution to the problem; for example, using an Airbus A350 in 
place of two Airbus A321s on medium-haul flights (leaving aside the fact that A350 aircraft 
are not designed to accomplish these missions) would reduce the number of slots occupied, 
but the size of the Airbus A350 would require the occupation of a much larger apron (ICAO 
Aerodrome Code C: max wingspan 36 metres vs. [CAO Aerodrome Code E: max wingspan 
65 m) and taxiway. It is clear that this operation does not bring any advantages. On the 
contrary, it would imply further disadvantages, considering the fact that the majority of 
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worldwide airport slots are represented by ICAO C aprons, while there are significantly 
fewer ICAO E aprons [50].Wingspan definitely plays a key role in ground operations, as it 
is constrained by the standardised dimensions of the apron and taxiways. 

The general solutions proposed to the three challenges presented above indicate that 
a technological step forward is necessary for commercial aviation aircraft: it is difficult to 
try to solve these problems through the adoption of aircraft with conventional architecture 
and/or conventional propulsion system [51]. 

In this work, a suitable solution to the presented challenges is proposed through 
the introduction of the non-conventional PrandtlPlane architecture, while aspects related 
to innovative propulsion systems are not considered; in particular, a detailed analysis is 
conducted to answer the question: “can the introduction of the PrandtlPlane configuration 
be a possible solution to face the challenges of commercial aviation in the near future?” 

The PrandtlPlane configuration is a candidate to satisfy the three challenges previously 
presented; in particular, considering the short/medium-haul market sector, it allows: 


1. The transport of a larger number of passengers compared to the present aircraft 
operating on medium routes, by exploiting the increased lifting capacity of the box- 
wing system, and using a new fuselage design, in order to increase the number of 
travellers without increasing the number of flights; 

2. The exploitation of the increased lifting capacity of the box-wing to design an aircraft 
with the same overall dimensions, in particular of wingspan, of present aircraft 
operating on short/medium routes, while improving the passenger capability, in 
order to avoid the increase in required apron space; namely, this requires the design 
of a box-wing with wingspan limited to the standard related to short-to-medium- 
route aircraft, i.e., maximum 36 m (ICAO Aerodrome Code C), to be compliant with 
the airport infrastructure (aprons, taxiways), but at the same time, to transport a 
larger payload; 

3. The exploitation of the box-wing architecture designed according to the “Best Wing 
System” theory to maximise the aerodynamic efficiency; indeed, a properly designed 
box-wing aircraft allows one to theoretically minimise the induced drag, by exhibiting 
Oswald factor efficiencies larger than 1 [45], as well as to increase the lifting capability 
(i.e., to trim a larger weight), and thus to reduce the fuel consumption per passenger— 
kilometre compared to conventional aircraft. 


In this paper, first of all, the multidisciplinary design methodology adopted to develop 
a short/medium-haul PrandtlPlane aircraft is presented by describing the multi-fidelity 
design and analysis strategy implemented in detail. The main focus is given to the aero- 
dynamic and structural design of the innovative box-wing lifting system, together with a 
detailed assessment of the mission performance. Then, the PARSIFAL project case study is 
presented and described in detail, starting from the definition of the design requirements to 
the aerodynamic and mission performance evaluation. The box-wing aircraft performance 
are then compared with those of a reference conventional tube-and-wing benchmark air- 
craft; the performance assessment methods used are the same in both cases, to provide 
direct and fair comparisons between two totally different aircraft architectures. 


2. Design Methodology Description 


The design strategy adopted for the development of the box-wing aircraft is schemati- 
cally summarised in Figure 2; three main different stages were performed: a conceptual 
design initialisation, a reference aerodynamic layout definition through aerodynamic op- 
timisation and a subsequent model refinement and performance evaluation. These three 
main blocks of the design process are detailed in the following. 
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Figure 2. Multi-level aircraft design and analysis diagram (abbreviations used: VLM = Vortex Lattice 
Method, FEM = Finite Element Model, CFD = Computational Fluid Dynamics). 


2.1. Conceptual Design Initialisation 


The first step of the initial conceptual sizing is represented by the definition of the 
fuselage layout, its overall size, and the characterisation of the cabin section layout and 
the internal arrangement; this is fundamental, because one of the main requirements to 
be satisfied is represented by the maximum number of passengers: to accomplish this 
specification, a proper fuselage design is needed to start the following aircraft design 
process. Then, the conceptual design initialisation mainly consists of a rough evaluation 
of the aircraft aerodynamic performance and the related estimation of the initial guess 
of the maximum take-off weight. The inputs of the procedure are the Top Level Aircraft 
Requirements (TLARs), namely the design payload and range, the operating conditions 
and the size constraints. As far as the aerodynamic performance is considered, a simplified 
model to estimate the aircraft drag was used by separately calculating the induced drag 
coefficient, Cp;, the parasite drag coefficient, Cpg, and the wave drag coefficient, Cowave, of 
the aircraft (Equation (1)). 


Cp tot — Chi F Cpo =i Cp wave (1) 


The induced drag coefficient, Cp; , is estimated by considering only the lifting compo- 
nents (Equation (2)); namely, the contribution of the front and rear wings are considered, as 
the contributions to the induced drag of the vertical wings and the fuselage are neglected 
at this stage. 


Front R Vertic Fy 
Coi= Cpi + Cp + Cpe + Cog (2) 


The induced drag coefficient, Cp;“"”"S , of each lifting component is evaluated by using 
the relations of Equations (3)—(4), where Cy, represents the wing lift coefficient, AR is the 
related aspect ratio, and e is the Oswald factor. As the considered configuration is a box- 
wing designed according to the Prandtl’s theory on the best wing system [44], the Oswald 
factor, e, used in this phase is fixed equal to 1.46 [52]. Since the aircraft geometry does 
not exist at this initialisation stage, some reference parameters must be set reasonably, as 
indicated in the reference literature. This is the case for the Oswald factor, as, although it 
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depends on the wings’ height to span ratio, h/b [44], it is set as equal to the theoretical 
maximum value for a box lifting system [52]. This does not represent a misalignment, 
since in the following design phases (Section 2.2) the Oswald factor is calculated for 
each configuration considered (taking into account its geometrical and lift distribution 
characteristics) by means of the VLM aerodynamic solver. 


Cp" = kCy? (3) 


k = 1/(7 ARe) (4) 


The parasite drag coefficient, Cpg, is evaluated by the equivalent skin friction model 
described in [53]; namely, the Cpg tor of the aircraft is the sum of each aircraft component, 
Coo", evaluated by means of the formula of Equation (5). 


Coo"? — Cre (S Oe a S ref ) (5) 


Cpotr= Cpe (6) 


where Cy, is the equivalent skin friction coefficient, equal to 0.0030 for civil transport aircraft 
from the approximation reported in [53]; the same reference presents a simplified way to 
estimate the components’ wetted surface, Swe", both in the case of a wing component 
(i.e., lifting surfaces, vertical tip-wing, tail, fin, Equation (7)) and cylindrical components 
(i.e., fuselage, nacelle, Equation (8)): 


Sue = S exposed LOOT 0.52 (t/c)| (7) 


Suet! = oF doy! ] cyl (8) 


where Sexposed is the planform area of the wing exposed to the flow, (t/c) is the wing 
thickness-to-chord ratio, and d! and I are the diameter and the length, respectively, of 
the considered cylinder. For aircraft cruising in transonic flight, the wave drag coefficient, 
Cpwave, 1S fixed as constant and equal to 0.002, according to the Boeing definition of transonic 
drag rise reported in [53]. 

Several simplified models are presented in the literature to evaluate the aircraft gross 
weight in the very initial stage of the aircraft design. In this phase, to initialise the design 
process, the aircraft maximum take-off weight, MTOW, is estimated by individually evalu- 
ating the contribution of the payload weight, Wpay, the fuel weight, W,,1, and the operating 
empty weight, Woe (Equation (9)): 


MTOW = Woay+W fuel + W oe (9) 


where, concerning the payload, 95 kg per passenger is assumed ([54]), the fuel weight, 
Wruel, is evaluated by using the well-known Breguet formula [55], and the operating empty 
weight, Woe, is estimated by the statistic interpolation of historical trends proposed in [53]. 


2.2. Optimisation Driven Preliminary Design 


The main information obtained from the conceptual design initialisation is the first 
estimation of the aircraft maximum take-off weight, which is useful to start the second 
level of the aircraft design process (Figure 2). In particular, in this stage, the aerodynamic 
design is driven by an optimisation procedure implemented in an in-house tool called 
AEROSTATE [56-58]; the schematic representation of the AEROSTATE code is reported 
in Figure 3. 
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Figure 3. AEROSTATE workflow. 


The general constrained aerodynamic optimisation problem addressed in this tool is 
summarised in the following Equations (10)—(12): 


e Objective function 


min(—L(x)/D(x)) (10) 
e Set of inequality constraints 
g(x) = 0 (11) 
e Design space 
lb <x < ub (12) 


where the objective function to be minimised is the -L/D ratio, which is equivalent to the 
maximisation of the aerodynamic efficiency (lift-to-drag ratio) evaluated in a reference 
point of the cruise (e.g., the 25% of the stage length); the optimisation is constrained by 
a set of inequality constraints 9(x); the design variables x, limited by lower and upper 
boundaries (Ib and ub, respectively), are reported in Figure 4. 
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Figure 4. Design variables of the aerodynamic optimisation. 


Aerospace 2021, 5, 292 


7 of 37 


The optimisation tool allows one to introduce constraints to search feasible configura- 
tions, from aerodynamic and aeromechanic standpoints; the set of inequality constraints 
are summarised in Equations (13)-(19): 


e =6© Vertical trim 


Waes —€L < L(x) = W des FEL (13) 


e Longitudinal static stability 


SSMmin < SSM(x) < SSMmax (14) 


e Pitch trim 


IA 


—em < Cm(x) < emu (15) 


e Max local lift coefficient 


max(ci(y, *)) S< Cir (16) 
e Wing loading constraints 
(L/S)enin < (E(®) wing/S(®) wing) < (L/S) max (17) 
e Taper ratio 
A; (x) oe (18) 
e Relative wings position 
ront wing ront wing rear Wing 
cA tip +cy> < Xrp tip +Ax (19) 


where W,,, is the input design weight; SSM is the longitudinal static stability margin; Cjy is 
the pitching moment coefficient; c;(y) is the spanwise distribution of local lift coefficient and 
C) ry is the related threshold value; (L,,,,,, el Cue) is the i-th lifting surface wing loading; A; 
is the j-th wing bay taper ratio; x; F jj, is the longitudinal position of the wing tip leading 
edge; Ct, is the wing tip chord; ¢ is a tolerance. Several details about the optimisation 
procedure and algorithms implemented in the AEROSTATE tool are described in [58-60]. 
This tool evaluates the lift-to-drag ratio at the reference design point by estimating the 
overall drag, as reported in Equation (1); nevertheless, in this case, the geometry of the 
lifting system is known, as it is computed at each iteration of the optimisation process, 
so the parasite drag coefficient of the wing components, Cpo'"8 (i.e., front/rear wing, 
vertical tip-wing, fin), is computed by integrating spanwise the airfoil drag contribution, as 
expressed in Equation (20): 





in i : 
Coo = — | Co gaily) e(y) dy (20) 
Sef —b 


where the function Cpyoi = f(Re, M, Crit) is calculated by Xfoil [61] and is provided as 
input to the procedure, S,,.¢ is the wing reference surface, and c(y) is the spanwise chord 
distribution. The Cpo'“S is computed by means of the component build-up method described 
in [53]; in particular: 


S 
Ca = OC (21) 
I Sref 





where the friction coefficient, Cy, the form factor, FF, and the interference factor, Q, are 
reported in [53]. Concerning the lift induced drag, once the design weight is defined as 
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an input, the optimiser uses a Vortex Lattice Solver to evaluate the vertical trim condition 
(Equation (13)), and so extracting the aircraft Cy, tj, and Cp; values. The solver used is the 
Vortex-Lattice-based AVL [62], and a typical representation of a box-wing configuration 
within this solver is reported in Figure 5. 


——— 
————SSS——— 
eee ee” 
—————— 
——————— 







































































Figure 5. AVL model for a box-wing configuration. 


As it is chosen to represent the fuselage as a flat plate in AVL, the vertical trim is 
evaluated in AVL with the fuselage always aligned to the flow; namely, the trim is solved 
with the angle of attack « = 0°, and by acting on the wings’ twist distributions, so as to not 
introduce the lift-induced drag coming from the simplified fuselage modelling. 

Additionally, in this second block of the design process, the maximum take-off weight 
of the aircraft is evaluated by means of the breakdown described by Equation (9); again, the 
payload and the fuel weight are calculated as described in Section 2.1, while the operating 
empty weight is subdivided as follows: 


Woe = wWinss 4 Wwrus 4 WweEns ae weys 4 Wwper (22) 


where W''"85 represents the weight of the lifting system, W'“’ is the fuselage weight, 
WEr8 is the engine weight, W°Y’ represents the weight of the on-board systems, and W0P? 
represents the weight of the aircraft operating items. As the geometry of each configuration 
designed in this phase is known, single contributions can be estimated by means of slightly 
more accurate methods than those used in the initialisation phase; in particular, the wing 
weight and the systems and operating weights are estimated by means of the method 
proposed in [63], and the fuselage weight and the engine weight are evaluated by means of 
the methodology used in [64]. During the optimisation, these estimations are also useful to 
evaluate the centre of gravity position, and so to evaluate the constraints concerning the 
pitch trim and the longitudinal stability. 

At the end of the process, it is possible to select an aerodynamic layout provided in 
output from the optimisation process, and to re-estimate the maximum take-off weight by 
means of Equations (9) and (22); if this value differs from the input provided by the starting 
initialisation within a prescribed tolerance, the procedure is iterated until the convergence 
is reached (Figure 6). 


Transonic Aerodynamic Assessment: Level 1 and 1/2 


If the design process involves transonic aircraft, such as commercial medium/long 
range airliners, it is necessary to provide some information about the transonic behaviour 
of the lifting system from the conceptual stage of the design. This aspect is relevant as, after 
a specific condition in transonic flight called drag rise, the aircraft performance deteriorates 
severely; in particular, for certain combinations of cruise Mach number, trim lift coefficient 
and wings geometry, a massive drag increase occurs; this is caused by the rise in shock 
waves (wave drag contribution), and by the interaction of these waves with the boundary 
layer, that thickens until local separation is reached, causing increases in pressure drag. 
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The drag rise condition needs to be avoided in the whole standard operating envelope of an 
aircraft. However, the transonic phenomena are very difficult to be detected and properly 
described by the low-fidelity methods that are typically used in the conceptual stages of 
the design process to provide fast and adequately reliable assessments. This because the 
transonic phenomena strictly depend on the local geometry of the wings, on the airfoil 
shape, on the joints and on the fillets’ design; moreover, transonic fluid physics is strongly 
three-dimensional and exhibits strong coupling and interaction between the boundary 
layer and the sonic phenomena such as shock waves. It is evident that potential methods, 
such as the VLM, cannot predict any of these physic phenomena, as they are not designed 
to take the fluid compressibility (and so the shock waves) and the fluid friction actions 
(and so the boundary layer) into account. In the literature, some simplified methods are 
provided that take transonic effect into account in the early stages of the design process: 
in [65], a quasi-3D model is presented, which couples a VLM solver with a strip model 
considering 2D airfoil transonic performance computed by means of a CFD solver; in [66], 
the wave drag coefficient is estimated by means of semi-empiric formulae provided by the 
Korn—Locke method. Although these methods are useful to identify relevant correlations 
between the aerodynamic performance and some macro-parameters (such as the airfoil 
shape, the wings’ sweep angles, the thickness distribution and the local lift coefficient), 
they may fail to provide reliable absolute performance prediction to be used with enough 
confidence in the conceptual design phase. In other words, with these models, we cannot 
ensure that a specific design coming from the Level 1 stage of the design process (Figure 2) 
is far from the drag rise conditions; this may compromise the development of the design, as 
it is possible to specifically evaluate the transonic performance only at Level 2 stage of the 
design process (Figure 2), where higher fidelity methods as CFD RANS solver are used. To 
avoid that, the design of configurations at Level 1 may present possible transonic issues at 
Level 2, and so to prevent restarting the design process from the beginning, a preliminary 
transonic assessment was performed for the specific design problem described in Section 3. 
In particular, a CFD campaign was carried out to improve the knowledge on the transonic 
behaviour of the box-wing lifting system, and to obtain information and data to calibrate 
the Level 1 aerodynamic optimisation set-up (design variables boundaries, aerodynamic 
constraints), with the aim of designing configurations that are confidently far from the drag 
rise condition; as this high-fidelity assessment was useful in properly tuning the Level 1 
design stage, this procedure can be used as a Level 1 and 1/2 block, parallel to the Level 1 
design block. 
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Starting geometry 


Aerodynamic 
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Update MTOW 
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Figure 6. Level 1 sizing procedure scheme. 
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Steady compressible RANS simulations were used to investigate the box-wing in 
transonic flight, and the performance correlation with some relevant parameters; in par- 
ticular, a focus on some macro parameters together with an analysis on some local effects 
have been provided, as described in [67]. The airfoil that was selected to be used for the 
box-wing is a supercritical profile and it is not considered as a design variable. Concerning 
the macro parameters, major effort was given to the role of the wing loading (Figure 7) and 
to the wings’ sweep angles; the effect of cruise Mach was investigated as well. 
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Figure 7. Wing loading variation effect on flow over the front wing (Mach contours @Mo = 0.79). 


Some local geometrical parameters were identified as more influential with respect 
to the transonic drag increase; in particular, the twist distribution in the tip zone of the 
front wing proved to be crucial to be properly set in order to avoid detrimental shock- 
induced drag increase (Figure 8 (left)). The load of the front wing tip region is enhanced 
by the aerodynamic behaviour of a positive swept wing; moreover, differently from a 
conventional wing in which the tip is free, for the box-wing, there are flow accelerations on 
the fillet with the vertical wing, which ease the flow to reach sonic conditions. The opposite 
happens for the rear wing, as the aerodynamic behaviour of the forward swept wing leads 
to a spanwise enhancement of the lift coefficient in the root zone; moreover, in the root 
zone, there is a significant straightening of the isobars, as the rear wing is not directly 
constrained to the fuselage: this effect reduces the aerodynamic sweep of this portion of 
wing. These two features of the forward swept wing can lead to transonic issues in the root 
zone (Figure 8 (left)). 

This qualitative information was then used to properly tune the Level 1 optimisation 
framework for the preliminary aerodynamic design. A conservative design approach 
was chosen for this stage; in particular, a set of constraints and boundaries for the de- 
sign variables was identified to allow the optimiser to design configurations that are 
confidently far from the detrimental transonic drag increase effects. For the design case 
described in Section 3, concerning a medium range transport box-wing aircraft, the bound- 
aries/constraints selected are summarised in Table 2. 
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Figure 8. Local transonic critical issues (left); box-wing without transonic issues (right). Mach 
contours referring to freestream Mo = 0.79. 


Table 2. Level 1 and 1/2 calibration of the preliminary aerodynamic optimisation. 














Front wing loading, Léont/Sfront <600 kg/m? 
Rear wing loading, Lrear/Syear <600 kg/ m? 
Front wing sweep angle, Agron; >35° 
Rear wing sweep angle, Ayear free 
Cruise Mach <0.79 
Front wing tip twist angle, Crone <—1° 
Rear wing root twist angle, ne <+1° 
Spanwise local lift coefficient, c)(y) <0.7 


2.3. High-Fidelity Performance Assessment 
2.3.1. Aerodynamic Performance 


In the third phase (Level 2, Figure 2), the aerodynamic layout provided by the op- 
timisation process was subject to refined analysis, concerning aerodynamics, structural 
sizing, weight estimation and performance assessment. First of all, the information gained 
through the design process allowed one to generate the 3D CAD model of the box-wing 
aircraft (Figure 9). 





Figure 9. CAD model of the reference PrandtlPlane. 


The designed configuration was available to be processed for detailed aerodynamic 
analysis. The high-fidelity aerodynamic assessment was performed by means of CFD 
simulations; drag polar curves were obtained by simulating the steady aerodynamics for 
different flight conditions, representative of the relevant operating conditions during the 
mission. The solver adopted was a steady compressible RANS with the k-e turbulence 
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model; the volume mesh was composed of a prism layer of 20 layers with a growth rate 
of 1.1, and tetrahedral cells outside; the volume mesh size typically adopted was about 
20 million cells for the half model, as a symmetry condition with respect to the longitudinal 
plane was applied. The software used to perform the computations was ANSYS Fluent [68]. 
In Figure 10, an example of the surface mesh is represented, and in Figure 11, an example 
of postprocessing of results is proposed. 





Figure 10. Surface mesh details (aircraft half-model). 
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Figure 11. Example of postprocessing. 


To reduce the computational cost, the component modelled was the wing—body 
assembly, while the drag of the vertical tails and engine nacelles were evaluated with the 
simplified models described in [53]. The solution chosen for the vertical tail assembly was 
a V-tail configuration; this has been shown to be necessary to alleviate potential flutter 
issues [69], as reported in [70-72]. This solution proved to be effective in meeting the 
flutter requirements in the case of the aircraft described in this article, as reported in [73]. 
The sizing of the V-tail was carried out with the aim of satisfying the directional stability 
requirements, as described in [74]; the flutter assessments were carried out subsequent to 
this sizing. 


2.3.2. Structural Design and Mass Estimation 


The mass breakdown was updated in this phase; in particular, the refinements regard 
the structural mass evaluation. As the simplified methods used to estimate the operating 
empty weight (and so also the structural weight) in the first two design stages are mainly 
referred to as statistical extrapolations concerning the structural data of classical aircraft 
architecture (tube-wing-tail), it was necessary to introduce specific physics-based models 
to evaluate the box-wing aircraft structural mass. The structural behaviour of a box-wing 
is very different from that of a conventional aircraft, concerning both the lifting system 
and fuselage structures. The main difference is related to the fact that the box-wing lifting 
system is over-constrained to the fuselage, as the main wing of a conventional aircraft is a 
statically determinate cantilevered lifting surface. The fuselage is a doubly supported beam 
in the case of a box-wing, borne in correspondence of the front and rear wing connections, 
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differently from the case of the conventional monoplane. To properly size and design 
the aircraft structural components, an in house-tool called WAGNER [75] was developed; 
the main goal of the WAGNER tool is to provide the preliminary structural design of 
the aircraft concerning the wings, the fuselage and the tailplanes’ structural components. 
The sizing procedure relies on an FEM-based solver, namely the ABAQUS commercial 
software [76]; the WAGNER tool provides the automatic generation of the 3D geometry of 
the aircraft structures and the mesh of the structural components, such as stringers, frames, 
ribs, pressurisation bulkheads, floor beams and struts (Figure 12); the load cases’ setting 
and computation through an ABAQUS solver (Figure 15); and the postprocessing of the 
results (Figure 16). 
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Figure 12. Typical Finite Element mesh of a PrP configuration (aircraft half-model) with some 
structural details. 


The WAGNER code follows the sizing scheme reported in Figure 13; the main input is 
represented by the geometry from the Level 1 block of the design process. 
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Figure 13. WAGNER structural sizing workflow. 


The first step of the structural design workflow of Figure 13 provides the fuselage’s 
structural sizing: the structure is modelled with shell elements for the skin, whereas 
stringers and frames are modelled with beam elements. The FE model of the fuselage is 
delimited by the end bulkheads, modelled as flat panels stiffened by beams in both radial 
(stringers) and circumferential (frame or ring) directions. Secondary structural parts of 
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the fuselage are considered as concentrated masses linked to the bulkhead structure (e.g., 
tail-cone) by means of multipoint constraints. Concentrated masses are also used to model 
landing gears and fuselage-mounted engines, connected to fuselage main frames. Finally, 
non-structural masses are spread over the mesh elements to account for the weight of pay- 
load, systems and operatings [63], as well as for doors, windows and special reinforcement 
structures [77,78]. 

Concerning the loading conditions for fuselage structural sizing, combined pressurisa- 
tion and inertial loads are considered; limit and/or ultimate pressure are superimposed to 
gravity load related to load factor, nz, ranging from —1 to +2.5. The sizing criterion is that 
the equivalent stress should be less than the allowable stress, with the latter being given 
by the yielding stress of the material divided by an appropriate safety factor, which usual 
standards require to be at least equal to 1.5. Figure 14 shows the global FE mesh of a half 
fuselage together with some example contour plots of stress and deformation induced by 
pressurisation and gravity loads. 


S, Mises 
Multiple section points 
(Avg: 75%) 











Figure 14. FE mesh of half a fuselage and example of contour plots of stress and deformation induced 
by ultimate pressurisation and gravity load with nz = +2.5 (deformed scale factor: 20). 


The second step concerns the assembly of the fuselage with the lifting system; the 
WAGNER software models the wings’ geometry and meshes the structural components; 
the arrangement of structural components, such as frames and stringers pitch, is provided 
in input by reference cases available in literature for transport box-wing [79,80]; referring 
to Figure 12, the front wing is fixed to the fuselage by means of surface-based constraints; 
in particular, the spar webs (front and rear) are fixed to the main frames of the fuselage, 
while the ribs of the front wing are fixed to the floor beams. The same type of connection is 
used to attach the fin structure to the fuselage. Then, the lifting system structural sizing is 
performed; in particular, the wings’ structural design is carried out by only considering 
static loads, according to a loading condition corresponding to a pull-up manoeuvre in 
cruise (M = 0.79, h = 11,000m, nz = 2.5, W = MTOW); qualitatively, two constraints are 
imposed into the sizing procedure: (1) the wing-tip deflection, intended as the maximum 
deflection occurring in one of the main wings, should be such that large displacements 
are avoided; (2) the structure must support the limit load (given as a combination of 
aerodynamic and gravity loads) within the elastic field, so the equivalent Von Mises stress 
in both wings is constrained to be lower than the yielding stress of the material (divided 
by a safety factor) [81]. The aerodynamic loads, such as the spanwise lift distribution on 
the lifting surfaces, are provided as inputs from the AVL code. The wing structural sizing 
process is integrated into an optimisation procedure; the objective function is the structural 
mass of the lifting systems, W“""8 s'™ucture (Equation (23)); the constraints are those related 
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to structural stiffness, expressed as a maximum limitation on wings’ tip displacement, 
gi? (Equation (24)), and structural strength, expressed as maximum limitation on the 
equivalent stress, a (Equation (25)), for the i-th wing; the design variable is the vector t of 
the thicknesses of the k-th structural component of the i-th wing t;, (Equation (26)). 


e Objective function 


N ’ 
min()~ wens structure ( t) ) (23) 
=] 
e Stiffness constraint 
SP) < oP = Ip” 24 
i (t) = Ax. tip 5 ( ) 
e Strength constraint 
a! (t) < omax/ksp (25) 
e Design variables 
=f (26) 
e Design space 
tmin< tik < tmax (27) 


In Equation (24), the maximum wing tip displacement, ee is defined as a fraction, 
Kip, of the half-wingspan, b/2; ktjy is considered equal to 0.1 for all the analyses carried 
out. In Equation (25), o\y,x is fixed as equal to the yielding stress of the material, and 
Ksp is the safety factor; ksp is set to 1.5. The wings are divided in a primary structure, Le., 
the structure in which all the loads are concentrated, namely the wing-box consisting of 
upper/lower stiffened panels, spars, ribs, and in a secondary structure, mainly composed 
of components related to fixed leading/trailing structure and movables installation. The 
FEM-based sizing only involves the primary wing structure, as the secondary structures 
are evaluated by the method presented in [82] and are modelled in the FEM as distributed 
mass properly placed spanwise (Figures 15 and 16). 





Figure 15. Secondary structures of the box-wing lifting system and aerodynamic forces (aircraft 
half-model). 
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S, Mises 
SPOS, (fraction = 1.0) 
(Avg: 75%) 





Figure 16. Example of deformed shape of the PrP configuration and Von Mises stress on both wings 
(aircraft half-model). 


The output of the procedure is a complete structural mass breakdown of the main 
structures of the aircraft; the evaluation of non-structural components, i.e., cabin furnish- 
ings, on-board systems, etc., is performed in WAGNER by using the model described 
in [63]. The final output of the WAGNER tool is represented by the total mass of the 
configuration and the related components’ breakdown, the coordinates of centre of gravity 
and the components of the inertia tensor. 


2.3.3. Mission Simulation 


The information collected in the Level 2 stage, concerning the refinement of the weight 
estimation and the aerodynamic performance assessment, are useful to integrate in an 
in-house-developed mission simulation tool; this simulator allows one to estimate rele- 
vant output regarding the aircraft mission performance; in particular, the main outputs 
achievable from the simulation are the range and the related fuel burnt, for every mission 
inside the payload-range envelope of the aircraft. The mission performance is computed by 
integrating the equation of motion of the aircraft [55,83]; the reference mission is divided 
into taxing, take-off, climb, cruise, descent and diversion (Figure 17). 


Ruler integration 
t+ At 







Simulation output: 
Flight parameters (Ah, Ax, V, etc.) 
Forces update (W, D, L, T) 
Fuel flow 






Mission 


ee ee ee 





Cruise : 7, 
Diversion 


Landing 


Figure 17. Mission simulation scheme. 
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Each mission segment is subdivided in timesteps At; the equations of motion are inte- 
grated by using the Euler algorithm, as reported in Equation (28) for a generic parameter k: 


k(t + At) = k(t) +k (£)At (28) 


The inputs for the simulation process are, for each time step: the aircraft’s aerody- 
namic performance, expressed as response surfaces of the drag as a function of lift and 
flight condition (Cp = f(C,, M), derived by the CFD simulations results); the initial mass 
breakdown, depending on the payload and on the initial fuel mass; the flight programme 
for each mission segment; the main performance of the propulsion system, namely the 
thrust specific fuel consumption, TSFC. The simulation models used for each mission stage 
are briefly described in the following: 


e The taxing fuel consumption is extrapolated by the data reported in [84]. The take- 
off phase was simulated by integrating the equation of motion of the aircraft in the 
longitudinal plane, also considering its pitch dynamics; the take-off simulation and 
analysis procedure used in this work is widely described in [85]; 

e The climb phase was simulated by integrating the equation of motion of the aircraft 
considered as a point mass in the longitudinal plane: 





tT — D 
tany = (29) 
L 
x = Vcosy (30) 
z = Vsiny (31) 
W = —TSFEC T (32) 


The climb programme was set as suggested in [86], namely by dividing the climb 
into segments flown at constant [AS or Mach number; the reference climb programme 
selected is “250 kt/300 kt/Moruise”. According to this programme: (1) the aircraft flies at 
IAS equal to 250 kt until it reaches an altitude of 10,000 ft; (2) the aircraft accelerates in an 
almost level flight until the IAS is 300 kt, and then flies at this speed until the crossover 
altitude (the crossover is defined as the altitude where the current Mach reaches the target 
cruise Mach); (3) the aircraft flies at a constant Mach number until it reaches the cruise 
altitude. During the climb, the aircraft changes its altitude, so its TAS increases as depicted 
in Figure 18 (left); an example of the aircraft climb trajectory is depicted in Figure 18 (right). 
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Figure 18. Reference climb programme. 
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e The cruise phase was simulated by integrating the equations of steady and level flights 
for the aircraft point mass model (Equations (33)—(35)); a constant altitude (z ) = 0 and 
constant speed (x = V = constant) flight programme was considered; the simulation, 
taking aircraft aerodynamic performance and cruise length into account, provided 
stepped cruise programmes if it resulted in performance gains. 


L = Weg (33) 
T=D (34) 
W = —TSFEC T (35) 


e The descent starts at the cruise altitude and ends at an altitude of 1500 ft; the equations 
of motion are obtained for the climb, and in the same manner, the flight programme for 
the descent is made by segments at constant [AS or Mach number [85]; the reference 
selected programme is “Meryise/300 kt/250 kt”, namely: (1) the aircraft flies at a 
constant Mach number from the cruise altitude to the crossover altitude; (2) the 
aircraft flies at IAS = 300 kt from the crossover altitude to an altitude of 10,000 ft; 
(3) the aircraft decelerates in an almost level flight until the IAS is 250 kt, and then flies 
at this IAS until an altitude of 1500 ft is reached. 

e Concerning diversion and loiter, analogous considerations about climb, cruise and 
descent were implemented. 


3. Results of the Design Process 
3.1. Input Data 


The design methodology described in Section 2 was applied to the design case of a box- 
wing transport aircraft; in particular, the specific case of the PARSIFAL project [74,87-89] 
is presented. 

To initialise the whole design process, it is necessary to define the Top Level Aircraft 
Requirements (TLARs). For the case study, the main TLARs are summarised in Table 3. 


Table 3. PrandtlPlane main TLARs. 








Max n° of passengers 310 
Design range 5000 km 
Cruise Mach 0.79 
Initial Cruise Altitude 11,000 m 
Max Wingspan 36 m 


The selection of TLARs for the PARSIFAL case is widely described in [90]; in particular, 
the aim of introducing the PrandtlPlane as an alternative to conventional tube-and-wing 
aircraft follows two main design drivers, provided by the forecasts for the commercial air 
traffic scenario of the future: 


e A huge increase in air traffic, especially on short/medium routes (up to 5000 km) [8] 
is expected; 
e At the same time, relevant problems of airport saturation are forecasted [11]. 


To address these contrasting objectives, the expected higher lifting capability of a 
box-wing architecture is exploited; in particular, the aim of this design study is to develop 
an aircraft with an increased payload (up to 310 passengers), but the with same overall 
dimensions (compatibly with airport apron constraints) of the conventional competitor 
currently operating in the sector of the medium routes (e.g., Airbus A320 family or Boeing 
737 family). The maximum wingspan is limited to 36 metres, compatibly with the apron 
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constraint provided by the ICAO Aerodrome Reference Code “C” [91]. Moreover, as the 
main objective for the commercial aviation of the future is to drastically cut the noxious 
emissions, the result of this design will be satisfactory only if a significant pollutant emission 
reduction is provided; this is theoretically possible by exploiting the expected better 
aerodynamic performance of the box-wing system. To properly assess the performance 
comparison with a conventional aircraft, it is necessary to identify a proper benchmark; 
the conventional benchmark selected is the CeRAS CSR-O01 [92], a short/medium-range 
aircraft developed in the framework of the project CeRAS (“Central Reference Aircraft data 
System”, [93]); the CeRAS project provides an open-access platform in which designs, 
geometries, data and procedures relevant to the CSR-01 configuration (CeRAS Short Range 
version 01, Figure 19) are collected. The CeRAS CSR-01 is a virtual simulacrum of an aircraft 
similar to the Airbus A320 [94] in terms of key features (Table 4), dimensions, payload 
and utilisation. 





Figure 19. CeRAS CSR-01 configuration. 
Table 4. Main features of the CeRAS CSR-01. 





Ref. wing area 122.4 (+32.2 *) m2 
Design range 5000 km. 
Max n° pax 186 
Wingspan 34.1 m 
Fuselage length 37.5 m. 
Cruise altitude 11,000 m 





Cruise Mach 0.79 


* Horizontal tail. 


3.2. Conceptual Design and Reference Layout Selection 


As the TLARs for the box-wing design have been defined, the first level of the design 
process is performed, providing the fuselage design and the general information to ini- 
tialise the aerodynamic optimisation in the Level 1 block of the design process (Figure 2). 
Concerning the fuselage design, the focus is mainly on the external shape, the cabin section 
layout and the internal arrangement, as the main design driver is represented by the target 
number of passengers. For the PARSIFAL case study, the fuselage of the box-wing aircraft 
was selected with the following main features: a double-aisle cabin section was selected, 
with a 2-4-2 seats layout (Figure 20(left)); the fuselage length is set equal to 44.2 m, similar 
to the longest fuselage of a competitor operating in the medium range sector (Airbus A321); 
the cabin section external shape is near-elliptical, following the studies provided in [95]. 
This fuselage, arranged in a high-density layout, is able to transport a number of passengers 
equal to 308 (Figure 20 (right)). 

Previous designs of wide body fuselage for box-wing aircraft were proposed 
in [40,96,97]; however, in the case of PARSIFAL, some actions to prevent penalisations 
in terms of turnaround time with respect to conventional competitors were provided, such 
as the introduction of enlarged aisles and of an extra boarding/de-boarding path; both 
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these aspects are highlighted in Figure 21, and their effects on turnaround time are widely 
discussed in [98,99]. In particular, the data provided by the simulations described in [98,99] 
show that the introduction of the PrP in the same sector as the A320-type aircraft implies a 
slight increase in turnaround time (+11% outstation, +25% full service) against a relevant 
increase in payload capacity (+66% passengers, +71% containers); these results were ob- 
tained by considering the same apron space (ICAO Code “C”) for the two aircraft. The 
abovementioned simulations take all the main aspects of airport operations into account, 
from passenger arrival at the terminal to aircraft take-off. 
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Figure 21. Cabin layout comparison between the PARSIFAL solution and a conventional single- 
aisle aircraft. 


Once the fuselage is defined, the next step required to initialise the design process 
(Level 0) is to identify the main characteristics of the lifting system to apply simplified 
models for the estimation of the aerodynamic performance; to define an initial study case, 
a wing loading similarity was set between the box-wing and the CeRAS: in particular, to 
define the initial lifting system surface, the same wing-loading of the CeRAS was selected 
for the box-wing. By using this assumption in the relations of Equations (1), (4) and (6) for 
the initial aerodynamic evaluations and Equation (9) for the rough weight estimation, as 
described in Section 2.1, the data were obtained as input to initialise the second stage of 
the design process (Level 1). In this phase, several layouts of the box-wing aircraft were 
evaluated by means of the aerodynamic optimisation described in Section 2.2; at the end of 
the process, the output selected for the final performance analysis is represented by the 
configuration reported in Figure 22, whose main features are described in Table 5. 

As detailed aspects related to the aeromechanical analysis of the box-wing aircraft 
may be out of the context of this paper, it is worth noting that such aspects have always 
been taken into account during the development of the aircraft. Specifically, general 
aspects of aeromechanical behaviour in the longitudinal plane of box-wing aircraft are 
extensively reported in [87]; a detailed summary of the design of the lifting system, with 
a focus on stability and controllability requirements, is given in [56]. The methods used 
to assess the aeromechanical characteristics of box-wing aircraft are extensively reported 


Aerospace 2021, 5, 292 


21 of 37 


in [87]; the approach used in the study involves a multi-fidelity approach, starting from the 
use of literature models [100], passing through extensive analysis campaigns with VLM 
solvers [58,74], up to the verification of requirements with CFD [74,88]. 





y [m] 


Figure 22. Reference box-wing layout. 


Table 5. Main features of the reference box-wing. 

















@ Design Point 
M079 
L/D 21.62 
Ch 0.4473 
Cp 0.02068 
(L/S) front 604 kg/m? 
(L/S) rear 477 kg/m? 
(L/S) rear L/S) front 0.789 
SSM 0.10 


It is worth pointing out the result concerning the wing loadings ratio (L/S)rear/(L/S) front 
reported in Table 5: this value is in agreement with what is reported in [58]: the parameter 
(L/S)rearL/S) front is the key lever to provide the maximum aerodynamic efficiency for box- 
wing aircraft compatibly with the constraints of static longitudinal stability and pitch trim; 
the result is also in agreement with what is described in [101], i.e., it provides the maximum 
performance in terms of Crying in clean conditions (unflapped) for box-wing lifting systems. 


3.3. Mission Performance Analysis 


The selected configuration was analysed through the Level 2 block of the workflow 
of Figure 2, by evaluating the aerodynamic behaviour and the structural mass in a more 
refined way to properly analyse the mission performance. 

The CFD simulation model presented in Section 2.3.1 was used to evaluate the box- 
wing aircraft aerodynamics in different flight conditions; in Figure 23, the wing—body 
aerodynamic efficiency curves (namely E = lift-to-drag ratio, L/D) versus the lift force 
generated are reported for four different flight conditions. 
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Figure 23. PrandtlPlane aerodynamic efficiency for different flight conditions. 

Details of the aerodynamic performance in the cruise condition (M = 0.79 and h = 11,000 m) 
is reported in Figure 24; in particular, Figure 24 (left) reports the cruise aerodynamic 
efficiency for the wing—body assembly (CFD results) and for the whole configuration, 
including the estimation of the parasite drag of vertical tail planes and nacelles ((53] 
model); in this graph, some reference points are reported corresponding to the required lift 
to trim the aircraft in some reference conditions, to highlight the aerodynamic efficiency 
envelope of this configuration varying weight through the mission (MTOW = maximum 
take-off weight; ZFW = zero fuel weight, WOE = operating empty weight, LF = cabin 
load factor). In Figure 24 (centre) the Cp-C;, polar curve is reported, together with the 
breakdown of the Cp contribution of each aircraft component; Figure 24 (right) reports the 
viscous pressure breakdown of the polar curve. 
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Figure 24. Details of the aerodynamic performance for PrandtlPlane in cruise condition; lift-to-drag 
curves (left), Cp breakdown -aircraft components- (centre), Cp breakdown -viscous/pressure- (right). 


The second main input to the mission simulation platform is the mass breakdown 
of the aircraft; in the Level 2 stage, the structural sizing has been performed through the 
FEM-based design and optimisation tool (Section 2.3.2); Table 6 reports the structural mass 
breakdown for the reference box-wing configuration. 
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Table 6. PrandtlPlane structural mass breakdown. 











Component Mass (kg) 
Front Wing 7166 
Rear Wing 6614 
Vertical Tip-Wing 460 
Fuselage 11,230 


Vertical Tail Plane 1026 


By calculating the other mass components as described in Section 2.3.2, the total mass 
evaluation of the reference box-wing conduces to a MTOW = 12,5126 kg, and a percentage 
repartition in components mass as reported in Figure 25; following the [63] model, the 
Systems category includes hydraulic, electric, pneumatic, fuel, engine, air conditioning, 
anti-ice, flight control, load systems, APU, avionics and instruments; Operatings includes 
cabin furnishings, cabin crew, crew seats, passengers seats, catering, emergency equipment, 
toilet fluid, oil and documentation. 
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Figure 25. PrandtlPlane mass breakdown: general (left), components detail (right). 


The information on the aerodynamic performance and weight of the aircraft are neces- 
sary to perform the mission simulations by means of the model described in Section 2.3.3. 
The simulations allow one to compute the mission performance of the reference box-wing; 
in particular, fuel burnt, range, flight time and other mission parameters such as time 
evolution of instant fuel consumption, speed, height, weight, etc., can be extracted from 
the simulation. Considering the box-wing harmonic mission, namely the mission with the 
maximum range at the maximum payload, the main results obtained are summarised in 
Table 7, and the main mission parameters are reported in Figure 26. 


Table 7. Box-wing harmonic mission performance. 











Number of passengers 308 
Mission range 5722 km 
Mission time 415 min 
Mission fuel 21,844 kg 
Total fuel 26,937 kg 


Mission fuel per pax/km 0.01239 kg/km pax 
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Figure 26. Box-wing harmonic mission height, aerodynamic efficiency and fuel flow evolution. 


3.4. Box-Wing Performance Comparison with Respect to the Conventional Benchmark 


The most relevant outcome of this work is the performance comparison between the 
unconventional PrandtlPlane aircraft and a conventional tube-and-wing benchmark. As 
reported in Section 3.1, the benchmark monoplane selected is the CeRAS CSR-01, a virtual 
mock-up of an Airbus A320-like aircraft, thus operating in the same routes of the reference 
box-wing; a top/front view comparison of the two aircraft is reported in Figure 27. 





Figure 27. Views comparison between the PARSIFAL PrandtlPlane and the CeRAS conventional aircraft. 


As far as the multilevel aircraft design and analysis diagram of Figure 2 is concerned, 
the activity on the CeRAS CSR-01 was limited to the ANALYSIS block, as the design and 
the geometrical characteristics of the aircraft are provided by the CeRAS database. Indeed, 
to make fair comparison between the box-wing aircraft and the benchmark monoplane, 
it is necessary to perform the aerodynamic, structural and flight analysis assessments 
with the same methods, tools and degree of fidelity. So, the same CFD analyses, the FEM 
based-structural sizing, the mass breakdown estimation and the mission simulation carried 
out for the box-wing aircraft were also performed for the CeRAS CSR-01. Figure 28 reports 
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the comparison of the aerodynamic performance of the two configurations: the wing—body 
aerodynamic efficiency (lift-to-drag ratio) curve trend with respect to the generated lift is 
reported in four flight conditions. Focusing on the cruise condition, namely at M = 0.79 and 
h = 11,000 m, it emerges that the box-wing aircraft is able to generate a larger amount of lift 
to trim larger weights, with a relevant increment in aerodynamic efficiency. Considering 
the point of maximum aerodynamic efficiency, the box-wing aircraft exhibits an increase 
in trim lift of +66% and a related increase in aerodynamic efficiency equal to +24%, while 
satisfying the same constraint on the maximum wingspan of the monoplane competitor. 
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Figure 28. Comparison of aerodynamic efficiency curves: PrandtlPlane vs. CeRAS. 


Figure 29 shows the detail of the operating points (from MTOW to Woe) reported on the 
aerodynamic efficiency curves for both the configurations; the dashed curves represent the 
result corrected with the additional drag of fins and nacelles calculated by the approximate 
method [53]. 
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Figure 29. Detail of cruise aerodynamic efficiency for the PrandtlPlane and the CeRAS. 


Figure 30 and Table 8 report the relevant data of the mass breakdown for the two 
aircraft; it is worth underlining that the fraction of the fuselage and lifting systems mass, 
together with the overall operating empty mass fractions, are very similar for the two 
configurations; the box-wing architecture does not introduce any significant penalisation 
in terms of structural mass with respect the conventional tube-and-wing solution. 
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Figure 30. Components mass breakdown comparison. 


Table 8. Mass breakdown comparison. 


PrandtlPlane MTOW = 125126 kg 








CeRAS CSR-01 PrandtlPlane 
Woe (kg) 42.054 68,866 
Woe /MTOW 54.7% 55.0% 
Wruel (kg) 17,100 27,000 
Wruel /MTOW 22.3% 21.6% 
Wray (kg) 17,670 29,260 
Wray /MTOW 23.0% 23.4% 
MTOW (kg) 76,824 12,5126 


The information on aerodynamic performance and weights were used as input for the 
simulation of the harmonic mission of the CeRAS CSR-01 configuration; Table 9 reports the 
main data relevant to the CeRAS harmonic mission together with the mission data for the 
same range for the PrandtlPlane. 


Table 9. Main mission data comparison. 








CeRAS CSR-01 PrandtlPlane 
Number of passengers 186 308 
Mission range 4790 km 4790 km 
Mission fuel 13,670 kg 18,108 kg 
Mission fuel per pax/km 0.01537 kg/km pax 0.01227 kg/km pax 


The main outcome of this comparison is that the PrandtlPlane solution allows one 
to reduce the fuel consumption per passenger of a quantity equal to 20% with respect 
to CeRAS, if the CeRAS harmonic range is considered. It is worth stressing the fact that 
a conventional aircraft needs a larger wingspan to satisfy increases in MTOW as that 
exhibited by the PrP. On the other hand, the PrP, while having a higher payload capability 
and consequently a larger MTOW, has an advantage in fuel consumption per passenger— 
kilometre, although it does not demand more apron space than the conventional plane 
(CeRAS CSR-01). 

Some relevant outputs of this mission simulation comparison, such as the flight 
trajectory, the mission aerodynamic efficiency and the instant fuel flow, are reported in 
Figure 31. 
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Figure 31. Comparison of flight trajectory, aerodynamic efficiency and fuel flow. 


The same comparisons were performed by considering different missions, by varying 
the initial condition in terms of payload and target mission range. In particular, several 
missions were simulated by considering a variation of the mission range in the interval of 
[500 km—4790 km] (the upper bound is fixed equal to the CeRAS harmonic range) and by 
varying the cabin load factor (namely the number of boarded passengers divided by the 
maximum number of passengers) in the interval [0.5-1]. Figure 32 shows the comparison, 
in terms of percentage variation in the averaged mission aerodynamic efficiency (namely, 
lift-to-drag ratio); the gains for the PrandtlPlane with respect to the CeRAS are relevant 
in the whole considered envelope, as the increase in mean aerodynamic efficiency ranges 
from 15% up to 22%. 
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Figure 32. Comparison in terms of mean mission aerodynamic efficiency. 


3.5. Box-Wing Operating Performance 

In this section, a comparison between the PARSIFAL PrandtlPlane and the CeRAS 
CSR-01 conventional aircraft is presented in terms of a payload-range diagram. The 
payload-range diagram and the mission performance inside this envelope were calculated 
with the mission simulation tool described in Section 2.3.3. In Figure 33, the payload-range 
diagrams of the two aircraft are superimposed: it emerges that the PARSIFAL box-wing 
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aircraft covers a wider market space, in terms of the maximum number of passengers 
(+66%) and maximum achievable ranges, despite having the same wingspan constraints 
of the CeRAS CSR-01 aircraft (compatibility with ICAO Reference Aerodrome Code “C”), 
and so operating from the same airports and aprons. 
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Figure 33. Comparison between the payload-range diagrams (PARSIFAL PrandtlPlane vs. CeRAS 
CSR-01). 


In the PARSIFAL payload-range diagram reported in Figure 34, different and flexible 
utilisations of the box-wing aircraft are qualitatively described by highlighting the different 
covered market sectors; in particular, the introduction of the PrandtlPlane allows one to 
open a new market sector related to the short /medium-haul with an increased payload (the 
‘PARSIFAL area’); up today, this market sector is not covered by any conventional aircraft 
compliant with the ICAO Reference Code “C” wingspan constraint, and the upper limit 
to the maximum payload for this interval of mission lengths is represented by the Airbus 
A321. Moreover, the same box-wing aircraft satisfies the requirements of the so-called 
Middle of the Market sector [102] by transporting about 200 (250) passengers for 8900 (7400) 
kilometres. In addition, it is possible to transport a similar number of passengers of the 
CeRAS CSR-O01 (configured in a high-density layout) for long range routes, up to 9350 km. 
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Figure 34. PARSIFAL PrandtlPlane payload-range diagram. 
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3.6. Discussion of the Performance Comparison between the Box-Wing and the 
Conventional Competitor 

The information presented in Sections 3.4 and 3.5 can be combined in order to compare 
the two reference aircraft in terms of utilisation and mission performance. In particular, the 
superimposition of the contour map of the percentage variation of the mission fuel burnt 
per pax km with the payload-range diagram is shown in Figure 35. Keeping in mind the 
difference in terms of payload capacity between the two configurations (as described in 
Figure 33), in the following discussion, the presentation of the performance comparisons 
refers to the cabin load factor range graphs, to simplify the representation of the data 
presented, as reported in Figure 35. 

The combined analysis of Figures 33 and 35 provides a complete overview of the 
performance comparison; it shows that: 


e The PARSIFAL PrandtlPlane has a larger pax-range envelope with respect to the 
CeRAS CSR-01 monoplane; in particular, at the harmonic point, the PrandtlPlane 
presents +66% more passengers and +19% longer range. Both the aircraft are compliant 
with the ICAO Aerodrome Reference Code “C” constraint (max wingspan equal to 
36 m); 

e The PARSIFAL PrandtlPlane can transport the same maximum number of passengers 
of the CeRAS CSR-01 (186 pax in high density) for about 9350 km, 95% more than the 
reference aircraft. 

e The PARSIPAL PrandtlPlane exhibits a gain in terms of mission fuel per pax km in the 
relevant area of the pax-range diagram, up to the harmonic range of the CeRAS CSR-01; 
considering the harmonic ranges, the PrandtlPlane needs 19% less fuel per passenger-— 
kilometre; the reduction in fuel per passenger is relevant, from —13% up to —22%, 
in the whole operating space considered; this also reflects the aircraft environmental 
performance: the introduction of the PrandtlPlane configuration allows a reduction in 
pollutant and greenhouse gas emissions per passenger, as widely discussed in [99,103]; 
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Figure 35. Comparison between the cabin load factor range diagrams (PARSIFAL PrandtlPlane vs. 
CeRAS CSR-01) and fuel consumption per pax km. 


e The diagram in Figure 36 shows the contour maps of the percentage difference of 
the fuel consumption per passenger—kilometre for the two aircraft, both for the zone 
inside the CeRAS CSR-01 envelope and outside this limit, up to the CeRAS CSR-01 
ferry range. It is clear that direct comparisons (i.e., with the same cabin load factor 
and range) can only be made within the limits of the CeRAS envelope; beyond this 
limit, the PARSIFAL PrandtlPlane can fly longer distances with the same cabin load 
factor, or have higher cabin load factors for the same range, with respect to the CeRAS 
competitor. In this area, the comparison in terms of fuel burnt per passenger—kilometre 
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cannot be made considering the same range and cabin load factor for the two aircraft; 
thus, the values obtained for PARSIFAL are compared with those relevant to the best 
performance of CeRAS, i.e., the missions at the border of the envelope (maximum 
cabin load factor for each considered range). As a result, the comparisons are carried 
out considering a same range for the two aircraft, but with different cabin load factors. 
In this zone of the diagram, the higher fuel efficiency of the PARSIFAL PrandtlPlane is 
combined with the capability to fly with higher cabin load factors for the considered 
ranges, and therefore the reduction in fuel consumption per passenger—kilometre 
increases very sensitively as the range increases, as shown in Figure 36. 
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Figure 36. Overview of the comparison (PARSIFAL PrandtlPlane vs. CeRAS CSR-01) of the fuel per 
passenger-kilometre and cabin load factor range diagrams, up to the CeRAS ferry range. 
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In the area of the cabin load factor range diagram beyond the ferry range of the 
CeRAS CSR-01 configuration, it is not possible to make comparisons in terms of fuel 
consumption. This region is highlighted in amaranth in the diagram of Figure 37. 
The PARSIFAL PrandtlPlane, with the same constraints on maximum wingspan of 
the reference monoplane competitor, is able to fly longer routes with a number of 
passengers comparable to the CeRAS CSR-01, thus offering an additional advantage 
in terms of operational flexibility. 
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Figure 37. Overview of the comparison (PARSIFAL PrandtlPlane vs. CeRAS CSR-01) of the fuel per 
passenger-kilometre and cabin load factor range diagrams, beyond the CeRAS ferry range. 
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4. Conclusions 


This paper described the potential effects of the introduction of PrandtlPlane aircraft in 
the short/medium-haul air traffic sector. The characteristics of this unconventional configu- 
ration, which is based on a box-wing lifting system, allow an increase in payload capability, 
without the need to increase the overall size of the aircraft compared to competitors in the 
same sector. Furthermore, the features of the lifting system enable the enhancement of the 
aerodynamic performance, in particular in terms of the lift-to-drag ratio, and therefore to 
achieve benefits in terms of reduction in fuel consumption per passenger. The introduction 
of this disruptive configuration is thus a possible solution to the main challenges of aviation 
in the near future, namely: the need to reduce pollutant and climate-changing emissions; 
the need to satisfy an expected growth in air traffic demand, especially on short/medium 
routes; and the need to alleviate problems of congestion and saturation of airport areas. 

Since the aerodynamic benefits may be reduced when a new configuration is assessed 
from the point of view of other disciplines, a multidisciplinary design and analysis method- 
ology was developed and used to develop a reference box-wing aircraft, and to analyse 
its operational potential and its performance. This methodology was structured with an 
increasing fidelity strategy, by using progressively more reliable design and analysis tools, 
with the aim of obtaining the best preliminary performance assessment. A conventional 
benchmark aircraft was identified to carry out the performance comparison; the selected 
aircraft was the CeRAS CSR-01, a virtual aircraft mock-up with characteristics similar to an 
Airbus A320. 

The analysis and performance comparison was carried out with the same procedures 
and tools for the two considered configurations. The results obtained proved that the 
PrandtlPlane aircraft is able to carry 66% more passengers than the conventional bench- 
mark, while maintaining the same dimensional constraints, hence keeping the possibility 
of operating from the same airport aprons. In addition, the PrandtlPlane is able to carry 
the same maximum payload as the CeRAS CSR-01 for up to 95% longer routes and, more 
importantly, it is able to reduce the fuel consumption per passenger over the whole enve- 
lope of the operational missions of interest. In particular, the maximum reduction in fuel 
consumed per passenger transported was quantified as 22%, when considering the typical 
operational missions of the CeRAS CSR-O01. Finally, this paper shows how the box-wing 
architecture confers, to aircraft with the typical dimensions of short-to-medium range, 
the capability to cover additional sectors of the air transport market, such as Middle of the 
Market and long range. 

The multidisciplinary design and analysis methodology presented in this paper can 
be applied to box-wing aircraft of various categories; further investigations will concern 
the development and the performance analysis of box-wing aircraft for the regional and 
commuter classes. In addition, the modularity and the flexibility of the design workflow 
allows the efficient integration of additional design elements, such as the introduction 
of hybrid-electric powertrains, or the inclusion of analysis and design methodologies of 
increasing fidelity, such as CFD-driven shape optimisation. Studies concerning the manu- 
facturability and assembly of the box-wing aircraft are also worthy of further investigation: 
a step forward has been made concerning the actual benefits deriving from the operability 
of the box-wing, but there is still uncertainty regarding the possibility of setting up effective 
assembly lines. Investigating this topic may represent a further step towards the actual 
introduction of box-wing aircraft in the transport aviation sector. 
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Abbreviations 


The following abbreviations have been used in the paper: 


APU Auxiliar Power Unit 

AR Aspect Ratio 

AVL Athena Vortex Lattice 

CAD Computer Aided Design 

CeRAS Central Reference Aircraft Data System 

CFD Computational Fluid Dynamics 

CSR CeRAS Short Range 

FEM Finite Element Method 

IAS Indicated Air Speed 

ICAO International Civil Aviation Organization 

LF Load Factor (passenger cabin) 

MTOW Maximum Take-Off Weight 

PrP PrandtlPlane 

RANS _ Reynolds Averaged Navier-Stokes equations 

SSM Static Stability Margin 

TLARs Top Level Aircraft Requirements 

TSFC Thrust Specific Fuel Consumption 

VLM Vortex Lattice Method 

Symbols 

The following symbols have been used in the paper: 

b Wingspan m 
C Chord m 
Cp Drag coefficient 

Cpo Parasite drag coefficient 

Cpj Induced drag coefficient 

Ciinane Wave drag coefficient 

Cp foil Airfoil drag coefficient 

Cp tot Total drag coefficient 

Cr Friction coefficient 

Cre Equivalent skin friction coefficient 

Gi. Lift coefficient 

C Section lift coefficient 

Cm Pitch moment coefficient 

d Diameter m 
D Drag N 
e Oswald factor 

E Aerodynamic efficiency (Lift to Drag ratio) 

FE Form factor 

g Inequality constraint 

g Gravity acceleration m/s* 
h Altitude m 
h/b Wings height to span ratio 

k Polar drag coefficient 


Ktip 


Stiffness constraint factor 
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Nz 

Q 

S exposed 
S ref 

S wet 


Other subscripts: 


comp 
cruise 
cyl 
eng 
front 
fus 
fuel 


max 
min 
oper 
pay 
rear 
root 
sys 

Wel 

tip 
trim 
vertical 
WINg 


Strength constraint safety factor 
Length 

Lift 

Lower boundary 

Lifting surface wing loading 
Mach number 

Vertical load factor 

Interference factor 

Planform area of the wing exposed to the flow 
Reference surface 

Wetted surface 

Vector of thicknesses of structural wingbox components 
Thickness to chord ratio 

Thrust 

Upper boundary 

Speed 

Weight 

Design weight 

Operating empty weight 
Design variables vector 

Aircraft longitudinal position 
Longitudinal leading edge coordinate 
Spanwise coordinate 

Aircraft vertical position 

Angle of attack 

Trajectory slope 

Wing tip displacement 
Tolerance 

Section twist 

Taper ratio 

Sweep angle 

Equivalent stress 


Component 
Cruise 

Cylinder 
Engine 

Front wing 
Fuselage 

Fuel 

Maximum 
Minimum 
Operating items 
Payload 

Rear wing 

Root section 

On board systems 
Threshold 

Tip section 
Trim condition 
Vertical tail 
Wing 


kg/m? 


MPa 
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Abstract: The feasibility of regional electric aviation to reduce environmental impact highly depends 
on technological advancements of energy storage techniques, available battery energy density, and 
high-power electric motor technologies. However, novel airframe technologies also strongly affect 
the feasibility of a regional electric aircraft. In this paper, the influence of novel technologies on 
the feasibility of regional electric aviation was investigated. Three game-changing technologies 
were applied to a novel all-electric regional aircraft: active flow control, active load alleviation, 
and novel materials and structure concepts. Initial conceptual design and mission analysis of the 
aircraft was performed using the aircraft design framework SUAVE, and the sensitivity of the most 
important technologies on the aircraft characteristics and performance were studied. Obtained results 
were compared against a reference ATR-72 aircraft. Results showed that an all-electric aircraft with 
airframe technologies might be designed with the maximum take-off weight increase of 50% starting 
from the battery pack energy density of 700 Wh/kg. The overall emission level of an all-electric 
aircraft with novel technologies is reduced by 81% compared to the ATR-72. On the other hand, 
novel technologies do not contribute to the reduction in Direct Operating Costs (DOC) starting from 
700 Wh/kg if compared to an all-electric aircraft without technologies. An increase in DOC ranges 
from 43% to 30% depending on the battery energy density which creates a significant market obstacle 
for such type of airplanes. In addition, the aircraft shows high levels of energy consumption which 
concerns its energy efficiency. Finally, the sensitivity of DOC to novel technologies and sensitivities 
of aircraft characteristics to each technology were assessed. 


Keywords: aircraft design; airframe technologies; aircraft sizing; all-electric aircraft; multi-disciplinary 
design optimization 


1. Introduction 


Significant climate changes and potential environmental impact due to increased 
transportation in the near future have motivated many industries to focus on reducing CO? 
and NOx emissions. As a major transportation method, the aviation industry also follows 
the trend to reduce the emission of new generations of aircraft. Improvements in airframe 
and engine technologies increase aircraft efficiency and reduce their emission. However, 
a potential increase in air transportation may still lead to an increase in overall CO and 
NOx emissions. Under Flightpath 2050 [1], The European Commission has set a future 
challenge for the new generation of aircraft to reduce their total emission. Figure 1 shows 
three schematic trends: if no advancements in aircraft technologies are present, if currently 
feasible technology advancements are achieved, and if novel technologies are introduced. 
This challenge leads to developing alternative environmentally-friendly energy sources 
as one of the main solutions for environmental impact reduction. Aircraft electrification 
has become one of the most popular approaches to reduce aircraft emissions. Today, many 
companies work in various directions to make electric flights available: improve battery 
energy capacity, modify and develop new propulsion systems, and introduce new aircraft 
configurations more applicable for future electrification. 


https://www.mdpi.com/journal/aerospace 


Aerospace 2021, 5, 163 


2 of 29 


Impact | 
CO 

(CO2) Frozen tech 

(no measures) 


SS ——————E 


: 
} es 
} 
} 
} 
i 
1 
i 
; 






Ref. _ 
100% 





2010 2020 2030 2040 2050 


Figure 1. Forecast of CO2 emission impact due to increased transportation and emission reduction 
goals [1]. 


At the moment, research is highly focused on hybrid- and all-electric propulsion 
systems and aircraft configurations to reduce the environmental impact produced by 
aircraft. Friedrich introduced a simulation technique for a single-seat hybrid-electric 
demonstrator, designed the aircraft, and performed a scaling analysis to determine a large 
fuel savings impact for the small- and mid-scale sector of aircraft [2]. Hamilton investigated 
the effect of hybrid-electric aircraft from the operational perspective to reach an optimal 
aircraft operation with battery energy density constraints [3]. A lot of work is also dedicated 
to the introduction of distributed electric propulsion (DEP) energy networks and aircraft 
concepts featuring DEP. Kim summarized major contributions towards the development of 
DEP aircraft and relative technologies [4]. Finger focused on sizing methodologies, and 
aircraft design of future general aviation aircraft that could be not only hybrid-electric 
but also fully electric [5-7]. De Vries also developed initial sizing methodologies for 
hybrid-electric aircraft for conceptual aircraft design but introduced effects of DEP at early 
conceptual design stages [8,9]. Pornet also introduced a sizing methodology for hybrid- 
electric aircraft and compared performance between the conventional reference aircraft and 
its hybrid version [10]. Sgueglia introduced an MDO framework for hybrid-electric aircraft 
with coupled derivatives and performed various optimization analyses to demonstrate 
the capabilities of the framework [11]. The work of Hepperle [12] addressed the potentials 
and limitations of all-electric aircraft specifically and introduced aircraft modifications that 
could potentially improve the aircraft performance at given battery energy densities. The 
ultimate goal of aircraft electrification is to achieve a fully-electric flight at adequate aircraft 
weight characteristics to maximize the emission reduction at the aircraft level. However, 
the availability of all-electric aircraft is limited to the General Aviation sector due to low 
battery energy density compared to the Jet-A fuel and relatively low maximum power 
capabilities of modern electric motors. Multiple all-electric concepts have already been 
introduced or are being developed at the moment. From existing airplanes, Pipistrel has 
already certified a twin-seat all-electric aircraft that has a maximum speed of 185 km/h 
and a maximum range of 139 km [13]. Bye Aerospace is currently performing flight tests 
on a one-seat general aviation aircraft that will be certified under the FAR Part 23 category. 
The eFlyer has a maximum cruise speed of 250 km/h and endurance of three hours [14]. In 
2016, Siemens tested an all-electric energy network by retrofitting the Extra 300 aerobatic 
aircraft and set two world records [15]. The Equation Aircraft has also demonstrated an 
all-electric twin-seat amphibious aircraft with the empennage-mounted engine that can fly 
up to 240 km/h and up to 200 km [16] Airbus has also developed several versions of an 
experimental all-electric E-fan aircraft to test the capabilities of an all-electric aircraft [17]. 
The E-Fan concept had multiple versions, from a technology demonstrator to a production 
variant. MagniX has introduced an all-electric propulsion system and first retrofitted De 
Havilland Beaver and then 208B Cessna Grand Caravan. The electrified Caravan became 
the largest all-electric airplane until today [18]. It must be noted that the list of existing 
all-electric aircraft companies is incomplete and more companies exist and are at different 
design stages. 
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To achieve all-electric regional aviation faster, not only the energy source and propul- 
sion system technologies must be improved, but also advanced airframe technologies 
need to be considered. Such novel technologies may significantly increase the aircraft 
performance characteristics and overall efficiency that will enable earlier integration of all- 
electric commercial aircraft into the market before. Liu introduced initial estimations of the 
impact of novel technologies for a range of aircraft from the short-range to the long-range 
[19]. Results showed that novel airframe technologies may significantly improve aircraft 
energy efficiency. 

Under the Excellence Cluster Se*A (Sustainable and Energy Efficient Aviation), three 
energy-efficient aircraft are to be designed to cover the majority of commercial aircraft 
operations. Figure 2 shows three sample new energy-efficient airplanes: a short-range 
propeller aircraft that shall be applicable to the ATR-72, the medium-range aircraft having 
similar mission requirements to the Airbus A320, and the long-range aircraft to cover 
ranges similar to the Boeing B777. 





Figure 2. A family of energy-efficient aircraft under the Se*A cluster. 


The initial cluster requirement for the short-range aircraft was to investigate the 
feasibility and availability of an all-electric short-range aircraft if novel airframe and 
propulsion technologies are introduced. A few novel technologies have been considered 
for the design of future regional aircraft. Hybrid laminar flow control (HLFC), load 
alleviation, and advanced materials and structure concepts. In addition, high-temperature 
superconducting motors and high-energy capacity batteries have been considered for the 
present research. An aircraft with similar top-level requirements as ATR-72 is designed 
considering the mentioned technologies and a battery-based full electric propulsion system. 

The present work is divided into multiple sections. Section 2 describes novel tech- 
nologies implemented in the new short-range aircraft. Section 3 describes methodologies, 
models, assumptions, and software used to perform novel technology assessments to de- 
sign the aircraft. Section 4 describes the aircraft’s top-level requirements (TLRs). Section 5 
describes the aircraft concept selection and its initial sizing. Section 6 presents multidis- 
ciplinary design optimization studies of the SE*A aircraft to refine the design and make 
important design decisions. Finally, Section 7 assesses the influence of novel technologies 
on aircraft availability and the strength of each technology impact. 


2. Novel Airframe and Propulsion Technologies 
2.1. Hybrid Laminar Flow Control 


The generation and extension of laminar flow over the aircraft surface significantly 
affect the overall aircraft drag, reducing aircraft weights, fuel burn, and operating costs. 
Preliminary estimations of an aircraft that features extended laminar flow along the wing, 
empennage, and fuselage demonstrated significant reduction (up to 50%) in overall drag 
and proved the importance of laminar flow control (LFC) [20]. Figure 3 shows the effects of 
the new wing design featuring the LFC and compares overall drag to the reference aircraft. 

However, natural laminar flow (NLF) is limited along the shape, so to extend it even 
more, hybrid laminar flow control (HLFC) is required: the laminar boundary layer is not 
only extended using NLF design but also using the active suction technology. In this 
technology, the air is sucked from the aircraft’s outer surface to delay the transition of 
the boundary layer and enable substantially higher percent laminarization compared to 
conventional surfaces. The skin of each aircraft component is split into two segments: 
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a porous sheet and an inner sheet that supports the outer sheet. The inner sheet has 
orifices that suck the air from the boundary layer and delay transition. In each chamber, 
an individual pressure is adjusted by the throttle orifices, so that the pressure difference 
between the outside and the chamber delivers the locally desired amount of mass flow 
through the surface. Figure 4 shows a schematic image of the skin layout and the HLFC 
system for a wing section. The applied technology in this project is based on [20,21], which 
describes numerical approaches with active laminar flow control and also describes current 
progress in this technology. 
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Figure 3. Effect of laminar flow control on overall aircraft drag [20]. 
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Figure 4. Schematic views of the active suction system [20]. 


2.2. Load Alleviation 


Load alleviation introduces various techniques to reduce the bending moment expe- 
rienced by the aircraft and that have a passive or active nature. Reduction of maximum 
bending moment enables the design of lighter wings for lower limit load factors, which 
will improve aircraft fuel efficiency. 

Passive load alleviation solutions consider nonlinear stiffness material design [22], 
viscoelastic damping design [23], new structural concepts [24,25], and local morphing 
structures [22]. Nonlinear stiffness materials may improve the load distribution on the 
wing under low load cases (0.5 g to 1.5 g) and improve performance efficiency at those load 
cases. Finally, both aeroelastic tailoring and local morphing structures aim to extend the 
aeroelastic design space. Morphing is considered into two scenarios: deliberate structural 
non-linearity that affects the wing deformation and reduces effective angle-of-attack of the 
wing section under the load and change in airfoil shape to achieve load reduction [22]. 

The wing active load alleviation uses different types of flow control over the wing to 
achieve a more favorable wing load distribution and reduce the wing bending moment. 
Previously, researchers have approached the design of active load alleviation systems 
in different ways. Rossow et al. [26] investigated the feasibility of an aircraft featuring 
load alleviation technologies; Sun et al. [27] looked at active load alleviation from the 
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control system perspective; and Zing et al. [28] demonstrated the experimental results 
of the high aspect-ratio wing wind-tunnel test with active load alleviation while using 
piezoelectric control. The work of Fezans et al. focused on novel sensors and control 
systems technologies to enable rapid and robust active load alleviation [29-32]. Finally, 
active load alleviation could also be reached by fluidic or micro-mechanical flow actuators 
that change the load distribution along the wing and reduce the bending moment [33]. 


2.3. New Materials and Structure Concepts 


Novel structural concepts and materials are being developed to improve the aircraft 
structure in terms of stiffness and weight. Bishara et al. [34] describe advanced structural 
design with the integration of active flow control, which is directly applicable to the 
presented research. Under the excellence cluster, the reduction in the airframe weight 
is assumed to be reached by application of CFRP thin ply laminates. In addition, new 
structural concepts must be applied to satisty HLFC requirements [22]. 


2.4. High-Power Superconducting Motors 


One of the major problems in hybrid- and all-electric commercial aircraft deals with 
limitations of maximum available electric motor power. The motor resistance increases 
rapidly with the increase of its power, which rapidly diminishes the engine efficiency. One 
solution to preserve the motor efficiency and even increase it is to use superconducting 
motors. Such motors can be synchronous, trapped flux, and fully superconducting. They 
may also operate at low and high temperatures. To maintain required temperatures inside 
the motor, a cryogenic cooling system is required. High-temperature superconducting 
(HTS) motors can achieve high power-to-weight ratios with high power densities without 
major weight losses compared to low-temperature superconducting motors. Ref. [35] 
provides a descriptive survey of HTS motors and their future trends. For the present 
research, a high-temperature, fully superconducting motor was considered. 


2.5. High Energy Capacity Batteries 


Currently, the gravimetric energy density of modern batteries generally does not 
exceed 300 Wh/kg [36]. Such energy densities may be suitable for ultra-light or limited 
general aviation all-electric aircraft. However, they are insufficient to enable not only 
all-electric aircraft of the size of ATR-42 and larger, but even the hybrid-electric version of 
such aircraft is infeasible at the moment since the energy capacities are around 50 times 
less compared to kerosene-based fuels. On the other hand, the research of high-capacity 
batteries makes progress. The latest laboratory research demonstrated battery energy 
capacities of 800 Wh/kg on the cell level and 500 Wh/kg on the pack level for military 
applications [36]. In addition, the National Academy of Sciences predicts battery pack 
energy density to achieve 400-600 Wh/kg at Technology Readiness Level (TRL) 6 in the next 
20 years (roughly, by 2035) and be commercially available in 30 years [37]. If the prediction 
becomes real and the linear trend remains, then batteries may reach TRL 6 with energy 
densities of 550-850 Wh/kg and be commercially available by 2060. To investigate potential 
scenarios for the far future and examine battery energy impact on the aircraft feasibility, it 
is assumed that the pack energy densities may range between 500 and 1100 Wh/kg. 


3. Implementation of Novel Technologies in Aircraft Analysis and Design 
3.1. Initial Aircraft Sizing 


The conceptual design was performed using various tools. OpenVSP [38] and CATIA 
were used for the aircraft geometric modeling, and SUAVE [39] was used for more defined 
aircraft sizing, performance, and mission analysis. To have more capabilities from the air- 
craft performance analysis standpoint, SUAVE was extended to have classical performance 
analysis methods used in [40-43]. 

The aircraft sizing within SUAVE is performed iteratively. First, initial geometric 
specifications based on TRLs and any available information are input into the SUAVE. 
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The information includes the initial wing planform, its section properties, flaps charac- 
teristics, empennage and fuselage geometries, the definition of the propulsion system, 
and all required components such as batteries, gearboxes, power management systems 
(PMAD), cables, etc. Finally, initial guesses of aircraft weight are input. Then, the first 
iteration of constraint analysis using equations provided by Gudmundsson [41] is per- 
formed. The constraint diagram block includes rapid estimation methods for gas-turbine, 
piston, and electric aircraft power lapse with altitude and has options to either estimate 
required aerodynamic properties for the constraint diagram or have them as fixed in- 
puts based on historical trends. Based on the design point selection, the wing loading 
and power-to-weight ratio are used to run the SUAVE mission analyses to estimate the 
aircraft performance and its required weight. Then, obtained maximum take-off mass 
(MTOM) is compared to the initial guess and updated if the tolerance is not reached. In 
addition, parameters such as the minimum drag coefficient (Cp jn), Oswald efficiency 
(ce), and maximum lift coefficient (Cina) for clean and flapped configurations are input 
into the constraint diagram again to update all constraint curves and run the loop again. 
When the tolerance is reached, the program moves to the aircraft performance block to 
obtain performance plots for the given aircraft. Particular care was taken to size the wing 
flaps to meet take-off and landing field length requirements. Methods of Torenbeek [42] 
and Roskam [43] were implemented within SUAVE to analyze various types of leading- 
and trailing-edge devices. The empennage sizing within SUAVE is based on the fixed tail 
volume ratio and then updated and corrected separately based on the desired CG envelope. 
Performance analyses within SUAVE included take-off, all engines operative (AEO), and 
one engine inoperative (OEI) climb, cruise, descent, and landing. Figure 5 shows the sizing 
process within SUAVE. 
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Figure 5. Initial aircraft sizing framework using SUAVE. 


3.2. Energy Network 


In the SUAVE analysis framework, a propeller electric energy network with HTS 
motors was implemented. Figure 6 shows the all-electric energy network layout. 
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Figure 6. Battery-electric energy network layout. 


Propeller modeling was performed using the cubic spline method described by Gud- 
mundsson [41]. Propeller thrust can be described with a cubic spline based on design 
operating conditions. The thrust is defined by: 


T= AV? +BV*+CV+D (1) 


where A, B,C, D are the coefficients required to fit the curve accurately. To obtain those 
coefficients, a system of four equations is established. The first equation is the equation 
for the static thrust. The second equation is derived for the cruise speed and the desired 
propeller efficiency. The third equation is the derivative of the curve at the cruise speed, 
which is equal to zero. Finally, the last equation is derived for the desired maximum cruise 
speed. The final system of equations for the thrust as a function of cruise speed becomes: 


0 0 O 1)/(A Te 
3 2 
V3. V2 Vy 1) (dD Ty 


Solving Equation (2), a general thrust expression is found, which then is used to find 
the propulsive efficiency of the propeller. Then, knowing efficiencies of components of 
the energy network shown in Figure 6, the total power at the battery is found to be used 
during the SUAVE mission analysis. 

The propeller was sized using a combination of propeller diameter trends provided 
by [44] and the sizing method described by Torenbeek [42] where the required propeller 
RPM depends on the tip Mach number was used to estimate the maximum propeller RPM. 

The assumption of the battery gravimetric energy density remains one of the most 
influential assumptions for aircraft design. Since the availability of a full-electric regional 
aircraft is challenging to achieve, an optimistic scenario of the energy density at TRL 6 
was assumed. This way, if the aircraft is feasible, the task will be to investigate more 
realistic bounds of minimum energy densities. Otherwise, the availability of the all-electric 
airplane is highly unlikely. Consequently, the average optimistic pack energy density (wz) 
of 700 Wh/kg was assumed. To ensure safe battery operation, the state-of-charge of 20% 
was assumed. 

High-temperature, fully superconducting motors were used for the airplane since 
the required power frequently reaches the order of Megawatts. All superconducting 
components also have cooling systems to enable desired operating temperatures. The 
cooling system power loss was estimated using the survey provided in [45] and assuming 
the Carnot efficiency (€) of 30%. Table 1 summarizes all assumptions considered for the 
gas turbine engine and an electric propulsion network. 
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Table 1. SE*A SR energy network assumptions. 
Component Parameter Value Reference 
High-power electric motor Efficiency ye = 99.5% [35] 
Motor cryo-cooling system Power loss Piosg = 10 KW [45] 
Carnot efficiency e=03 [45] 
PMAD Efficiency YpoMApD = 99% [46] 
W 
Battery Energy density Wp, = 700 ie [35] 
Electric cables Densit = 3.9 kg [46] 
- Pcable — m 
Efficiency Neable = 29-970 [46] 


3.3. Weights Estimation 


Weights estimation analysis was performed using the FLOPS [47] method programmed 
within SUAVE. To improve the accuracy of high aspect ratio wing weight estimation, 
FLOPS was coupled with the physics-based wing weight estimation tool EMWET [48]. The 
EMWET analysis was performed for a metal wing structure and was corrected by a factor 
based on the composite wing technology level. Loads used for the EMWET weight estima- 
tion were obtained within SUAVE using coupling with AVL [49] Vortex-Lattice method 
software. In the case where the aircraft wing span reaches Class-C airport regulation limits, 
an additional weight penalty function due to the wing folding was introduced, so the air- 
craft can maneuver at the airport and maximize its flight efficiency. The estimation method 
is based on [50] and accounts for the weight of the insert, the folding mechanism, and the 
pin to fix the wing. The wing weight corrected by the weight of the folding mechanism is 
defined by 


Id Ida | wi : _ pi 
mp = my 1 + mfsia + Mygid + Mora) (3) 


where mo is the mass of the wing with the folding mechanism, mons “'@ is the mass ratio 


without the folding mechanism, 7'"S is the mass ratio of the folding insert, mee is the 
mass ratio of the folding mechanism inside the insert, and mo) , is the mass ratio of the 


folding pin. The insert mass ratio was defined based on the curve-fitting of three separate 
formulations available from [50]: folding at 32%, 48%, and 64% of the span 


0.24(mpws — 0.924)* — 0.007, 32% span 
Mold = \ 0.45(titpws — 0.916)? — 0.02, 48% span (4) 
0.84(mpws —0.98)* — 0.107, 64% span 


where 7ipws is the mass ratio of the section outboard of the folding pivot point defined by 


(Cy ai C;) (0.5b — y)(tyCy =F trCy) 


mpws(y) = (GAGE p ©) 


where C is the chord length, b is the wing span, t is the wing thickness-to-chord ratio, and y 
is the spanwise location of the folding insert. Subscripts r, t, and y correspond to the wing 
root, tip, and the folding insert. 

Finally, folding mechanism and pin masses are defined by 


mye = —3.76(Mpws — 0.2)* + 0.08 (6) 
miner = —0.6(mMpws — 0.2)? +0.02 (7) 


Figure 7 shows the wing weight penalty as a function of the wing joint location along 
the semi-span. 
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Figure 7. Wing folding penalty as a function of the wing semi-span [50]. 


HTS motor weight estimation was based on relations defined by Stuckl [35] and 
account for both the motor and the cooling system weight. The motor power-to-weight 
ratio is defined by 


iP 
WwW 2.63(Paot)*”” (8) 


where Por is the motor power in kW and the power-to-weight ratio P/W is defined 
in kW/kg. 
The cryogenic cooling system weight is defined by 


Mcooling = 27.9 Pioss ° e~ 1225(l0810(Fioss)) (9) 


where Pj,., is the cooling system input power loss in kW and mMoo}ing is the cooling system 
weight in kg. 

Propeller weight was estimated using a relation described in [51] where the propeller 
weight depends on the number of blades, the activity factor, the motor power, propeller 
maximum RPM, and its material. A composite propeller with six blades was considered 
for the present design. 

The gearbox weight was estimated using relations provided by [52] and is defined by 

po-76 . RPM?:13 


gearbox motor 


Wn — 
gearbox 0.89 
RPM sropeller 


(10) 


where P is the power at the gearbox, K is the gearbox technology factor equal to 26 for 
future gearboxes, RPMinotor is the motor power assumed at 4800 RPM, the RPMnotor, and 
the Mearbox is the gearbox mass in kg. 


3.4. Novel Airframe Technologies Assumptions 


At the conceptual design stage, novel technologies have been integrated into the 
design process in terms of coefficients and correction factors. Laminar flow control is 
estimated by modifying the flow transition location along the wing profile. The desired 
goal within the SE*A for the flow transition is set to 70% for the wing if no disturbances 
such as the propeller wake or the location of other components is present. Otherwise, 
specific transition assumptions are prescribed. Wing folding also affects the application 
of the HLFC. The folding location has four general locations. If the wingspan is less than 
36 m (Class C airport limitation), then no folding is applied, and the wing features the 
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HLFEFC along the complete span. If the wingspan is longer than 36 m but not longer than 
110% of the Part C requirement, then the tip folding is applied, the wing has HLFC until 
the wing folding location, and the folding tips feature natural laminar flow (NLF). If the 
folding location occurs at any position which splits the aileron, then the folding position is 
automatically moved to the position where the flap ends and the aileron start, so the aileron 
design is significantly simplified. The portion of the wing after the folding joint does not 
feature HLFC and only has NLF. Finally, if the wingspan is such that the aileron can be 
located at the folding portion of the wing, the folding joint is fixed to 18 m, and the flap 
is limited to the outboard position of 18 m. HLFC is also applied before the folding joint. 
The transition location of the NLF portion of the wing is assumed at 45% [53]. The fuselage 
laminar flow shall be preserved until the wing-body fairing, which is specified based on 
the wing position, which is determined by the center-of-gravity (CG) envelope and the 
empennage characteristics. The active load alleviation effect is estimated by reducing the 
aircraft limit load factor from the required 2.5 [54] to 2.0. For the load alleviation, it is also 
assumed that regulations of CS-25 may allow lower limit loads if sufficient maturity of the 
technology is reached. Advanced structure configurations and material effects are assumed 
to reduce the empty weight by 19% compared to the metal structure. 


4. Top-Level Requirements of the SE7A SR Aircraft 


Top-level requirements have been set to match the reference ATR-72 aircraft. Table 2 
shows the summary of the SE*A SR TLRs. 


Table 2. SE*A SR top-level requirements. 


Requirement Value Units 
Design range for maximum payload 926 km 
Maximum payload 7500 ke 
Cruise Mach number 0.42 

Service ceiling 7620 m 
Take-off field length 1400 m 
Landing distance 1100 m 
Certification CS-25 [54] 


Figure 8 shows a mission profile considered for the design. The mission includes 
the main flight segment, where the airplane cruises at 7300 m, and the required reserve 
segment which consists of the divert segment and the 30 min hold. In addition, 5% battery 
energy contingency was assumed after the reserve mission. 


7000 


6000 


ui 
>) 
oO 
© 


Altitude (m) 
W t+. 
5 8 
oO © 


0 200 400 600 9800 1000 1200 1400 
Range (km) 


Figure 8. SE*A SR mission profile. 
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5. Initial Aircraft Design 
Concept Selection 


In this work, two concepts were considered for the initial aircraft design. The first 
concept was the top-wing conventional configuration with wing-mounted engines. The 
battery was split between the fuselage belly and the engine pylons: one segment was 
located before the landing gear, and the other one was located after the gear. The battery 
distribution was based on the CG envelope to satisfy all possible passenger loading cases. 
Battery modules feature a ’quick-swap mechanism, so a mechanic can switch batteries 
between the flights using a mechanism similar to a high loader. The second configuration 
is the low-wing configuration with aft-mounted engines. Such configuration, unlike the 
first configuration, has a fully clear wing which enables laminar flow control along the 
entire wingspan. Engines are mounted at the aft of the fuselage. The vertical location of 
the engines approached their maximum to minimize the adverse effect from the wing at 
high angles of attack. Batteries also feature the ’quick-swap mechanism. Figure 9 shows 
the described configurations schematically modeled in OpenVSP. 





Figure 9. OpenVSP models of two configurations for the selection. 


Each configuration has its assumptions, features, and limitations. Assumptions are 
based on research goals of the SE*A research groups that focus on the research and devel- 
opment of each technology separately. The first concept has a limited range of laminar flow 
along the wingspan. The region which is subjected to the propeller wash is assumed to 
have a fully turbulent flow. The remaining wing portion has 70% laminar flow due to the 
application of the HLFC. To have a favorable CG envelope, the location of the wing was set 
to 43% of the fuselage length. The second concept has limited laminar flow control before 
the landing gear at the wing pressure side. The remaining pressure side and the suction 
side are 70% laminar. It was also possible to shift the wing further aft to 53% of the fuselage 
length so that extra 10% could be laminarized compared to the conventional configuration. 

Comparison among the concepts was performed using the sizing methodology as 
was described in Section 2. Due to the presence of various novel technologies and, as a 
consequence, increased aircraft development, and operating costs, the concept selection 
was based on the magnitude of Direct Operating Costs (DOC). DOC was calculated using 
the method presented in [55] where the DOC is divided into energy, maintenance, capital, 
crew costs, and fees. 2121 flight cycles per year were assumed for both aircraft. For the 
battery pack, 2000 cycles were assumed. The electricity price was taken from [56] for the 
year 2050 and is equal to 0.118 EUR/kWh that year. The airframe price was assumed similar 
to the ATR-72 [55] and is equal to 1595 EUR/kg. The remaining parameters for all required 
cost components were taken from [55]. A major concern is related to the modification of 
DOC due to novel technologies. It is unknown how much the DOC will increase due to an 
early stage of all technologies presented in this work. However, the complexity related to 
the introduction of the porous skin, the suction system and its structure, which is subject 
to maintenance complexities, NLF front portion, integration of various load alleviation 
technologies, and the structure to support the wing laminarization and load alleviation 
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will increase maintenance costs by an order of magnitude. As an optimistic assumption, a 
scaling factor of 2 was used for the maintenance DOC gain. After the introduction of the 
maintenance DOC gain, the total DOC formulation becomes 


DOC Total — DOCEnerey + DOC Crew + keDOCMa + DOCcap + DOC Fees (11) 


where DOCEnerey are costs of energy, DOCcyew are crew costs, DOC, are maintenance 
costs, DOCcap are capital costs, DOC Fees are costs of fees, and kg is the maintenance cost 
gain factor. Finally, all costs were calculated in 2020 EUR. 

Due to the presence of an all-electric propulsion system, a conventional assumption of 
the best design point at the lowest power-to-weight ratio may not hold, as was previously 
shown by Finger [5] where the best design points were sensitive to the energy type of 
networks and could show points of minimum take-off or fuel weight at wing loadings larger 
than the one with minimum power-to-weight ratio. To ensure that a good initial guess for 
each concept is selected, an aspect ratio and wing loading trade study were performed 
using the initial sizing procedure described in Section 3. The procedure has multiple steps. 
First, the initial aircraft sizing is performed to find the design point corresponding to the 
maximum possible wing loading. Assuming that the selection of an arbitrary design point 
will not change the constraint diagram significantly, a set of sample wing loadings ranging 
from 60% of the maximum wing loading to the maximum possible wing loading were 
selected. Figure 10 shows the sweep of possible design points along the constraint diagram 
boundary. The aircraft is then sized for each selected wing loading using SUAVE, and the 
DOC is estimated. This approach is repeated for a selected set of aspect ratios to determine 
the combination of the wing loading and aspect ratio that may give the lowest desired 
aircraft characteristics important for the decision making. A sample set of weights and 
DOC for a low-wing configuration as a function of the wing loading and the AR is shown 
in Figure 11. Multiple trends can be observed from the Figure 11. First of all, the battery 
weight generally reduces with the AR, while the operating empty weight (OEW) and 
maximum take-off mass (MTOM) trends do not follow such trends. Moreover, DOC and 
weights show discontinuities for some aspect ratios and particular wing loadings. Points 
of the curve discontinuity correspond to the folding penalty, which is introduced starting 
from a specific wing loading. 
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Figure 10. A sample constraint diagram and the range of possible design point or the initial sizing. 


Table 3 shows the summary of the high-wing and the low-wing configurations that 
have the lowest DOC. Based on the given results, the high-wing configuration has larger 
weights due to more battery required to complete the mission with less laminar flow. 
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Consequently, the airframe with less possible drag is more favorable. Consequently, the 
low-wing configuration is more efficient from the performance standpoint. 


Table 3. Comparison between the high-wing and low-wing concepts. 





Parameter High wing Low wing Units 
MTOM (Wo) 38,238 35,369 ke 
OEW (W,) 30,738 27 ,869 ke 
Battery weight (W,) 14,629 12,783 ke 
AR 11.0 11.0 - 
Wing span (b) 40.2 38.7 m 
taper ratio (A) 0.5 0.5 - 
Wing area (Sw) 147.0 136.0 m2 
Power-to-weight 
144 14 
ratio (P/W) : va BENS 
DOC 10,074 9603 EUR/flight 
— AR=9.0 — AR=11.0 — AR=13.0 —— AR=15.0 —— AR=17.0 — AR=18.0 
—— AR=10.0 — AR=12.0 — AR=14.0 — AR=16.0 
MTOM Empty weight 
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Figure 11. Sensitivity of weights and DOC to the aircraft AR and wing loading for the low-wing configuration (pack energy 


density of 700 Wh/kg). 


Due to an unconventional configuration of the low-wing configuration, a concern 
related to the aircraft stability and control may occur. Coe [57] has investigated the sta- 
bility and control characteristics of an aircraft with aft-mounted propeller engines. Re- 
sults demonstrated a reduction of longitudinal stability for this configuration due to the 
propeller—pylon interaction. However, minor destabilization can be compensated by a 
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minor increase of the horizontal tail. Lateral stability is not significantly affected by aft- 
mounted propellers and may improve directional control due to the boundary layer control 
introduced by the propeller slipstream. Riley [538] also provided a detailed description 
of business/commuter turboprop aircraft configuration with aft engines. Although the 
configuration is a twin-pusher, general information about handling qualities during the ap- 
proach can be extracted. Generally, the aircraft can have acceptable handling qualities and 
can be certified, but more care to achieve proper dynamic stability and desired handling 
characteristics may be required. 


6. SE*A SR Design Refinement Using MDO 


To refine the initial design, a set of multi-disciplinary design optimizations (MDO) 
were performed. The main goal of the refinement is to reduce the aircraft DOC further. 
Multiple cases were considered. The first MDO featured all prescribed assumptions. The 
second MDO assumed that no folding penalty is present to estimate the influence of the 
folding penalty of the aircraft. This case is used only as a comparative study to investigate 
the folding mechanism effect on the optimal design. Finally, a special case of battery weight 
minimization was calculated to investigate the reduction of the battery weight and increase 
of DOC with respect to the aircraft optimized for the DOC. This case also directly influences 
the aircraft CO» emission. 

Production of electricity also involves emission CO, and other greenhouse gases, 
since modern electricity generation is not a product of renewable energy sources. A major 
portion of electricity today comes from fossil fuels, consequently, an all-electric aircraft 
cannot be considered as a zero-emission one. Since the goal of the study was to design and 
assess the aircraft from the environmental perspective as well as performance and costs, 
proper metrics of preliminary emission rates must be introduced. Scholz describes a rapid 
and useful procedure of comparing kerosene-based and battery-electric emission [59,60] to 
calculate the CO2 emission including additional factors for each energy source. For aircraft 
operating with kerosine-based fuels, additional chemicals other than CO? contribute to the 
emission level. In addition, the altitude at which the emission happens plays an important 
role. To account for secondary effects, the equivalent CO2 emission level per flight can be 
described by 


MCO2,eq — (Elco2fkm + ElnoxfkmCFmidpoint,NOx 2s CE midpoint, AIC) Rkm (12) 


where fi is the fuel flow per km, EI is the emission index for each species, CF is the 
characterization factor, and Rx, is the range in km. The Elco2 is equal to 3.16 kg CO, 
per kg fuel [60]. The EIyo,x is found using the Boeing Method 2 [61] where the index is 
computed based on existing data from the ICAO Engine Emissions Databank [62] and is 
corrected with respect to the flight altitude. Characterization factors are calculated by 


CE. _ SGTPo35,100 SGT Po3z,100 . SGT PcH4,100 . (13) 
midpoint, NOx SGTPco2.100 O3s SGTPcop 10 O3L SGTPco2,100 CH4 
SGTP trails.100 SGTP.; 100 
CF midpoint, AIC a SCTProp ig Sonia SCTPoop oo “cits (14) 


where SGTP is the sustained global temperature potential, summarized in Table 4 and s is 
the forcing factor which depends on the altitude and is shown in Figure 12 for each species. 
The formulation given in Equation (12) was used for each flight segment of the 
reference ATR-72 and then added to obtain the total mission equivalent emission level. 
Segments with variable altitude such as climb and descent used average values of forcing 
factors s and emission index EIyo,. If the value of the forcing factor was for the altitude 
lower than the available data, the minimum value represented in the Figure 12 was used. 
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Table 4. SGTP values [60]. 

Species SGTP; 100 Units 
CO> 3.58 x 10714 K/kgco, 
Or, 7.79 x 10-2 K/keno, 
O3, —914 x 107-8 K/keno, 
CH, —3.90 x 10-4 K/keno, 
Contrails 1.37 x 10-8 K/km 
Cirrus 4.12 x 10~}8 K/km 
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Figure 12. Forcing factor as a function of altitude [60]. 


Battery emission, on the other hand, depends on the amount of energy used to charge 
them. The emission then depends on the share of fossil fuel to non-fossil fuel sources 
participating in the electricity production and the primary efficiency factor—the factor 
which describes the ratio between the energy required to produce electricity and the final 
output energy amount. Including all factors described above, the equivalent emission per 
flight becomes 
kpEFMpatChat 


(15) 
NchargeC Jet—A 


Mco2,eq = Elco2Xx fF 
where x ¢¢ is the fossil fuel share ratio equal to 0.39 according to the forecast of 2050 [63], 
kper is the primary efficiency factor assumed equal to 2 [59], 1charge is the battery charging 
efficiency equal to 0.9 [59], and Cyg; and Cje+__4 are battery and Jet-A fuel energy densities, 
respectively. From Equation (15), the CO emission is directly proportional to the battery 
mass, so the optimization for the minimum battery weight corresponds to the minimum 
emission case. 
All problems have similar design variables and constraints, as shown in Table 5. 
Due to the application of new airframe technologies, requirements for the sea-level climb, 
take-off, and landing are more significant compared to cruise and OEI requirements, as 
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shown in Figure 10. Such behavior occurs due to a substantial reduction of the cruise 
drag, which strongly diminishes the required cruise power. Consequently, optimization 
problem constraints can be limited to these three main performance metrics. Additional 
geometric constraints of 51.8 m wingspan (similar to the wingspan of the NASA N+3 
strut-braced wing concept [64]) and the wingtip of no larger than 1.4 m were introduced 
to limit further increase in the wing slenderness. Table 5 describes the formulation of the 
optimization problem. There, C, and C, are the wing root and tip chords, respectively, t/c 
is the wing thickness, and 7,,,, is the maximum throttle during the flight. To perform the 
MDO, SUAVE was coupled with MATLAB, and the Genetic algorithm was used to find 
optimal solutions. 


Table 5. SE*A SR aircraft optimization problem definition. 


Lower Upper Units 
minimize 1. DOC 
2. Wp 
wrt AR 8.00 16.00 
A 0.25 0.75 
o 3.00 7.00 m 
ae ee 12.00 18.00 %o 
t/C\riy 12.00 18.00 % 
P/W 0.05 0.30 
subject to Take-off field length (TOFL) 1400.0 m 
Landing field length (LFL) 1100.0 m 
Nmax 1.0 
b 50.0 m 
C, 1.4 m 


To account for the constraints, the penalty function similar to the one defined in [65] 
was used. The penalty function is defined by 


y \7 

ip uly y)(“) (16) 
Yc 

where y is the design variable, y, is the design variable constraint, 1 is the unit step function 

equal to zero for y < y-, and y = 3. With the introduction of the penalty function, the 

objective function becomes 


bat eae Bhi: (17) 


where N is the total number of design variables. Sixty species per generation were set 
to have sufficient population size without major accuracy losses. Take-off and landing 
constraints were calculated using physics-based time-dependent performance formulations 
presented in [41]. The climb power-to-weight ratio constraint was calculated using a similar 
formulation used for the constraint analysis. 

Figure 13 shows the aircraft planforms for each optimization problem and Table 6 
summarizes important properties of each optimized configuration. Limiting constraints 
for each optimization were similar to the constraint diagram: the solution was trying to 
minimize the power-to-weight ratio which reached the constraint limit of either the climb 
or the take-off. In addition, the top-of-climb throttle constraint was an additional constraint 
which dominated the others. 
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Figure 13. Wing planforms for selected optimization problems. 
Table 6. Comparison between optimized concepts. 


DOC- DOC-Optimal _ Battery- 


Initial 


Parameter Design Optimal Design (No Optimal ATR-72 Units 
Design Folding) Design 
Wo 35,369 34,441 33,649 34,783 23,000 ke 
We 27,869 26,941 26,149 27,283 13,311 ke 
W, 12,783 12,262 11,527 10,379 - ke 
W fuel,mission ~ 7 ~ ~ 1913 kg 
AR 11.0 10.86 12.23 16.0 12.0 : 
b 38.7 39.6 41.3 51.7 27.05 m 
A 0.5 0.31 0.3 0.285 0.55 - 
c, 4.68 5.83 5.23 5.00 2.74 m 
ui 136.0 145.4 139.6 165.5 61.0 m? 
ated Dee 0.18 0.18 0.18 0.17 0.18 : 
t/C\ tin 0.13 0.13 0.13 0.12 0.13 - 
P/W 0.143 0.125 0.13 0.11 0.16 kW/ke 
DOC 9603 9479 9397 9838 6740 a. 
ight 
MCO, 2045 1961 1844 1660 11,838 ke 
TOFL 1338 1381 1368 1360 1367 m 
LFL 1025 972 985 934 1207 m 


As shown in Table 6, the optimized aircraft with all imposed assumptions reduced all 
weights, slightly reduced the aspect ratio, increased the wingspan, reduced the taper ratio, 
and reduced the aircraft power-to-weight ratio. Minor differences in the field performance 
occurred. The DOC for the optimized aircraft reduced by 124 EUR/flight. In addition, the 
difference in emission between the initial and DOC-optimized designs is equal to 84 kg 
per flight which is not substantial. If the folding penalty is removed, then the wing has 
more freedom to increase its span and aspect ratio due to the absence of the folding weight 
penalty. The optimized aspect ratio increased by 1.23, the wingspan increased by 2.6 m 
compared to the initial design. However, the power-to-weight ratio increased due to a 
reduced wing planform area. All weights were further reduced compared to the aircraft 
configuration with the folding mechanism. The DOC was reduced by 206 EUR/flight 
compared to the initial configuration, which indicates a relatively strong effect of the 
folding mechanism on the aircraft. The emission level reduced by 201 kg compared to 
the initial design which is more substantial. Comparing the DOC optimal solutions with 
and without the folding mechanism, it becomes evident that wing folding affects the final 
design and cannot be neglected. If the battery minimization strategy is considered, then 
the wing planform approaches design limits. Due to a stronger reduction of the battery 
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weight compared to the increase in the wing weight, minor reductions of the MTOM and 
empty weights are observed. Finally, DOC increased by 235 EUR/flight compared to 
the initial design or by 359 EUR/flight compared to the optimized configuration. The 
battery weight reduced by 2404 kg and 1883 kg compared to the initial and DOC-optimized 
aircraft, respectively. The emission level, however, reduced by 385 kg compared to the 
initial configuration. In addition, the wing span of the aircraft increases to 51.7 m and the 
folding portion from each side of the wing becomes 7.85 m. Such a large folding portion 
may affect the operational convenience at the airport and will make the aircraft taller when 
the wings are folded. That may create additional difficulties related to the aircraft storage. 

If DOC and emission levels are compared against the reference aircraft, the increase 
in DOC demonstrates a substantial increase compared to the reference ATR-72. The most 
cost-effective aircraft version has 41% increase in DOC which is a substantial cost growth 
from the market standpoint. The DOC for the battery-optimal configuration increases by 
46%. On the other hand, if emission levels are compared, the most DOC-efficient version 
has a reduction in equivalent CO2 emission of 83% compared to the ATR-72, and the 
battery-optimal configuration has 86% emission reduction. From the comparison, the 
emission level of the all-electric aircraft reduces substantially, but such reduction of the 
emission comes at a significant cost increase. Moreover, DOC shows mote sensitivity with 
respect to the configuration compared to the emission level. To minimize costs as much as 
possible, the DOC-optimal configuration was selected for future analyses. 

Observing the wingspan deviations from the Part C airport requirements of 36 m, 
aircraft configurations minimized for DOC did not show major increases in wingspans. 
Consequently, it is possible to design an aircraft that will have no folding mechanism and 
will not have a significant DOC penalty so that the wing design complexity may be reduced 
substantially. An additional MDO study with a fixed wingspan constraint was performed 
to investigate the sensitivity of DOC to the configuration with a fixed span of 36 m. Table 7 
demonstrates a comparison between the initial design, the optimal design without the span 
constraint, the optimal aircraft with a constrained span, and a reference ATR-72 to compare 
optimal models to their reference. 


Table 7. Comparison between optimized concepts with and without the wing span constraint. 


Initial DOC- DOC-Optimal 

Parameter Design Optimal Design ATR-72 Units 
Design (Fixed Span) 

Wo 35,369 34,441 35,745 23,000 keg 
We 27,869 26,941 28,244 13,311 kg 
W, 12,783 12,262 13,861 : ke 
W fuel,mission 7 7 ~ 1913 kg 
AR 11.0 10.86 10.0 12.0 - 
b 38.7 39.6 36.0 27.05 m 
A 0.5 0.31 0.29 0.55 - 
oe 4.68 5.83 0.07 2.74 m 
ae 136.0 145.4 129.3 61.0 m? 
cl ee 0.18 0.18 0.18 0.18 - 
t/c fe 0.13 0.13 0.12 0.13 - 
P/W 0.143 0.125 0.14 0.16 kW/keg 
DOC 9603 9479 9573 6740 EUR/flight 
MCO» 2045 1961 2217 11,838 ke 
TOFL 1338 1381 1392 1367 m 
LFL 1025 972 1027 1207 m 


The aircraft with a fixed span has a higher weight due to a restricted planform and, as 
a consequence, more battery is required for the mission. The battery weight increased by 
1600 kg compared to the optimal design without the span constraint. The power-to-weight 
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ratio also increased to satisfy the take-off requirement. However, the wing area is reduced, 
so it is less expensive to increase the motor power than increase the wing area. The DOC 
increased by 94 EUR/flight, which is not a large increase. At the same time, the wing has 
no folding mechanism and is easier to design and maintain. However, the emission level 
increases by a relatively small margin compared to the initial DOC-optimal configuration. 
At this moment, the configuration without folding is selected as the final decision. 

The geometric summary of the selected aircraft is shown in Table 8. Single-slotted 
Fowler flaps were used for this aircraft. The empennage was initially sized within SUAVE 
using a constant tail volume fraction similar to existing aircraft in the class using data 
from Raymer [40] and then refined using AVL to ensure sufficient stability and trim during 
critical flight cases for the complete CG envelope obtained separately. The vertical tail was 
sized based both on the one-engine inoperative (OEI) case and the lateral stability condition 
of C Ng > 0.01573 [66]. The propeller was positioned vertically as far as possible to reduce 
the possibility of the wing wake impinging on the propeller blades. The angle between the 
wing trailing edge and the propeller’s lowest blade tip is equal to 14 deg, which should be 
sufficient for the majority of operational cases. The propeller features six blades and has a 
diameter of 3.56 m. 


Table 8. Geometric properties of the SE*A SR aircraft. 


Parameter Wing Horizontal Tail Vertical Tail Units 
AR 10.0 5.07 1.25 - 
b 36.0 10.0 3.5 m 
A 0.29 0:55 0.8 - 
Ac/4 0.0 15.0 26.0 deg 
r 4.0 0.0 0.0 deg 
C, 5.58 2.5 3.13 m 
Sree 129.3 17.8 9.80 m? 
aes ae 0.18 0.10 0.10 - 
t/C\tiy 0.12 0.10 0.10 - 
Flap span ratio 0.6 - - - 
Flap chord ratio 0.2 7 - = 
Aileron span ratio 0.25 = : : 
Aileron chord ratio 0.20 ; - - 
Elevator span ratio - 1.0 - - 
Elevator chord ratio - 0.25 - - 
Rudder span ratio - 1.0 0.8 - 
Rudder chord ratio - 0.25 0.25 - 
Vut - 0.59 = = 
Wr - : 0.024 - 


The CG envelope was created using available information about the stick-fixed and 
stick-free neutral points, payload clouds, and multiple baggage arrangement cases. The 
CG-range is equal to 21.5% mean aerodynamic chord (MAC) starting at 12% and ending at 
33.5% MAC. 

Figure 14 shows the SE7A SR payload-range diagram compared to the reference ATR- 
72 [67]. In addition, the obtained payload-range diagram was validated using a Breguet 
range equation for an all-electric aircraft defined by 


1 L Mp 
K= g Total Coat 5 a 





a (18) 


where mmo is the maximum take-off mass, g is the gravitational acceleration, and L/D is the 
cruise lift-to-drag ratio. 
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Figure 14. Payload-range diagrams for the SE7A SR and ATR-72. 


For the range analysis, multiple options were considered. The first option is the case 
when the battery mass remains constant, independent of changes in payload weight. This 
way, the battery volume is fixed and equal to the one used for the maximum passenger case. 
The second option is to have sufficient battery volume to replace missing passengers with 
an extra battery. Such an option is more problematic to achieve due to the increased battery 
size and additional growth of the already large belly fairing but is useful to analyze for the 
range sensitivity analysis. The comparison shows that although the SE*A SR satisfies the 
harmonic range requirement, which matches the ATR-72, the ferry range is significantly 
limited due to the battery energy density and weight limitations. The absence of the 
payload extends the aircraft range to 1415 km, which is substantially smaller than the ATR- 
72 with its ferry range of 3380 km. Introduction of more battery which replaces the payload 
increases the ferry range until 1870 km, which is 455 km longer than the fixed battery 
weight case. However, this ferry range still does not introduce major range improvements 
compared to the ATR-72. Comparison between the SUAVE payload-range and Breguet 
payload-range diagrams show a difference of both harmonic and ferry ranges of 100 km. 
Such tendency happens due to the higher fidelity of SUAVE which includes all mission 
segments and calculates incremental energy consumption and aerodynamic characteristics, 
unlike the Breguet equation which considers a constant lift-to-drag ratio. On the other 
hand, minor deviations are expected and are satisfactory for the selected analysis. 

To summarize current design outcomes, an all-electric SE7A SR aircraft is capable of 
achieving the harmonic range similar to the reference ATR-72 and satisfies all prescribed 
TLRs. Moreover, its emission level is reduced by 81% compared to the reference aircraft. 
On the other hand, the ferry range of the all-electric aircraft is limited to no more than 
1870 km compared to 3380 km for the ATR-72 and has a 42% increase in DOC with respect 
to the reference which will create substantial market application problems for this type 
of airplane. 

Figure 15 shows a rendered image of the optimized SE*A SR aircraft. 
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Figure 15. SEA SR aircraft geometry. 


7. Sensitivity Analysis of the Aircraft Characteristics to Technology Advancements 


After performing a conceptual design of a regional all-electric aircraft and analyzing 
obtained results, several questions still remain: What is the influence of novel technologies 
on the aircraft configuration? How much do the deviations of novel technology assump- 
tions affect the airplane configuration? What technologies affect the aircraft geometry 
the most? 

A sensitivity analysis of the SE*A SR aircraft to the battery energy density with the 
absence of all novel airframe technologies was performed to determine what potential 
energy density could make all-electric regional aviation possible and how the aircraft 
characteristics are affected. The aircraft without novel airframe technologies assumed 5% 
laminar flow on the wing and empennage, 5% laminar flow for the fuselage, a limit load 
factor of 2.5, and no weight reduction due to composite materials. Technologies related 
to the propulsion system such as HTS motors remained since their absence immediately 
leads to the aircraft’s infeasibility due to the Megawatt-level power required per motor. 
Furthermore, the performance characteristics of the aircraft without novel technologies 
were compared against the case if all technologies met the desired requirements. The 
sizing procedure for all configurations was performed for a similar mission profile and 
used similar methods to those described in Section 3 to ensure the satisfaction of all TLRs. 
Similar constraints based on TLRs were also implemented. As for the initial sizing, wing 
loading and aspect ratio sensitivity analyses were performed to determine a sweep of all 
possible aircraft that satisfy TLRs. For all possible configurations, the ones with minimum 
DOC were selected for each selected battery energy density. Finally, all characteristics 
were compared against the reference ATR-72 aircraft. Figure 16 shows the effect of battery 
energy density on the aircraft gross weight. 

Results demonstrate that the absence of novel technologies significantly limits the 
feasibility of all-electric regional aviation unless substantially higher battery cell energy 
densities are reached. For instance, if the aircraft with all available airframe technologies 
can reach the MTOM of 40,000 kg at the pack energy density of 600 Wh/kg, the aircraft 
without airframe technologies can reach similar weights at 900 Wh/kg. In addition, 
weights approach an asymptotic value at high energy densities, so the effect of airframe 
technologies becomes more significant for high energy densities. For instance, if the pack 
energy density reaches 1000 Wh/kg, then the aircraft MTOW with all technologies becomes 
slightly less than 30,000 kg while similar weight is achieved by the aircraft without airframe 
technologies only at 1500 Wh/kg. 
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Influence of battery energy density on aircraft characteristics 
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Figure 16. Sensitivity to the battery energy density based on minimum DOC concepts. 


To compare the DOC values, the aircraft without technologies did not include an 
additional maintenance cost penalty introduced for the aircraft with all airframe technolo- 
gies. Moreover, a sensitivity analysis of the maintenance cost gain factor was performed 
to investigate its effect on the overall DOC and predict potential future DOC behavior. 
If DOC is compared, two trends can be observed. First, both aircraft have substantially 
higher DOC compared to the reference ATR-72, and the difference remains large even for 
high energy densities. On the other hand, the application of technologies affects the DOC. 
Comparing aircraft with and without technologies, there exists a point after which the 
DOC increases due to the airframe technologies and their maintenance complications and 
is equal to 980 Wh/kg. After this energy density, an aircraft without airframe technologies 
becomes cheaper to operate. Such a trend is a consequence of the maintenance cost gain 
and redistribution of dominant cost drivers, as shown in Table 9. For low battery energy 
densities, the most dominant DOC component is the capital cost which depends on the 
airframe and propulsion system weights. Due to significantly larger weights, these costs 
contribute to total DOC such that benefits of maintenance costs disappear. Moreover, more 
battery energy is required for the aircraft without technologies and slightly higher fees also 
contribute to a larger total DOC. On the other hand, with the increase in energy density, air- 
craft weights reduce substantially, so capital costs reduce as well. The overall distribution 
of costs becomes more evenly distributed for the aircraft without technologies, so lower 
maintenance cost effects become more significant. On the other hand, maintenance costs 
for the aircraft with technologies start dominating. That relationship leads to higher overall 
DOC for the aircraft with technologies at high battery energy densities. The present aircraft 
was designed having 700 Wh/kg of pack energy density which is still more beneficial than 
the aircraft without technologies. However, the maintenance penalty factor of 2 may be 
too optimistic. If the maintenance costs increase more, then the equilibrium point between 
the two DOCs will shift towards lower energy densities, as shown in Figure 16. The only 
possibility of approaching the reference aircraft DOC exists if maintenance cost gain is 
equal to 1. The maintenance cost of novel airframe technologies shall be considered later to 
ensure the financial applicability of novel technologies further. 
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Table 9. Comparison of DOC components depending on the battery energy density for the SE7A SR 
aircraft in EUR/flight. 


DOC 700 Wh/kg 700 Wh/kg 1300 Wh/kg 1300 Wh/kg 
Component with without with . without 
Technologies Technologies Technologies Technologies 
Energy 808 1438 700 1037 
Crew 1369 1369 1369 1369 
Maintenance 3199 1891 2937 1596 
Capital 2233 3554 1928 2514 
Fees 1994 2350 1799 1928 
Total 9603 10,604 8734 8446 


The effect of technology deviations on the aircraft weights and, as a consequence, 
applicability of all-electric aircraft is shown in Figures 17-19. In addition, points for 
extreme technology cases are summarized in Tables 10-13. A similar sizing approach, 
as discussed in Section 3, was used for the sizing with technology deviations. The wing 
loading, however, was fixed to the one that occurred after the MDO of the aircraft. In 
addition, to simplify the sensitivity and account for significant complexities to design an 
appropriate suction system for the fuselage, the fuselage was assumed fully turbulent for 
all technology sensitivities. 
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Figure 17. Sensitivity of aircraft characteristics to changes in technology levels at Structural weight 
reduction of 0.0. 
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Structural weight reduction = 0.1 
AR = 10, Cruise altitude = 7300m; Mach = 0.42 
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Figure 18. Sensitivity of aircraft characteristics to changes in technology levels at Structural weight 
reduction of 0.1. 
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Figure 19. Sensitivity of aircraft characteristics to changes in technology levels at Structural weight 
reduction of 0.19. 
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Table 10. MTOM values for the airframe weight reduction of 0%. 


70% Laminar Flow 45% Laminar Flow Units 
limit load 1.5 39,665 42,500 ke 
limit load 2.5 42,100 46,800 ke 


Table 11. MTOM values for the airframe weight reduction of 19%. 


70% Laminar Flow 45% Laminar Flow Units 
limit load 1.5 37,300 39,800 ke 
limit load 2.5 38,900 41,800 ke 


Table 12. MTOM values for the limit load factor of 1.5. 


0% Airframe 19% Airframe Units 
Weight Reduction Weight Reduction 
70% laminar flow 39,665 37,300 ke 
45% laminar flow 42,500 39,800 ke 
Table 13. MTOM values for the limit load factor of 2.5. 
0% Airframe 19% Airframe Units 
Weight Reduction Weight Reduction 
70% laminar flow 42,100 42,500 ke 
45% laminar flow 46,800 41,800 ke 


Based on obtained results, if the weight sensitivity on the technology is compared for 
two other technologies where one technology has a constant value and the other technology 
is fixed at either their minimum or maximum, the laminar flow control shows the largest 
influence on the MTOM. Structural weight reduction plays the second most important 
role, and the load alleviation affects the aircraft MTOM the least compared to the other 
two technologies. For instance, for the airframe reduction of 0%, if the maximum possible 
boundary layer transition is possible, then the change in MTOW due to the load alleviation 
is equal to 2435 kg. However, if for the same airframe reduction factor, the load alleviation 
has its minimum value, then the change in MTOM is equal to 2835 kg due to the HLFC 
technology. If technologies are combined, then the effect of the MTOM is magnified. 


8. Conclusions 


The present manuscript investigated the influence of novel technologies on the fea- 
sibility of all-electric short-range aircraft with passenger capacity and performance char- 
acteristics similar to ATR-72. Three novel airframe technologies were considered: hybrid 
laminar flow control, active load alleviation, and advanced materials and structure con- 
cepts. In addition, superconducting electric motors and high energy density batteries were 
considered for the aircraft energy system. 

Multiple tools were used and improved to assess the feasibility question. SUAVE 
was used to perform the initial sizing and sensitivity analyses and was also coupled with 
MATLAB to perform MDO studies. EMWET and FLOPS weight estimation methods and a 
weight penalty function for folding wings were integrated into SUAVE. The cubic spline 
method for a propeller was used to enable the aircraft sizing. AVL was used to assess the 
aircraft stability and control. 

Two aircraft configurations were considered for the design: the conventional con- 
figuration with a high wing and wing-mounted propellers and the low wing with pro- 
pellers mounted at the aft of the fuselage. The low-wing configuration was selected 
because it had lower DOC, lower weights, and is possible to have certifiable stability and 
control characteristics. 
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To refine the selected concept and determine the aircraft characteristics’ sensitivity to 
the folding penalty and the objective function, a set of MDO studies were performed. The 
configuration without folding wings was selected due to a relatively small increase in DOC 
due to the fixed wingspan and reduced complexity of the wing design and maintenance 
due to the absence of the wing folding mechanism. 

The importance of airframe technologies was assessed by comparing aircraft con- 
figurations with and without them. The presence of airframe technologies may reduce 
the required battery pack energy density by 300 Wh/kg compared to their absence. The 
difference in required energy density also increases with the increase of the battery energy 
density. However, the application of technologies becomes more expensive than their 
absence if the battery pack energy density exceeds 980 Wh/kg for a given maintenance 
cost gain factor of 2. 

Comparison of the all-electric aircraft to the reference ATR-72 showed substantial 
reduction of CO emission by 81% of the all-electric aircraft. On the other hand, DOC 
increase of 42% presents a significant limit to market success of such aircraft. Moreover, 
neither increase in battery energy density, nor reduction of maintenance costs change 
aircraft DOC significantly. The only possibility to have a comparable DOC level may appear 
if maintenance cost gain is absent and battery pack energy density exceeds 1000 Wh/kg. 
Moreover, the all-electric aircraft consumes a lot of energy to complete the mission which 
creates major uncertainties related to its energy-efficiency from the operational standpoint. 

Finally, a sensitivity analysis of technology availability on the aircraft characteristics 
was performed. Simulations were performed for each technology combination from the 
most optimistic to its absence. Results demonstrated high sensitivity of aircraft weights to 
technology levels where the HLFC has the most influence, followed by the aircraft weight 
reduction and finishing with the load alleviation. The combined effect of technologies 
magnifies the differences in weight even further. 

Future research steps will include a detailed conceptual design of the SE*A SR aircraft. 
First of all, the design with available models for the laminar flow required power and 
additional weight penalties due to the implementation of all technologies and additional 
hardware will be performed. A more sophisticated assessment of the aircraft stability 
and control must be performed to ensure adequate flying qualities. Finally, higher fi- 
delity analysis for the energy network and modeling of the motors and propellers will 
be considered. 
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Abstract: This paper proposes a feasibility study concerning a large turboprop aircraft to be used 
as a lower environmental impact solution to current regional jets operated on short/medium hauls. 
An overview of this market scenario highlights that this segment is evenly shared between regional 
turboprop and jet aircraft. Although regional jets ensure a large operative flexibility, they are usually 
not optimized for short missions with a negative effect on block fuel and environmental impact. 
Conversely, turboprops represent a greener solution but with reduced passenger capacity and speed. 
Those aspects highlight a slot for a new turboprop platform coupling higher seat capacity, cruise 
speed and design range with a reduced fuel consumption. This platform should operate on those 
ranges where neither jet aircraft nor existing turboprops are optimized. This work compares three 
different solutions: a high-wing layout with under-wing engines installation and both two- and 
three-lifting-surface configurations with low-wing and tail tips-mounted engines. For each concept, 
a multi-disciplinary optimization was performed targeting the minimum block fuel on a 1600 NM 
mission. Optimum solutions were compared with both a regional jet such as the Airbus A220-300 
operated on 1600 NM and with a jet aircraft specifically designed for this range. 
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ee From a preliminary overview, nowadays, the regional aircraft market segment is 


mainly influenced by a combination of three factors. A first key parameter to be considered 
is the environmental impact reduction, as can be seen by ever-more ambitious targets 
envisaged by international associations such as the Air Transport Action Group (ATAG) 
and the International Air Transport Association (IATA) or the Clean Sky 2 consortium. In 
2008, the ATAG board developed a set of environmental goals for the short, medium, and 
long term which were supported and reiterated by the [ATA Board and the association’s 
Annual General Meeting [1]. Those goals can be summarized in the followings: an average 
improvement in fuel efficiency of 1.5% per year from 2009 to 2020; a cap on net aviation 
CO emissions from 2020 (carbon-neutral growth); and a reduction in net aviation CO? 
emissions of 50% by 2050, relative to 2005 levels. However, as pointed out in [2], most 
of the targets forecasted by IATA and ATAG for 2020 have not been achieved, and the 
expected carbon-neutral growth appears to be quite far. This highlights the need to stress 
even more the reduction of the civil aviation environmental impact as the main driving 
5 parameter in the design of new and greener aircraft models. 

In terms of innovation aiming at reducing civil aviation environmental impact, the 
Clean Sky 2 Program of Horizon 2020 has accelerated by now the development and 
introduction of new technologies designed for entry into service in the 2025-2035 timeframe. 
conditions of the Creative Commons  Xccording to the Clean Sky 2 development plan, by 2050, 75% of the world’s fleet now in 
Attribution (CC BY) license (https://. Service (or on order) will be replaced by aircraft that can deploy Clean Sky 2 technologies [3]. 
era tivecommons-o8e/licenses/by/ High level objectives defined by Clean Sky 2 expects a reduction in CO2, NOx, and 
4.0/). environmental noise from —20% to —30% in the 2014-2024 timeframe. Furthermore, 
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reductions of —75% in CO», —90% in NOx, and —60% in environmental noise have been 
forecast by 2050. 

In addition to the environmental impact topic, a second main factor currently influ- 
encing almost all market segments is represented by the recent COVID-19 outbreak, which 
has caused severe damage to the worldwide aviation industry. 

Market forecasts made by the most important aircraft manufacturers at the beginning 
of 2020 have been completely changed by this unexpected event resulting in an overall 
reduction of air passengers (both international and domestic) of 50% in 2020 compared to 
2019 [4]. In addition, the latest estimates made by the International Civil Aviation Orga- 
nization (ICAO) [4] have highlighted a passenger demand drop from the 2020 originally 
planned forecasts by 2690 million (—60%), which has been related to a 370 billion US$ loss 
of gross passenger operating revenues. A partial recovery has been forecast by ICAO [4] for 
the first half of 2021 reducing the passenger revenue loss faced in 2020, equal to —369,768 
million US$, to a value between —178,373 and —148,719 million US$. 

However, when a full recovery from COVID-19 damages will be reached, most ma- 
jor airlines will still have to face the third main factor influencing the regional market 
scenario, which is represented by the need to replace several hundred heritage aircraft, 
especially in the segment from 20 to 150 seats, which are now close to the end of their 
useful commercial life. 

The combination of all these factors will deeply influence the current aviation industry 
shaping the new generation of aircraft. On the one side, ever more mutable and demanding 
requirements will be asked by the market, resulting in a very challenging design process to 
come up with new aircraft solutions. On the other side, in the near future, new ways and 
possibilities will be accessible to the aerospace research sector with innovative technologies 
and disruptive new aircraft concepts probably making their appearance on the upcoming 
market scenario. 

To have an idea of which could be possible future market needs in the post-pandemic 
scenario, especially in the regional market segment, some interesting information can be 
obtained from recent market analyses made by the major aircraft manufacturers concerning 
the small single-aisle aircraft segment up to 150 seats. 

Dealing with possible future aircraft requirements, Airbus Global Market Forecast 
(Airbus S.A.S., Blagnac, France) [5] provides a useful statistical trend which highlights a 
continuous increment in average seats capacity and mission ranges for the single-aisle jet 
fleet. Figure 1 shows that up to 2017 aircraft seats average has increased from 140 to 169, 
while average mission ranges have experienced an increment from 422 NM to 586 NM. 
Figure 1 also highlights that there is a significant variation around the mean which can be 
attributed to different aircraft sizes. In addition to showing the wide spectrum of operations 
for which the airlines use these aircraft, Figure 1 demonstrates why the range capability of 
aircraft products is an important consideration for airlines and manufacturers, a capability 
which also equates to flexibility. 


Aircraft seating flown 
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Figure 1. Distribution of single-aisle aircraft seats and ranges over the years. Reprinted from ref. [5]. 


As shown in Figure 2, up to 2016, both turboprop and regional jet aircraft were 
sharing the regional market evenly with a slightly increasing advantage of regional jets 
due to a higher seat capacity and longer mission distances than conventional regional 
turboprops. However, despite the regional jets success, turboprop engines are 10-30% more 
efficient than jet engines in terms of Specific Fuel Consumption (SFC), leading to a potential 
reduction of the amount of fuel used per mission as well as pollutant emissions. In fact, 
according to ATR forecasts [6], assuming all short haul flights operated by regional jets 
today are replaced by modern turboprops, 11% of overall regional aviation CO2 emissions 
could be saved. 

To design an advanced turboprop aircraft that could be competitive with respect to 
regional jets on short and medium haul, several key aspects must be assessed if compared 
to conventional turboprop designs. There is the need to increase the seats capacity, the 
aircraft must be designed for longer missions, and cabin noise and passenger comfort must 
be improved, as well as both the aerodynamic efficiency and the maximum cruise speed. 

The authors have already faced the design of an advanced turboprop coupling higher 
seat capacity, improved cruise speed, and longer design range with a reduced fuel consump- 
tion within the framework of the European funded project named Innovative turbopROp 
configuratioN (IRON) [7-12]. The project deals with the design of an innovative large 
turboprop aircraft with rear engine installation. 


100% 
90% 
80% 
70% 
60% 
50% 


40% 





o nw © a oO 
fo.) a ® ro>) ° 
ro) ro) o ro) o 
- - - - N 


Figure 2. In-service fleet share of turboprops and regional jets. Reprinted from ref. [13]. 
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Such a configuration shows several design issues dealing with the rear engine in- 
stallation. The most comprehensive reference dealing with all key design aspects for a 
high-capacity propeller driven aircraft, including the rear engines installation, is the study 
performed by Goldsmith in [14], where a feasibility studies related to a possible conversion 
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of the McDonnell Douglas (Saint Louis, MO, USA) DC-9 from turbofan to turbo/prop-fan 
power plant is discussed. 

Goldsmith’s work shows useful evaluations dealing with the aircraft configuration 
feasibility, aerodynamics, weights, balance, and performance identifying three promis- 
ing configurations with respect to the reference aircraft platform assumed as the DC-9 
Super 80. The configurations proposed by Goldsmith exploited the results of the research 
programs carried out by NASA to develop technologies needed to implement potential 
fuel savings by adopting advanced turboprop engine architectures. These works high- 
lighted 15-30% savings in aircraft block fuel compared to the state-of-the-art turbofan at 
the time by adopting advanced turboprop engines characterized by multi-bladed and high 
efficiently propellers that could be used at higher Mach numbers [15,16]. 

The selection of the most suitable configurations was made after a first qualitative 
comparison between many different possible layouts, followed by a multi-disciplinary 
quantitative analysis of three possible solutions. The first layout was an upper-wing- 
mounted propfan, indicated as Configuration 1 in Figure 3; the second one was a T-tail 
configuration with aft fuselage-mounted propfan, named Configuration 2 in Figure 3; 
and the last aircraft concept provided for a horizontal tail-mounted propfan, indicated 
as Configuration 3 in Figure 3. All modifications to the baseline platform were made by 
keeping constant both wing area and Maximum Take-Off Weight (MTOW) of the aircraft. 


Configuration 1 Configuration 2 
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Figure 3. Configuration selected by the NASA feasibility study. Reprinted from ref. [14]. 


As shown in [14], the first and third layouts were identified as promising solutions 
providing for a block fuel reduction of about 14% and a mean Direct Operating Cost (DOC) 
reduction ranging from 6% to 8% with respect to the reference aircraft. The rear-mounted 
engines installation provided a much larger center of gravity excursion, a higher hori- 
zontal tail weight, and higher values of the equivalent parasite area. However, while 
the T-tail configuration led to lower performance and DOC improvements than the wing- 
mounted engines installation, the horizontal tail-mounted engines layout mitigated the 
above-mentioned negative effects, allowing to reach similar, or even better, performance 
and DOC improvements, if compared with the conventional wing-mounted engine config- 
uration. Moreover, in the case of turboprop or open rotors, rear engine installation reduces 
the cabin noise providing for a better passenger’s comfort. 

The high propulsive efficiency of propellers nowadays still makes them quite an attrac- 
tive solution for regional aircraft flying up to high subsonic speeds. Different concepts for a 
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large turboprop aircraft can be found in the literature. A novel concept for highly efficient 
and ecological propeller driven aircraft is proposed in [17]. This aircraft, characterized by a 
high wing, T-tail, and four turboprop engines with large propeller diameters, has shown 
a potential Direct Operating Cost (DOC) reduction of about 11% while reducing trip fuel 
mass and therefore CO, emissions by about 14.9% compared to the reference aircraft Airbus 
A320 (Airbus S.A.S., Blagnac, France). A modified architecture of the concept proposed 
in [17] is shown in [18]. This solution, featured by a strut-braced wing with natural laminar 
flow, could potentially reduce fuel consumption and CO 2 emissions by about 36% with 
respect to a reference jet aircraft (A320). In addition, DOC could be potentially reduced by 
about 17% on a mission range of 755 NM. The improvements mainly come from the lower 
fuel consumption of turboprop engines compared to turbofan engines. Drawbacks are the 
21% increase in block time due to the lower cruise speed compared to a jet aircraft and the 
additional mass due to the higher landing gear lengths caused by high propeller diameters. 

Recent studies concerning innovative turboprop aircraft designs are shown in [19], 
where two configurations featuring tail-mounted propeller and fuselage-mounted ducted 
propellers have been assessed in terms of their key performance indicators assuming they 
are designed for a harmonic mission of 1530 km and 68 passengers and compared with 
a conventional aircraft with wing-mounted propellers similar to ATR-72-600 in terms of 
seometry and weight. 

Although previous research works have already dealt with the comparison between 
existing turbofan aircraft and future turboprop models to highlight possible improvements 
in environmental impact and DOC, none of them starts from the design of a modern 
turboprop aircraft based on a set of Top-Level Aircraft Requirements (TLARs) related to 
a short/medium haul (e.g., 1600 NM) coupled with a seat capacity in line with market 
analysis trends (e.g., 130-150 seats). Conversely, they start from an existing turbofan 
platform trying to make a conversion to a different propulsive system. In addition, most of 
the studies are focused or on aircraft designed for longer missions with a slightly higher 
passenger number, or they focus the attention on smaller aircraft with similar design ranges 
but with a lower passenger capacity. 

This work aims to investigate the combination of the selected mission range and seats 
capacity to design and optimize a possible modern turboprop platform suitable for an 
entry into service expected by 2035 and targeting the minimum environmental impact. 
Since there are no existing regional turbofans specifically designed for such a design range 
and seats capacity, this work does not limit its application in comparing the provided 
turboprop solutions with the current state of the art regional jet but also faces the design 
and optimization of a completely brand-new jet platform to perform fair comparison 
between a turboprop and jet driven aircraft both optimized for the same set of TLARs. 

Such a comparison could represent useful information for aircraft designers who want 
to develop new aircraft suitable for regional applications by providing a good starting 
point in the decision process related to the preliminary design phase. 


2. Aircraft Requirements and Critical Issues 


The feasibility study starts from the definition of a specific set of TLARs, suitable for 
modern regional transport applications and in line with the main aircraft manufacturer 
market forecast, which is better illustrated in the Introduction. The TLARs driving the 
design of a large capacity turboprop are summarized in Table 1. As can be appreciated 
by looking at Table 1, to design a turboprop that could be competitive with respect to 
regional jets, there is the need to increase the seats capacity, as well as reduce the specific 
fuel consumption and the associated emissions. The cabin comfort must be improved by 
reducing the perceived noise and maximum aerodynamic cruise efficiency, the cruise speed 
value must be higher than the currently flying turboprop, and maximum lift coefficients 
at take-off and landing must be increased to be attractive with respect to regional jets in 
terms of ground performance. 
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Table 1. Topic Level Aircraft Requirements (TLARs). 
Performance Requirements 
Passenger Capacity 130 pax at 32” 
1600 NM at Max payload 
Range (108 kg per pax), usual reserve (alternate + holding + 
reserves) 

Cruise Speed At FL300 ISA at 97% MTOW Mach (0.64—0.68) 

Time to Climb At MTOW-ISA-from 1500 ft to FL250 < 16 min 
Take off Field Length (TOFL) ISA-SL at MTOW, TOFL < 4600 ft (1400 m) 
Landing Field Length (LFL) ISA-SL at MLW, LFL < 4600 ft (1400 m) 


Starting from the work shown by Goldsmith [14] and exploiting the lesson learned 
within the frame of the European project IRON [7-12], the authors made several considera- 
tions about different aircraft configurations that can be drafted in facing the design of a 
large turboprop aircraft. 

The increase of the number of passengers from 70 to 130-150 leads to heavier aircraft, 
which means having a larger wing area and increased wingspan, affecting both the size and 
weight of the landing gears, as well as their positions to keep similar ground performance. 
According to the ICAO/Federal Aviation Administration (FAA) [20] span category, for 
aircraft in this class, the wingspan is limited to 36.0 m. Moreover, to ensure ground stability 
requirements due to a heavier aircraft with a larger wingspan, a larger wheel track is 
required as well. 

A conventional high-wing configuration can accommodate large propellers and could 
be a possible solution for a high-capacity turboprop. However, for the high-wing configura- 
tion, two landing gear installations are suitable, as shown in Figure 4: nacelle-mounted or 
fuselage-mounted within fuselage pods. These two solutions may lead to several issues if 
used for high-capacity turboprops. The nacelle-mounted landing gears would require very 
long and heavy legs, which may be difficult to retract inside the nacelles, while fuselage- 
mounted landing gear might require quite large and heavy pods to ensure an adequate 
wheel track. Moreover, for both these layouts, the length of the landing gear legs must be 
increased to achieve a sufficient value of the upsweep angle due to a longer fuselage. 





Figure 4. Landing gear installation solution for a large propeller driven aircraft. 


A low-wing configuration could solve this issue. However, as shown by the example 
in Figure 5, a high-capacity turboprop will need a large propeller diameter to comply with 
an increased thrust requirement due to a higher weight (about 12-14 ft according to the 
estimations made within the EU C52 IRON project). Thus, it would be impossible to ensure 
the required engine clearance from the ground. 

For these reasons, a new turboprop aircraft configuration should be characterized 
by a rear-mounted engine installation. This can lead to a more efficient wing, thanks to 
the possibility to extend the laminar flow region and reduce the cabin noise by installing 
engines far from the cabin. Conversely, such an engine installation would lead to a 
wing weight increment (no mass relief on the wing loading due to the weight of the 
engines). Furthermore, this configuration can provide for a very large center of gravity 
(CG) excursion, which can also affect aircraft performance. A wide range of the CG 
excursion could imply a very large horizontal tail to trim the aircraft in the most rearward 
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CG positions, resulting in a reduction of maximum lift capabilities, while, at the most 
forward CG position, the longitudinal Static Stability Margin (S.S.M.) could be very high, 
requiring for a large download on the horizontal tail to trim the aircraft. This would reduce 
the aerodynamic cruise efficiency, affecting the fuel consumption and the aircraft DOC. 
One possible solution to mitigate this aspect consists in a reasonable cut of the overall CG 
excursion (by limiting the aircraft operability with low passenger number), complying 
with typical aircraft missions. 





Figure 5. Large propeller installation issues for a low wing aircraft configuration. 


Given all the above considerations, by looking at the feasibility study performed by 
NASA [14], and thanks to the experience matured within the IRON project, the authors 
identified in the rear-engine installation the most promising layout. 

Following the work of NASA, the authors also investigated a conventional config- 
uration with engines installed on the wing, with the difference that, in this research, the 
configuration with the upper wing engine installation was discarded because of both 
propeller and nacelle sizes. This layout was replaced by a conventional high-wing configu- 
ration with engines installed under the wing, assuming that the main landing gear could 
be installed in fuselage pods. Another difference with the feasibility study performed by 
NASA is that the T-tail configuration with the rear engine installation was replaced by a 
three-lifting-surface layout. 

The T-tail layout could ensure an adequate longitudinal static stability, as well as 
cruise aerodynamic efficiency, while probably providing for poor high-lift capabilities and 
balance issues due to a very afterward mass concentration. Moreover, the installation of a 
large propeller (12-14 ft of diameter) would require a pylon with a non-negligible span, 
making the latter comparable to a tailplane. 

This configuration was discarded from the first stage of design and replaced by a three- 
lifting-surface layout which could ensure the required longitudinal stability potentially 
increasing as well as both aerodynamic efficiency and high-lift capabilities. 

The introduction of a third lifting surface could help reduce the tail download required 
to trim the aircraft. Thus, the trimmed maximum lift coefficient (Cy. max) could be increased 
while the trim drag contribution could be reduced, resulting in a higher value of the cruise 
aerodynamic efficiency. The potential increment of the CLmax could also lead to a reduced 
wing area providing for additional benefits on cruise aerodynamic efficiency as well as a 
reduction in aircraft weight (unless the third lifting surface would be heavier than the wing 
weight reduction). 

Furthermore, the third lifting surface could allow to shift forward the aircraft neutral 
point position, as well as giving the possibility to optimize both the center of gravity 
excursion and the neutral point position, at fixed static stability margin, by changing wing, 
horizontal tail, and canard sizes together and their relative positioning. In short, the three 
lifting surfaces would extend the design space, allowing to cope with both aerodynamic 
cruise efficiency and maximum lift capabilities of the aircraft. In designing a three-lifting- 
surface aircraft, particular attention must be paid to evaluating the negative effects in terms 
of downwash introduced by the canard that acts on both the main wing and the horizontal 
tail, as shown in [9,11,12]. 
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3. Investigation Workflow and Methods 


Starting from the TLARs described in Table 1, for each of the considered aircraft config- 
urations, a sensitivity analysis was carried out by changing the following design parameters: 


e Wing planform (span, aspect ratio, sweep angle, and longitudinal wing position) 
e Horizontal tail plane planform (span, aspect ratio, and sweep angle) 
e Canard planform if present (span, aspect ratio, and longitudinal position) 


The complete picture of the design space is shown in Table 2, where, for each of the 
previous design variables, the considered range of variation is indicated. A specific combi- 
nation of the design parameters shown in Table 2 generates an aircraft that underwent to a 
Multi-Disciplinary Analysis process (MDA). The MDA involves all major disciplines of the 
preliminary aircraft design (i.e., weight, balance, aerodynamics, performance, and DOC). 
The single MDA cycle shown in Figure 6 and the full factorial Design of Experiments (DOE) 
represented in Table 2 were performed by means of JPAD (Java tool chain of Programs 
for Aircraft Design), software developed at the Industrial Engineering Department of the 
University of Naples Federico II [21-23]. 

JPAD includes several interconnectable analysis modules which can be used both 
in standalone mode and in a complete MDA cycle. One of the most important analysis 
modules is the one dedicated to the evaluation of the most important aircraft performance, 
which implements a smart simulation-based approach to analyze the complete mission 
profile including on-ground phases such as take-off and landing [24]. 

The analysis starts with a first estimation of the amount of fuel needed for the specified 
mission. Then, a balance analysis is carried out to determine the center of gravity excursion. 
For each center of gravity, the aerodynamic and stability module estimates the trimmed 
aircraft lift curve and the trimmed drag polar curve according to the following flight 
conditions: take off, climb, cruise, and landing. Finally, the performance module uses these 
data to make a detailed simulation of the initial mission profile estimating a new amount 
of fuel required to fulfill the mission. Thus, an iterative process is carried out until the first 
estimated fuel mass is equal to the one calculated by the mission profile analysis. Within 
the fuel mass iterative loop shown in Figure 6, a second nested iteration is carried out to 
make the take-off field length and cruise Mach number comply with the assigned TLARs. 
This is obtained by scaling the reference static thrust of all engines with a step of £2.5%. 
Together with this thrust scaling, engine weight is updated as well so that an intermediate 
weight convergence loop must be carried out as well. Both fuel mass and performance 
iterative loops have a maximum number of iterations fixed to 50. All configurations that 
do not converge within this limitation are supposed to be unfeasible so that a penalty will 
be added to those response surface points during the optimization process. 

At the end of the convergence loop, other checks are performed on the analyzed 
configuration. These concern the operative maximum aft center of gravity position, which 
must be at least 5% of the wing Mean Aerodynamic Chord (MAC) more forward than the 
main landing gear position due to rotation issues, and the maximum fuel tank capacity, 
which must be greater than the estimated total mission fuel mass plus reserves. 

The multidisciplinary optimization was carried out by imposing the following constraints: 


e The minimum static stability margin (with the max aft. center of gravity position) 
must be greater than zero. 
Take-off and landing field lengths must be at most 1400 m. 
The time to climb from 1500 to 25,000 ft must be at most 16 min. 
The maximum cruise Mach number should be higher than 0.64 and should not exceed 
0.68 to consider for propeller operative limitation in terms of compressibility issues. 
e The distance between the main landing gear position (Xic) and the maximum aft cen- 
ter of gravity position (XcG max aft) must be greater than 5% of the mean aerodynamic 
chord to avoid ground stability issues. 
e The estimated total mission fuel mass must be lower than wing fuel tank 
maximum capacity. 
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Table 2. Variables design space for each aircraft configuration. 
Pacinicters Wing Mounted Rear Mounted Three Lifting 
Engines Configuration Engines Configuration Surfaces 
Wing Design Parameters 
Span bw (m) 32.0-35.5 32.0-35.5 32.0-35.5 
Longitudinal position XLE w (m) 13.0-16.0 18.0-22.0 19.0—22.0 
Aspect Ratio ARw 9.8-13.3 10.2-13.8 10.2-13.6 
Sweep at leading edge ALE w (°) 2.5-10.0 2.5-10.0 2.5-10.0 
Horizontal Tailplane Design Parameters 
Span by (m) 9.14-11.17 11.7-14.3 11.7-14.3 
Longitudinal position Xz y (m) 37.15 fixed value 30.0-32.0 30.0-32.0 
Aspect Ratio ARy 4.96 fixed value 4.40 fixed value 4.40 fixed value 
Sweep at leading edge ALEH (°) ALE,W + 5.0 ALE,W + 5.0 ALE,W +5.0 
Canard Design Parameters 
Span bc (m) n.a n.a 7.2-8.8 
Longitudinal position Xzc (m) n.a n.a 5.0—8.0 
Aspect Ratio ARc n.a n.a 5.97 fixed value 
Sweep at leading edge ALEc (°) n.a n.a 10.0 fixed value 
Total number of analyzed aircraft O72 6075 34,020 


Aircraft geometry 
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Figure 6. JPAD MDA workflow including thrust update inner loop. 





Mission fuel and static thrust iterative loop 






All analyses were performed under several assumptions concerning masses and 
positions of main components, on-board systems, the center of gravity position to be 
used for each mission phase, fuel tank systems, laminar flow, and downwash effects 
(especially in the case of three-lifting-surface configuration). Fuselage length and maximum 
diameter were kept constant for all analyzed aircraft. Horizontal tailplane sweep angles 
at leading edge were assumed to be 5° greater than the related wing sweep angle at 
leading edge. Components and on-board systems masses were calculated using equations 
provided by Torenbeek [25], which indicate several corrections to be made according 
to the specific aircraft configuration. Center of gravity positions of the main aircraft 
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components were assumed according to Torenbeek [25] as well. Some assumptions were 
made concerning major on-board systems longitudinal position, as shown in Table 3. An 
example of aircraft components and on-board systems center of gravity position is shown 
in Figure 7 considering the case of a three-lifting-surface configuration. 


Table 3. On-board systems assumptions for the MDA cycle of each aircraft configuration. 


Assumption Concerning Systems and Equipment Positions 


Auxiliary Power Unit (APU) 
Air Conditioning and Anti-Icing System 
Electrical System 
Furnishings and Equipment 


87% fuselage length 
30% MAC ahead the wing root leading edge 
42% fuselage length (half cabin) 
37% fuselage length 
Using lifting surfaces positions weighted average 
with components masses 
8% fuselage length (cockpit) 
21% fuselage length (close to the front galley) 


Control Surface Actuators 


Instruments and Navigation System 
Operating Items 


WING cg 
CANARD cg 
HORIZONTAL TAIL cg 
VERTICAL TAIL cg 
ENGINE cg 
FUSELAGE cg 
NOSE GEAR cg 
MAIN GEAR cg 
MEW cg 

OEW cg 

FUEL TANK cg 
---MAC wing 

-e*- Wing root 


FURNISHINGS and EQUIPMENT 


ELECTRICAL CONTROL SURFACES ACTUATORS 


° 


AUXILIARY POWER UNIT (APU) 





AIR CONDITIONING 
FUSELAGE LENGTH 





Figure 7. Example of aircraft main components and systems center of gravity positions for a three-lifting-surface configuration. 


To take into account the aircraft weight reduction due to fuel consumption during the 
design mission, a variable center of gravity position must be considered for each phase, 
affecting in this way the aircraft stability and providing for a variable trim drag value. 
This effect was simulated by considering different center of gravity positions from the 
boarding diagram. For take-off and climb conditions, the maximum take-off weight center 
of gravity was considered. The cruise phase was evaluated assuming a center of gravity 
position related to an average fuel weight condition (starting from the maximum take-off 
weight). The landing condition was investigated assuming a center of gravity position 
related to the maximum landing weight (about 90% of the maximum take-off weight). 
Finally, the longitudinal static stability margin was estimated considering the most aft 
center of gravity position. 

For each of the selected aircraft configurations, some assumptions and some technolo- 
gies were envisaged, as shown in Figure 8 where their effects related to both aerodynamics 
(in terms of drag, maximum lift capabilities, and downwash on the horizontal tail plane) 
and weight estimation are summarized. 
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Figure 8. Technological assumptions for each considered aircraft configuration. 


Wing-mounted engines configuration. This is the state-of-the-art of conventional 
turboprop configuration with a high-wing and T-tail layout. Thus, the horizontal tail-plane 
position was fixed according to vertical tail tip chord position with a fixed incidence angle 
equal to 0.0°. Furthermore, the horizontal tailplane aspect ratio was fixed to about 5.0 with 
the possibility to scale the span and the planform area. Engines positions were linked to 
the starting section of the wing outer panel, assuming laminar flow effects to be active 
only on the outer wing panel. Landing gears were assumed to be mounted in fuselage 
pods. Their positions were linked to the wing position assuming that wing and landing 
gears attachments are applied to the same fuselage frame. Furthermore, landing gear 
pods provide for a higher parasite drag, which was modeled considering an increment 
of about 15 drag counts. According to Torenbeek [25], wing mass were reduced by 5% 
due to wing-mounted engines. To account for landing gear fuselage pod installation, a 
further reduction of 5% of the wing mass was considered. Spoilers’ installation effect on 
wing weight was modeled assuming an increment of the component mass of 2% while, 
for the fuselage, Torenbeek [25] suggested increasing the estimated mass by 8% due to 
pressurization effects and by 7% due to landing gears pods. 

Rear-mounted engines configuration. This innovative configuration provides for 
a low-wing configuration with engines installed at the horizontal tailplane tips. This 
solution would avoid a large engine pylon needed to install a large turboprop engine with 
a propeller diameter of about 12 ft on the fuselage. 

A similar solution is also investigated in [26], where the main configuration criticalities 
dealing with aircraft weight and balance are highlighted. 

From preliminary calculations made by the authors, this configuration will require for 
a large engine pylon affecting the aircraft aerodynamics and interfering with the horizontal 
tail plane. A solution could be the adoption of a T-tail configuration. However, a large 
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pylon would provide an additional lift contribution at high aircraft attitudes leading to a 
larger longitudinal static stability margin with a detrimental effect on maximum trimmed 
lift coefficients. Moreover, a T-tail configuration would provide for a heavier vertical 
tail plane (due to the horizontal tail installation), heavy pylons to hold engines, and a 
heavier fuselage due to additional structural frames required to hold all the rear aircraft 
components. Thus, the authors selected a simplified configuration merging the horizontal 
tail with the pylon. This solution, already investigated within the frame of the European 
funded project named IRON, could yield to a reduced overall mass increment as well as to 
provide for a reasonable value of the aircraft longitudinal stability. The horizontal tail is 
fuselage mounted, and it was assumed to have a variable pitch angle (iz;). The horizontal 
tailplane Aspect Ratio (ARy) was fixed at 4.40 for structural considerations regarding 
the engine installation. Being the considered layout a low-wing configuration, landing 
gears were assumed to be wing-mounted. This will help to limit the required fuselage 
pods size to accommodate landing gear legs since only wheels and a small legs portion 
should be retracted inside them. Thus, the detrimental impact that large fuselage pods 
could have on fuselage weight and aerodynamic drag will be reduced. According to the 
selected main landing gear installation, its position was linked to the wing at 60% of the 
MAC. According to Torenbeek [25], wing mass was increased by 2% to consider for wing 
spoilers. Concerning the fuselage, Torenbeek [25] suggested increasing the estimated mass 
by 8% due to pressurization effects and by 4% due to rear engine installation. Moreover, 
the horizontal tailplane mass was increased by 65% on the base of preliminary evaluations 
of combined aerodynamic loading and engine inertial contributions. Unlike the wing, the 
horizontal tailplane usually works with negative aerodynamic loads, thus the engine mass 
does not provide for load relief effects. Thanks to the rear engine installation, the whole 
wing was supposed to work in laminar flow conditions, reducing the wing parasite drag of 
about 20 drag counts. 

Three-lifting-surface configuration. This configuration, as well as the rear-mounted 
engine aircraft, provides for a rear engine installation at the horizontal tailplane tips. 
As for the previous layout, the horizontal tail was supposed to have a variable pitch 
angle according to the specific flight phase (same values assumed for the rear-mounted 
configuration) and the tail aspect ratio was kept constant at 4.40 due to structural reasons 
linked to the engine installation. The same horizontal tailplane mass increment of 65% was 
assumed due to the same structural considerations made for the rear-mounted engines 
solution. The same assumptions also applied for the main landing gear position, horizontal 
tail mass increment, and laminar flow effects on the wing. The amount of storable fuel 
is linked to the estimated tank capacity calculated assuming standard spar positions and 
using the volume equation proposed by Torenbeek [25]. 


4. Results 


Once the MDA analysis cycle for each aircraft model investigated with the DOE 
was completed, results dealing with some of the most important design objectives were 
gathered. Figure 9a illustrates the complete cloud of points analyzed by the MDA process 
for the three-lifting-surface configuration in terms of block fuel versus maximum take-off 
weight and maximum aerodynamic cruise efficiency. In addition, Figure 9b,c show how 
the constraint on the static stability margin (S.S.M. > 0.0) reduces the number of possible 
candidate solutions in two different planes: block fuel versus maximum take-off weight 
and block fuel versus maximum aerodynamic cruise efficiency. Blue circles represent all 
solutions that are compliant with the imposed constraint on the S.5.M. The same applies 
for each additional constraint included in the optimization problem. 
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Figure 9. Results of design of experiments for the three-lifting-surface: (a) complete cloud of points; 
(b) effect of the 5.S.M constraint on possible candidate solutions in the plane block fuel versus MTOW; 
and (c) effect of the S.S.M constraint on possible candidate solutions in the plane block fuel versus 
cruise maximum efficiency. 
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The optimization problem defined in Equation (1) was mainly focused on the mini- 
mization of the environmental impact provided by each configuration. Thus, the block fuel 
was considered as the only objective function. 


Minimize the Block fuel on a mission range of 1600 NM 

with respect to lifting surfaces planform design variables and relative positioning 

subject to : 

5.S.M > 0.0 

Take — off field lenght & Landing field length < 1400 m (1) 
Climb time from 1500 ft to 25,000 ft < 16 min 

Cruise Mach number 0.64 < M,, < 0.68 


[Xtc — Xcomaxaft|] = 5%MAC 


Estimated mission fuel < Max storable fuel mass 


Although aircraft are typically designed aiming at minimizing the DOC, this quantity 
was not considered inside the objectives set. Main contributions to DOC are due to aircraft 
weight, cruise Mach number, and mission block fuel. However, concerning the aircraft 
weight for a fixed configuration, its variation inside the obtained response surface is mainly 
due to the mission block fuel. Thus, the effect of the aircraft maximum take-off weight on 
the DOC is strictly based on the block fuel variation. In addition, following the lessons 
learned from the EU project IRON, the response surface cruise Mach number variation 
was limited to a narrow range from 0.64 to 0.68. Those limits were driven by market 
requirements for the lower bound and by propeller tip compressibility issues for the upper 
bound. In this range, the effect of the cruise Mach number on the aircraft DOC resulted 
to be of second order with respect to the mission block fuel. For all these reasons and 
considering the additional limitations brought by the imposed set of constraints to the 
DOC value, the latter was not considered as an objective function. However, its value, 
in terms of both total and cash DOC, was monitored to make a comparison between the 
optimized configurations. 

To demonstrate the above, Figures 10—12 show the variation of the estimated total 
DOC with respect to the block fuel for each analyzed configuration, highlighting all feasible 
solutions compliant with the imposed set of constraints. It must be noted that the minimum 
feasible block fuel almost corresponds to the minimum DOC solution. In addition, the 
range of variation of the total DOC value for each aircraft configuration resulted to be 
very limited. 


0.136 


0.134 


>) 
= 
Ww 
NO 


© 
— 
Ww 


© 
= 
N 
(oe) 





© 
— 
NO 
o>) 








Total DOC (US$/seat*NM) 


0.124 











0126506 6000 6500 7000 7500 
Block Fuel (kg) 


Figure 10. Total DOC versus block fuel for the wing-mounted engines configuration. Blue circles 
represent the solutions compliant with the imposed set of constraints. 
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Figure 11. Total DOC versus block fuel for the rear-mounted engines configuration. Blue circles 
represent the solutions compliant with the imposed set of constraints. 
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Figure 12. Total DOC versus block fuel for the three-lifting-surface rear-mounted engines configura- 
tion. Blue circles represent the solutions compliant with the imposed set of constraints. 


Selected optima aircraft models are compared in Figure 13 including a visual represen- 
tation of their center of gravity excursion, while their main characteristics are summarized 
in Table 4. 

Both the wing- and rear-mounted engine aircraft are characterized by a wing area of 
about 105 m* and wing aspect ratio of about 12, while the three-lifting-surface aircraft has 
a reduced wing area (about 100 m7) thanks to the additional lift contribution provided by 
the canard. 

The two configurations with rear-mounted engines present a very large horizontal tail 
area (about 44% of the wing area) due to the need to balance a very large center of gravity 
excursion, as illustrated in the boarding diagrams from Figure 13. 

Nevertheless, the three-lifting-surface configuration presents a higher aerodynamic 
cruise efficiency thanks to the positive effect of the forward lifting surface (the canard) 
on the global trim drag contribution. It is worth noticing that the total aircraft wetted 
area (which is linked to the parasite drag) is slightly higher for the three-lifting-surface 
compared with the rear-mounted engine configuration. 

All three optimal configurations are stable with respect to the most aft center of gravity 
position with a reduced Static Stability Margin (S.S.M.), as shown in Table 4. In the authors’ 
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opinion, the obtained S.S.M. value is acceptable and compliant with an improved flight 
control system, which is forecasted for entry into service 2035. 

Using the main performance and DOC as rules of comparison of the three optimal 
large capacity turboprop configurations, the three-lifting-surface configurations was proven 
to be the best solution among the considered layouts. In fact, although with very similar 
values of ground performance, the three-lifting-surface aircraft provides for the lowest 
value of the block fuel, thanks to the improved cruise aerodynamic efficiency, as well as for 
the lowest maximum take-off weight. 

Despite the additional component, the main reason behind the three-lifting-surface 
configuration reduced max take-off weight is related to the different engine weight coupled 
with a lower amount of block fuel. Focusing on the engine mass, this value was linked to the 
value of the static thrust inside the workflow in Figure 6, where, to match both the required 
take-off field length and cruise Mach number, the static thrust was adjusted iteratively. 
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Figure 13. Selected optima solutions for each aircraft configuration: (a) optimized high wing with under wing-mounted 
engines configuration and its boarding diagram; (b) optimized low wing with rear-mounted engines configuration 
and its boarding diagram; and (c) optimized three-lifting-surface with rear-mounted engines configuration and its 


boarding diagram. 
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Wing Mounted Rear Mounted Three-Lifting- 
Parameters ; ; 

Engines Engines Surface 
Wing area, Sw (m7) 104.4 104.6 101.2 
Horizontal tail area, Sy (m7?) 16.83 46.5 44.82 
Vertical tail area, Sy (m7?) 25.0 25.0 25.0 
Canard area, Sc (m*) - - 9.37 
Wing aspect ratio, ARw 12.07 12.04 12:02 
Horizontal tail aspect ratio, AR 5.00 4.40 4.40 
Vertical tail aspect ratio, ARy 1.37 1.37 £37 
Canard aspect ratio, ARc - - 07 
Fuselage length, |; (m) 38.04 38.04 38.04 
Fuselage diameter, d¢ (m) 3.535 3.535 3-090 
Single engine static thrust, To (lbf) 23,603 26,054 23,027 
Single Engine Mass (kg) 3367 4122 3297 
Max Take-Off Weight (kg) 57 A419 58,794 56,640 
Operating Empty Weight (kg) 35,665 36,820 34,593 
Design Payload (kg) 14,040 14,040 14,040 
Max forward CG (% MAC) 22.7% 8.5% 1.2% 
Max afterward CG (% MAC) 50.6% 53.9% 23.3% 
Aerodynamic cruise efficiency (at 97% MTOW) 16.7 16.6 17.6 
Maximum aerodynamic cruise efficiency 18.4 17.8 18.8 
Static Stability Margin (%) 3.64 1.16 1.95 
Cre 1.63 1.57 1.79 
Ch max, TO 2:09 2.30 2.52 
Cl max, LND 207, 2.86 3.10 
Take-Off Field Length (m) 1396 1380 1380 
Landing Field Length (m) 1339 1384 1336 
Climb Time (1500 ft—25,000 ft at 190 kt) (min) 15.0 13.0 15.7 
Max Cruise Mach Number 0.64 0.67 0.66 
Block time—1600 NM (min) 239 234 257 
Block fuel—1600 NM (kg) 6259 6479 5958 
Total DOC—1600 NM (¢/seat NM) 13.02 12.95 12.83 
Cash DOC—1600 NM (¢/seat NM) 7.68 7.70 7 54 


By comparing the three-lifting-surface configuration with the rear-mounted engines 
configuration, the single engine mass is reduced by about 1000 kg (2000 keg for both en- 
eines), which justifies the difference in the maximum take-off weight and the lower value 
of the three-lifting-surface configuration block fuel. On the other hand, the comparison 
between the three-lifting-surface model and the high wing with wing-mounted engines 
configuration highlights similar static thrust and engine weight values. However, the re- 
duced aerodynamic efficiency of the wing-mounted engines configuration, due to the set of 
assumptions shown in Figure 8, coupled with a heavier fuselage led to higher values of both 
maximum take-off weight and block fuel with respect to the three-lifting-surface model. 

A reduced block fuel mass is also linked to a lower amount of pollutant emissions 
(at fixed engine database); thus, the three-lifting-surface configuration resulted to be the 
ereenest among the investigated configurations. Moreover, thanks to a lighter structure 
(lowest Operating Empty Weight (OQEW)) and the reduced amount of fuel needed for the 
design mission of 1600 NM, this aircraft also provided the lowest total DOC and cash 
DOC. Those are mainly influenced by the block fuel mass, the block time, and the aircraft 
utilization, expressed in terms of block hours per year and calculated as proposed in the 
book by Kundu [27]. 

The reference regional jet aircraft was modeled using the JPAD software considering a 
set of TLARs related to the Airbus A220-300. The main information concerning this aircraft 
was retrieved from several public sources, including the aircraft manual, the European 
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Aviation Safety Agency (EASA) type-certificate data sheets, public aircraft data archives, 
and the Base of Aircraft Data (BADA) database considering aircraft models similar to 
the A220-300 [28-32]. A data summary of the most important aircraft characteristics is 
provided in Table 5. 


Table 5. Main data concerning the A220-300. 





TLAR 

Accommodation (Typical-Full Economy) 135 
Design range (typical) 3100 NM 
Take-Off Field Length (Max Take-Off Weight, ISA conditions, Sea Level) 1890 m 
Landing Field Length (Max Take-Off Weight, ISA conditions, Sea Level) 1509 m 
Cruise Mach number (typical) 0.78 
Cruise altitude (typical) 37,000 ft 
Max cruise Mach number at 37,000 ft 0.82 
Max operating altitude 41,000 ft 
Alternate cruise range (assumed by authors) 200 NM 
Alternate cruise altitude (assumed by authors) 20,000 ft 
Holding duration (assumed by authors) 30 min 
Holding altitude (assumed by authors) 1500 ft/min 
Residual fuel reserve (assumed by authors) 5% 

Geometrical and Operational Data 
Wing area 112.3 m? 
Wingspan 35.1 m 
Wing aspect ratio 10.97 
Fuselage length 38.71 m 
Fuselage diameter 3.7m 
Single engine static thrust 24,400 Ibf 
Engine by-pass ratio 12:1 
Max Take-Off Weight 67,585 kg 
Max Landing Weight 58,740 kg 
Max Zero-Fuel Weight 55,792 kg 
Operating Empty Weight 37,081 kg 
Max Payload 18,711 ke 
Max Fuel Mass 17,726 ke 
BADA averaged climb speed (CAS) 271 kt 
BADA averaged rate of climb 1642 ft/min 
BADA maximum rate of climb 2862 ft/min 
BADA averaged descent speed (CAS) 218 kt 
BADA averaged rate of descent 2186 ft/min 
BADA maximum rate of descent 3700 ft/min 


A complete case study concerning this aircraft model and its implementation inside 
the JPAD software was carried out by authors and reported in [33]. Here, starting from 
publicly available data concerning the aircraft geometry and using the aircraft 3-views for 
the geometry digitization process, a parametric model was generated. The latter was used 
to perform a complete multi-disciplinary analysis cycle, as described in Figure 6, neglecting 
the static thrust update loop but including the mission fuel feedback loop. 

It must be noted that, in terms of DOC, the fuel price value needed to perform the 
analysis was assumed according to IATA fuel price monitor [34], while aircraft price was 
assumed considering reference data reported in [35]. Engine unit costs were deducted 
from [36] as well. 

The main output of the JPAD multi-disciplinary analysis was compared with aircraft 
data reported in Table 5 to validate the calculation case. A summary of this comparison is 
reported in Table 6, while a comparison in terms of payload—range diagram is provided 
in Figure 14, where the main differences between the two charts can be addressed to the 
gap between the values of the Max Zero-Fuel Weight (MZFW), the OEW, and the max fuel 
mass of the A220-300 with respect to the values calculated by JPAD. 
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Table 6. Comparison between JPAD output and A220-330 data in Table 5. 














Parameters JPAD A220-300 (Table 5) Difference (%) 
Max Take-Off weight (kg) 66,911 67,585 —1.00% 
Max Landing weight (kg) 56,875 58,740 —3.18% 
Max fuel mass (kg) 17,233 17,726 —2.78% 
Max Zero-Fuel Weight (kg) 55,017 55,792 —1.39% 
Operating Empty Weight (kg) 36,306 37,081 —2.09% 
Take-Off Field Length (m) 1814 1890 —4,02% 
Landing Field Length (m) 1575 1509 +4.37% 
60,000 
55,000 
“Bo JPAD 
£50,000 pa-=s====95555252252522522=5 } A220-300 
S | aan Design Condition JPAD 
S 45,000 | ----- Design Condition A220-300 
° 
o OEW (JPAD) = 36,004 kg 
~ 40,000 ! 
35,000 ' 


O 500 1000 1500 2000 2500 3000 3500 4000 
Range (NM) 


Figure 14. Payload—range chart calculated with JPAD compared with the A220-300 Payload-Range 
chart. Adapted from ref. [28]. 


Once the A220-300 JPAD model was validated, it was used to perform an additional 
simulation considering the design range from the set of TLARs in Table 1. In addition, to 
compare this reference regional jet with all three optimal turboprop aircraft layouts, its 
payload mass was modified to comply with the design payload of 14,040 kg used for each 
high-capacity turboprop configuration. 

Finally, a comparison between the three optimal configurations and a reference re- 
gional jet platform, assumed as the Airbus A220-300, is presented in Table 7 considering a 
mission of 1600 NM. 

As shown in Table 7, all high-capacity turboprop configurations provide a beneficial 
effect on the amount of block fuel related to the design mission of 1600 NM. However, it 
must be noted that the considered reference regional jet is designed for a mission range of 
3300 NM, thus the aerodynamic efficiency and the take-off weight are not optimized for a 
1600 NM mission range. Although this seems to be an unfair comparison, it must be high- 
lighted that this study aims to compare the benefit in terms of fuel consumption that a new 
large capacity turboprop, specifically designed for short/medium haul, would have with 
respect to the existing state-of-the-art RJ also currently operated on this kind of missions. 

In addition to this comparison, the authors also designed a possible new regional 
jet aircraft using the same set of TLARs of the turboprop aircraft. Only the cruise Mach 
number was shifted from 0.64 to 0.78 since this is a jet-driven aircraft. This comparison, 
shown in Table 7, was added to demonstrate that, even if a regional jet specifically designed 
on such a mission range existed, the proposed large turboprop aircraft would still provide 
benefits in terms of environmental impact. The major drawback of the turboprop solution 
is represented by a higher total DOC value mainly due to the reduced block speed. 
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Table 7. Comparison with the reference regional jet in terms of block time, block fuel, and DOC (percent differences, in 


brackets, with respect to the regional jet). 


Design Mission: 


Wing Mounted 


Rear Mounted 


Three-Lifting- 


RJ Designed 








1600 NM Engines Engines Surface at 1600 NM teen) 

Take-off weight (kg) 57,419 (—3.11%) 58,794 (—0.79%) 56,640 (—4.42%) 54,378 (—8.24%) 59,260 
Cruise Mach number 0.64 (—17.95%) 0.67 (—14.10%) 0.66 (—15.38%) 0.78 (+0.0%) 0.78 
Aerodynamic Cruise efficiency 16.7 (0.0%) 16.6 (—0.60%) 17.6 (+5.38%) 17.2 (+2.99%) 16.7 

Relative cruise SFC (w.r.t. RJ) —19.22% —17.02% —18.56% - - 

Utilization (h/year) 3353 (+0.60%) 3346 (+0.39%) 3350 (+0.51%) 3329 (—0.12%) 3333 
Block Time (min) 239 (+6.22%) 234 (+4.00%) 237 (+5.33%) 222 (—1.33%) 225 
Block Fuel (kg) 6259 (—10.40%) 6479 (—7.24%) 5958 (—17.24%) 6569 (—5.96%) 6985 

Total DOC (¢/seat*NM) 13.02 (—13.37%) 12.95 (—13.84%) 12.83 (—14.63%) 12.03 (—19.96%) 15.03 
Cash DOC (¢/seat*NM) 7.68 (—7.69%) 7.70 (—7.45%) 7.54 (—9.37%) 8.32 (+0.00%) 8.32 


The design of this new regional jet was accomplished using the same multi-disciplinary 
analysis and optimization workflow adopted for all proposed turboprop configurations. 
Starting from the TLARs shown in Table 1 and considering as baseline the JPAD paramet- 
ric model of the A220-300, the set and ranges of design parameters reported in Table 8 
were considered. 


Table 8. Variables design space for the regional jet model designed on a 1600 NM mission profile. 


Parameters Wing Design Parameters 
Span by (m) 32.0-35.1 
Longitudinal position XLz w (m) 10.0-13.0 
Aspect Ratio ARw 9.10-12.19 
Sweep at leading edge ALE w (°) 25.0-35.0 


Horizontal Tailplane Design Parameters 





Span by (m) 10.75-13.14 
Longitudinal position XLE Hy (m) 31.6 fixed value 
Aspect Ratio ARy 5.12 fixed value 


Sweep at leading edge ALE H (°) ALEw + 95.0 


The optimization problem carried out for this aircraft is represented in Equation (2) 
where the block fuel was assumed as the only objective function. It must be noted that the 
optimization problem is the same as the one considered for all turboprop aircraft, except for 
the time to climb condition and the cruise Mach number, which was considered to change 
in the range from 0.78 to 0.82 corresponding to the typical and maximum Mach numbers of 
the A220-300, respectively. The previous climb time condition was not included in the set 
of constraints since this aircraft is supposed to operate at a higher cruise altitude (37,000 ft 
instead of 30,000 ft). 


Minimize the Block fuel on a mission range of 1600 NM 

with respect to lifting surfaces planform design variables and relative positioning 

subject to : 

S.S.M > 0.0 

Take — off field lenght & Landing field length < 1400 m (2) 
Cruise Mach number 0.78 < M,, < 0.82 


[Xtc —Xcomaxaft|] = 5%MAC 
Estimated mission fuel < Max storable fuel mass 
The MDA workflow adopted for the analysis of each regional jet models included in 
the design space was the same as the one shown in Figure 6, where the engine static thrusts 


were modified iteratively to match the requirements in terms of cruise Mach number and 
take-off field length. However, the rubberized engine database used for the performance 


Aerospace 2021, 5, 132 


21 of 26 


simulation was the same used to analyze the A220-300 parametric model, thus the same 
SFC was considered at fixed altitude and cruise Mach number. 

To make a fair comparison with the optima turboprop solutions, the selection of the 
optimum regional jet aircraft was driven by the minimization of the block fuel targeting 
the lowest environmental impact. The DOC was estimated only for comparison purposes. 

A geometrical comparison between the baseline regional jet model and the selected 
optimum jet aircraft is shown in Figure 15. 


20.0 


“iy - A220-300 
16.0 - Optimized RJ at Range 1600 NM 


12.0 






8.0 


4.0 


X(m) 
0.0 


0.0 5.0 10.0 LU 20.0 25.0 30.0 35.0 40.0 45.0 


Figure 15. Comparison between top views of the baseline regional jet model (red) and the selected 
optimum regional jet aircraft (black). 


Starting from the results shown in Table 7, a first comparison that can be made is 
between the proposed turboprop solutions and the existing state of the art regional jet (the 
last column of Table 7) typically operated on short/medium haul with a seats capacity 
similar to the required value of 130. 

Focusing on each turboprop platform, the high-wing configuration with wing-mounted 
engines has the greatest advantage in terms of mean cruise SFC with a reduction of —19.22% 
and provides for the second best reduction in take-off weight (—3.11%). However, the 
cruise aerodynamic efficiency shows no improvements, with respect to the reference re- 
gional jet, and the cruise Mach number is the lowest among all analyzed turboprop aircraft 
models. Thus, this configuration is the worst in terms of block time providing for an 
additional little negative effect on the block fuel. 

From the combination of all these effects, this aircraft may be a good candidate to 
challenge a regional jet on this kind of missions with an overall block fuel reduction of 
0.40% and only 15 min of additional block time. In term of DOC, good results were reached 
as well, with a reduction of —13.37% in total DOC and —7.68% in cash DOC. However, 
despite the beneficial effect provided by the smaller block fuel mass, this aircraft is the 
worst in terms of both utilization and block time, leading to the lowest reduction in total 
DOC. On the other hand, being the effect of the block fuel much more effective on the cash 
DOC rather than on the total DOC, this configuration provides for a slightly lower cash 
DOC value with respect to the one with low-wing and rear-engine installation. 

The second high-capacity turboprop configuration showed the lowest reduction in 
term of block fuel (—7.24%), with respect to the reference regional jet, since the beneficial 
effect provided by the smaller value of the mean cruise SFC (—17.02%) is mitigated by 
a very little reduction in MTOW (—0.79%) and a lower value of the cruise aerodynamic 
efficiency (—0.60%). On the other hand, the cruise Mach number was the greatest among 
all analyzed turboprop configurations, leading to the lowest increment in block time with 
respect to the reference regional jet. However, being this effect very limited, the related 
effect on the block fuel is also very limited. In terms of DOC, comparing this configuration 
with all other turboprop platforms, the lower block time provides for a better value of the 
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utilization parameter which positively affects all costs contributions. Major drawbacks are 
provided by the increased amount of block fuel (and so of the fuel price) and by the take-off 
weight (the highest among all turboprop configurations) which influences maintenance 
costs. However, those effects are smaller than the one provided by the utilization, making 
this aircraft the second best solution in terms of total DOC reduction with respect to the 
reference regional jet. As for the previous case, the effect of the block fuel on the cash 
DOC is much more effective than on the total DOC, providing for the worst reduction of 
this value. 

As stated above, the three-lifting-surface configuration resulted to be the best tur- 
boprop configuration. This was further confirmed by the comparison with the reference 
regional jet, from which this innovative high-capacity turboprop aircraft has highlighted 
the greatest reduction both in term of block fuel and DOC. In particular, this aircraft pro- 
vided a mean cruise SFC reduction similar to the first configuration (—18.56% instead of 
—19.22%) together with the lowest value of the take-off weight and the highest value of 
the cruise aerodynamic efficiency (—4.42% and +5.38%, respectively, if compared to the 
reference regional jet). A cruise Mach number slightly lower than the second configuration 
led to intermediate values of both block time and utilization between the first two analyzed 
aircraft models. 

A second comparison that can be made is between the proposed turboprop solutions 
with the regional turboprop designed according to the same set of TLARs. 

Considering that this aircraft was optimized on the same set of TLARs adopted for 
each turboprop aircraft, the amount of block fuel needed to operate the 1600 NM design 
mission is obviously reduced with respect to the one related to the A220-300 parametric 
model. Although with block fuel values in line with some of the analyzed optima turboprop 
models, the three-lifting-surface configuration still provide for a beneficial effect in terms 
of environmental impact with respect to this possible regional jet model achieving a fuel 
saving of about 611 kg (—9.30%). 

Conversely, the new regional jet represents the best solution in terms of total DOC 
having both an improved block time with respect to the optima turboprop aircraft and the 
lowest take-off weight. In particular, by comparing the three-lifting-surface configuration 
with this newly designed regional jet, there is a total DOC increment of about +6.65%. 
However, dealing with the cash DOC, which is not influenced by the beneficial effect pro- 
vided by a reduced take-off weight, the three-lifting-surface configuration allows achieving 
the same cost saving shown with respect to the A220-300 model of about —9.37%. Thus, 
the proposed three-lifting-surface configuration represents the best solution in terms of 
both environmental impact and cash DOC. 


5. Conclusions 


This paper considers the possible design of a modern high-capacity turboprop aircraft 
that could potentially reduce the environmental impact of regional aircraft with respect to 
the current state-of-the-art regional jet widely adopted on short/medium hauls. The study 
investigated three possible architectural solutions for a large capacity turboprop aircraft. 
A multi-disciplinary analysis and optimization process was carried out to assess the best 
solution complying with the imposed design space boundaries and operative constraints. 

By comparing the results obtained for each optimum turboprop aircraft with respect 
to a reference regional jet model, represented by the Airbus A220-300 operated on a range 
of 1600 NM, this study highlighted a maximum potential fuel saving of about —17.24% 
considering a three-lifting-surface configuration. 

In addition, to make a fair comparison between the analyzed turboprop aircraft 
platforms and a regional jet, a brand-new turbofan aircraft was designed using the same set 
of TLARs used for the turboprop configurations. It must be noted that such a regional jet 
aircraft, coupling a design range of 1600 NM with a seat capacity of 130 passengers, does 
not exist in the current regional market scenario. By comparing the best analyzed turboprop 
configuration with this regional jet model, the initial block fuel saving of 7.24% reaches a 
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lower value of —5.96%, still proving that a modern optimized turboprop platform could 
provide for a beneficial effect in terms of environmental impact. 

For the sake of completeness, a comparison between all optimal turboprop configu- 
rations and both the analyzed regional jet models was carried out also in terms of DOC. 
Starting with the existing reference regional jet model, being the latter not optimized on 
a 1600 NM mission range, all turboprop models provide for a cost reduction in terms 
of both total and cash DOC (max reduction of about —14.63%). However, considering a 
regional jet specifically designed for such a mission range, it can be highlighted that the jet 
aircraft represents the best solution in terms of total DOC (—6.65% with respect to the best 
analyzed turboprop configuration) but not in terms of cash DOC due to the higher block 
fuel (+9.37% with respect to the best analyzed turboprop configuration). 

Since this work is not focused on the investigation of the best possible turboprop 
aircraft architecture, implementing future technologies or layouts, the performed investiga- 
tion only considered three possible configurations of a large capacity turboprop assumed 
to be compliant with an entry into service by the year 2035. A wider study, extended to a 
larger portfolio of innovative, or even disruptive, configurations should be performed to 
find a better solution with respect to the one proposed in this paper. 
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Abbreviations 
The following abbreviations are used in this manuscript. 


APU Auxiliary Power Unit 

ATAG _ Air Transport Action Group (www.atag.org) 

ATR Avions de Transport Régional (www.atr-aircraft.com) 
BADA | Base of Aircraft Data (www.eurocontrol.int/model/bada) 
CAS Calibrated Air Speed 

CG Center of Gravity 

CS2 Clean Sky 2 (www.cleansky.eu) 

DOC Direct Operating Cost 

DOE Design of Experiments 


EASA — European Aviation Safety Agency (www.easa.europa.eu) 

EU European Union 

FAA United States Federal Aviation Administration (www.faa.gov) 
FL Flight Level 
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IATA 
ICAO 
IRON 
ISA 
JPAD 
LEL 
LG 
MAC 
MDA 
MDAO 
MLW 
MTOW 
MZEW 
NASA 
NLF 
OEW 
RJ 
S.S.M. 
SFC 
TLARs 
TOFL 


Symbols 
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International Air Transport Association (www.iata.org) 
International Civil Aviation Organization (www.icao.int/Pages/default.aspx9) 
Innovative turbopROp configuration 

International Standard Atmosphere 

Java tool chain of Programs for Aircraft Design 
Landing Field Length 

Landing Gear 

Mean Aerodynamic Chord 

Multi-Disciplinary Analysis 

Multi-Disciplinary Analysis and Optimization 
Maximum Landing Weight 

Maximum Take-Off Weight 

Max Zero-Fuel Weight 

United States National Aeronautics and Space Administration (www.nasa.gov) 
Natural Laminar Flow 

Operating Empty Weight 

Regional Jet 

Static Stability Margin 

Specific Fuel Consumption 

Top-Level Aircraft Requirements 

Take-Off Field Length 


The following symbols are used in this manuscript. 


¢ 

ARec 
ARy 
ARy 
ARw 

be 

by 

bw 
CLoH 
Tmax 
CLmax,LND 
Chmax,TO 
Coz 

dcs 

de 


KeGmacat 
XLE,C 
XLE,H 
XLE,W 

XLG 


United States dollar cent 

canard aspect ratio 

horizontal tailplane aspect ratio 

vertical tail aspect ratio 

wing aspect ratio 

canard span 

horizontal tailplane span 

wingspan 

horizontal tail lift curve slope 

maximum lift coefficient 

maximum lift coefficient at landing 

maximum lift coefficient at take-off 

carbon dioxide 

drag counts 

fuselage diameter 

horizontal tailplane root incidence angle 
fuselage length 

Mach number 

cruise Mach number 

nautical mile 

nitrogen dioxide 

canard area 

horizontal tailplane area 

vertical tailplane area 

wing area 

single engine static thrust 

United States dollar 

canard sweep angle at the leading edge 
horizontal tailplane sweep angle at the leading edge 
wing sweep angle at the leading edge 

center of gravity longitudinal position 

max aft center of gravity longitudinal position 
canard root leading edge longitudinal position 
horizontal tail root leading edge longitudinal position 
wing root leading edge longitudinal position 
main landing gear longitudinal position 
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Abstract: This paper describes the relative frequency of reports of oil and hydraulic fluid fumes in the 
ventilation supply air (“fume events”) compared to other types of fumes and smoke reported by U.S. 
airlines over 10 years. The author reviewed and categorized 12,417 fume/smoke reports submitted 
to the aviation regulator to comply with the primary maintenance reporting regulation (14 CFR § 
121.703) from 2002-2011. The most commonly documented category of onboard fumes/smoke was 
electrical (37%). Combining the categories of “bleed-sourced”, “oil”, and “hydraulic fluid” created 
the second most prevalent category (26%). The remaining sources of onboard fumes/smoke are 
also reported. To put the data in context, the fume event reporting regulations are described, along 
with examples of ways in which certain events are underreported. These data were reported by U.S. 
airlines, but aviation regulations are harmonized globally, so the data likely also reflect onboard 
sources of fumes and smoke reported in other countries with equivalent aviation systems. The data 
provide insight into the relative frequency of the types of reported fumes and smoke on aircraft, 
which should drive design, operational, and maintenance actions to mitigate onboard exposure. The 
data also provide insight into how to improve current fume event reporting rules. 


Keywords: aircraft; fume event; reporting; engine oil; hydraulic fluid; electrical 


1. Introduction 


The majority of commercial and military aircraft are designed to extract (or “bleed”) 
the cabin and flight deck ventilation air from an engine compressor located upstream of the 
cabin and flight deck, whether the compressor is in the main aircraft engines or an auxiliary 
engine in the tail of the aircraft called the auxiliary power unit (APU). The problem with this 
“bleed air” design is that the oil that lubricates the engine compressors can contaminate the 
ventilation supply air. This can be caused by an actual mechanical failure or by the effects of 
engine parts subjected to physical and thermal stress during engine power setting changes, 
improper oil drainage/reingestion, or overfilling of an oil reservoir during servicing [1—4]. 
Additionally, small amounts of oil seep past the engine/ APU bearing chamber seals into 
the compressor air flow, which is considered normal and acceptable if the volume of oil 
loss over time is within defined limits [1-4]. When oil contaminates the bleed air stream, 
the fumes are delivered to the cabin and flight deck located downstream of the compressor. 
The APU (and, less often, the main engines) can also ingest hydraulic fluid that has spilled 
or leaked from a burst line, a worn part, or overfilling during servicing. The ingested 
hydraulic fluid is also heated in a compressor and, likewise, contaminates the ventilation 
air supplied to the flight deck and cabin. Engine oil and hydraulic fluid fumes contain 
a complex mixture of compounds, including base stocks, additives, and decomposition 
products generated upon exposure to heat and moisture [1,5]. Because current systems are 
not equipped with bleed air filters, these contaminants are mixed with the ventilation air 
and delivered to the occupied zones through the cabin and flight deck air supply vents. 

Of the compounds that can contaminate the aircraft air supply system, oil and hy- 
draulic fluid have received the most attention because of decades of reports of compromised 
flight safety and security when crewmembers are either impaired or incapacitated by acute 


https://www.mdpi.com/journal/aerospace 


Aerospace 2021, 5, 122 


2 of 14 


symptoms inflight [6-9]. Additionally, these “fume events” are associated with operational 
impacts, such as diversions, cancellations, and repeat events [10,11], as well as reports of 
chronic neurological and respiratory illness [12,13]. 

However, there are other sources of supply air contamination onboard aircraft, includ- 
ing engine exhaust, deicing fluid, birds, engine wash, and ambient pollution. Additionally, 
there are sources of airborne chemical contaminants inside the cabin and flight deck, in- 
cluding electrical systems (e.g., lights, heaters, doors, and entertainment systems), fans 
(e.g., ventilation recirculation, cooling blowers, etc.), ducting (blown or clogged), and 
defective batteries. 

Globally, within the aviation industry, there is discussion and debate over the rel- 
ative frequency and nature of the different types of onboard fumes/smoke reported by 
crewmembers to airlines and by airlines to regulatory agencies. The answers to these 
questions of “how often” and “what type” of onboard fumes are important because they 
directly influence the debate regarding what design, operational, and maintenance reg- 
ulations are needed (if any) to mitigate exposure to fumes. Some reports recommend 
engineering control measures and worker training programs to prevent onboard exposure 
to engine oil and hydraulic fluid fumes to mitigate the associated flight safety and crew 
health hazards [14-16]. Other reports downplay the frequency and seriousness of engine 
oil and hydraulic fluid fume events, acknowledging that they occur, but referring to them 
as “smell events” and claiming that they are either “rare” or “extremely rare” [17-20], 
rendering additional engineering and maintenance control measures unnecessary. 

The aviation regulator in the United States (Federal Aviation Administration, FAA) 
has acknowledged that airlines underreport fume events [21] and that some of the fume 
event data that airlines submit is not reported properly [22]. However, even the number 
of events that airlines document with the FAA significantly underrepresents the actual 
number of events because airlines are not required to report either events that occur during 
eround operations or events for which fumes or smoke were present, but a mechanical 
fault was not subsequently identified. Further, the FAA has significantly undercounted the 
number of oil and hydraulic fume event reports it received from U.S. airlines [23]. 

To identify the most effective means to mitigate exposure to onboard fumes and smoke, 
it is helpful to start by assessing the relative frequency of the sources of onboard fumes 
and smoke that are documented. 


2. Materials and Methods 


On 14 March 2014, the author filed a Freedom of Information Act (FOIA) request with 
the FAA, formally asking for a copy of reports of fumes and smoke in the cabin/flight 
deck that the FAA received from U.S. airlines occurring from 1 January 2002 through 
31 December 2011 and including one or more of the following search terms: fume, odor, 
smell, smoke, bleed air. On 19 May 2015, the FAA sent the author two spreadsheets, one 
of which contained 15,885 reports that met the author’s criteria and had been reported by 
airlines to comply with the agency’s Service Difficulty Reporting (SDR) regulation (14 CFR 
§ 121.703). The other spreadsheet contained 366 reports that also met the author’s criteria 
and had been reported by airlines to comply with the agency’s Accident and Incident Data 
System (AIDS) rule (Order 8020.11D, 2018). Other than a cursory review of the narrative 
fields, the AIDS data are not summarized in this study. 

For each SDR fume/smoke report, the author reviewed the narrative field (“remarks”) 
and other relevant details to determine the source of the fumes or smoke. The author 
excluded reports that did not involve actual fumes or smoke in the cabin or flight deck, 
and then categorized the remaining reports according to the specified /implied source of 
fumes (Table 1), based on a combination of the text in the narrative field and professional 
judgement. The categories were applied consistently throughout the data set according 
to objective criteria, wherever possible. For some reports, selecting the category involved 
a degree of personal judgement, for example, fumes from a failed or defective fan were 
generally classified as “fan”, even though some of the failures were electrical rather than 
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mechanical. For narratives that described an unfamiliar mechanical condition or fault, 
the author relied on the professional judgement of two senior airline mechanics and two 
commercial airline pilots, respectively. 


Table 1. Types/sources of onboard fumes/smoke. 


Types/Sources of Onboard Fumes/Smoke (N = 15) 


Air conditioning packs (“packs”) 

Battery (smoking, defective) 

Bird strike (into engine) 

Bleed source (specifics undefined) 

Deicing fluid 

Duct (blown, disconnected, clogged) and duct insulation 
Electrical item/system 

Engine oil 

Engine wash 

Fan 

Fuel/exhaust 

Hydraulic fluid 

Oven (in the galley) 

Other, not otherwise classified 

Source not identified /report too vague to classify 


When airline staff submit a fume/smoke report to the FAA’s SDR system, they cite 
the “maintenance manual” chapter assigned to the report internally. These chapters were 
created by the Air Transport Association (ATA) and define both how airline pilots classify 
abnormal conditions in their logbooks, and how maintenance workers classify aircraft 
system troubleshooting and repairs in the maintenance records, all of which are traceable 
to the aircraft tail number. As another source of insight into sources of documented fumes 
and smoke inflight, the author summarized the airline-defined ATA chapters for the whole 
data set and for the three leading sources of fumes. 


3. Results 


Of the 15,885 reports, the author excluded 3468 reports that did not involve an actual 
smoke or fume in the cabin or flight deck. Examples include “lens of smoke goggle was 
cracked” or “test of smoke detector failed”. Those narratives included the word ”smoke”, 
but it is clear that no smoke was present. 


3.1. Types of Onboard Fumes/Smoke Reported by U.S. Airlines 


A total of 12,417 reports remained, including 12,382 reports that were assigned to one 
source category based on the narrative description, and 35 reports that were assigned to 
two categories, creating a total of 12,452 source types (Figure 1). Electrical fumes were the 
most documented type of report (37%). The combination of “bleed-sourced”, “engine oil”, 
and “hydraulic fluid” represented the second most documented type of report (26%). The 
third most documented mechanical failure that generated smoke or fumes was fans (8.5%), 
followed by air conditioning packs (6.4%), and fuel/exhaust (5.2%). A total of 1.8% of the 
reports were “other, not otherwise classified”, which included sources that did not fit the 
other categories, such as engine failures, blown tires, and lightning strikes. Notably, for 
9.1% of the reports, either the mechanical fault that created the fumes/smoke was either 
not identified or the airline supplied insufficient information for the author to reliably 
classify the source (e.g., “fumes reported during climb in the aft cabin”). 
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Figure 1. Sources of onboard fumes/smoke reported by U.S. airlines, 2002-2011 (n = 12,452). 


3.2. ATA Chapters, Whole Dataset 


Airlines classified 89% of the 12,417 reports into the nine ATA chapters listed in 
Figure 2, led by air conditioning (45%), and followed by lights (9.1%), equipment/furnishings 
(8.5%), airborne APU (8.1%), and engine (6.1%). The “other” category (Figure 2) repre- 
sents reports that airlines classified into 29 additional ATA chapters, ranging from one to 
211 reports per chapter. 
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Figure 2. ATA chapters assigned by airlines to fumes/smoke reports, 2002-2011 (n = 12,417). 


3.3. ATA Chapters, Onboard Electrical Fumes/Smoke 


Of the 4573 reports that the author classified as electrical, the most commonly cited 
ATA chapters were lights (24%), air conditioning (15%), and equipment/furnishings (13%). 
The nine ATA chapters for 89.2% of these electrical reports are listed in Figure 3. Airlines 
classified the remaining 10.8% of the electrical reports into one of 24 additional ATA 
chapters (listed as “other” in Figure 3), reflecting the number and diversity of electrical 
systems onboard. 
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Figure 3. ATA chapters assigned by airlines to 4573 reports of electrical fumes /smoke. 
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3.4. ATA Chapters, Onboard Fume Events (Engine Oil, Hydraulic Fluid) 


Of the 3360 reports that the author classified as either engine oil, bleed-sourced, or 
hydraulic fluid, the most commonly cited ATA chapters were air conditioning (51%), 
engine (20%), APU (18%), hydraulic power (3.7%), pneumatic (1.6%), and bleed air (1.1%) 
(Figure 4). Airlines classified the remaining 3.9% of these reports in one of 20 additional 
ATA chapters (listed as “other” in Figure 4). Note that four of the 3360 fume reports were 
classified as both oil and hydraulic fluid because maintenance found both contaminants 
in the bleed system during troubleshooting, but each of those four reports only had one 
assigned ATA chapter. 
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Figure 4. ATA chapters assigned by airlines to 3360 fume events (oil, hydraulic, bleed). 
3.5. ATA Chapters, Onboard Fan-Sourced Fumes, and Smoke 


Failures or malfunctions in a total of 23 types of fans were referenced in the fume and 
smoke reports (Table 2), eight of which appear to be related to the air conditioning system 
and are indicated with bold text in Table 2. Of the 1058 reports that the author classified 
as “fans”, airlines classified the vast majority (96%) as air conditioning (ATA chapter 
21). Airlines classified the remaining 4% of fan failures into one of 11 ATA maintenance 
manual chapters. 
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Table 2. Types of fans referenced in fume/smoke reports. 


Types of Fans (N = 23) Identified as Sources of Fumes or Smoke in Airline Reports 





Air Conditioning Condenser Cooling Fan 
Advisory fan 

Aft rack cooling fan 

APU oil cooling fan 

Attitude heading reference unit fan 
Avionics cooling blower fan 
Avionics vent blower 

Cargo compartment fan 

Cargo SD fans 1 and 2 

Cargo smoke blower 

Chiller fan (forward and aft) 
Equipment cooling fan 

Exhaust fan 

Extract fan 

Flight deck recirculation air fan 
Galley oven fan 

Ground cooling fan 

Heat exchanger cooling fan 
Instrument cooling fan 

Inverter fan 

Lavatory fan 

Radio rack cooling fan 
Recirculated air fan 


3.6. Other Sources of Onboard Fumes/Smoke 


In all, 6.4% of the 12,452 source types for the 12,417 reports (n = 793) were classified 
as “air conditioning packs” (indicated as “packs” in Figure 1) and 13% (n = 1667) were 
classified as “bleed, specifics not defined”. There is more judgement involved with these 
categories than the others because, generally, the pack itself is not as much a source of 
contamination as it is a sink for sources of upstream contaminants. Those contaminants can 
be internal (e.g., leaking oil seal in the engine) or external (e.g., ingested exhaust fumes) 
and the soiling may build up gradually or be associated with a specific upstream failure. 
Additionally, the pack contains electrical components, which can fail. For these types of air 
supply system-sourced fumes, the author and, often, one of the selected experts (two pilots 
and two senior mechanics) had to make a judgement call as to whether the fumes were 
more likely sourced to the pack, bleed air upstream of the pack, or another external source. 
Generally, the “packs” category was chosen if a mechanical fault was later found in the 
pack. Such mechanical faults were typically not specified but were suggested because all 
or part of the pack was replaced. Conversely, the “bleed, specifics not defined” category 
was chosen if the fumes were reported during an engine power setting change (usually 
take-off/climb or descent), no mechanical fault was subsequently found (or inferred) in the 
pack, and an external source (such as exhaust fumes) was not referenced. A small subset of 
pack failures were classified as “engine oil” if the description of fumes was consistent with 
oil and the air cycle machine bearings in the pack would have been oil-lubricated. 

In all, 1.8% of the 12,452 source types (n = 224) were classified as “other, not otherwise 
classified”. The sources in this category included engine failures (unless oil fumes were 
specified), supply air from a ground cart, lightning strikes, blown tires, lavatory fluid, and 
maintenance products, such as grease or solvents. 

Only 10 of the 12,417 fume/smoke reports cited odors that were sourced to either 
lavatory fluid or waste (e.g., servicing issue, blocked drain, leaking lav seals), so these were 
not assigned an independent category. Additionally, despite anecdotal reports that nail 
polish is a significant source of onboard fumes, there were only six reports that mentioned 
the smell of either nail polish or acetone, and the maintenance response indicated that none 
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of those reports were sourced to either product. Presumably, lavatory odors and nail polish 
are easily identifiable and, generally, do not require maintenance intervention. 


4. Discussion 
4.1. Airline Reporting Rules for Onboard Fumes/Smoke 


In the U.S., commercial airlines are required to report each “failure, malfunction, or 
defect” that causes “smoke, vapor, toxic or noxious fumes” to accumulate or circulate in the 
flight deck or cabin during flight (defined as “wheels up”) (14 CFR § 121.703(a)(5)). Airlines 
are also required to report ground-based events if the airline is of the “opinion” that flight 
safety could have been “endangered” had the aircraft continued (14 CFR § 121.703(c)). 
There are comparable SDR rules for charter flights, too (14 CFR § 135.415). These “Service 
Difficulty Reporting” (SDR) rules will not capture fume events for which the causal defect 
is not apparent during a ground-based inspection. For example, oil fumes that contaminate 
the bleed air through a worn engine part subjected to thermal or physical stress during 
an engine power setting change inflight need not be reported if that same part does not 
leak under ground-based conditions during a maintenance inspection. As well, when 
maintenance workers overfill the oil or hydraulic fluid reservoir, the subsequent spillage 
into the bleed air is a recognized source of fumes [24—26], but spillage is not a mechanical 
defect, so does not need to be reported. 

In this dataset, 9.1% (n = 1135) of the SDR reports from airlines to the FAA do not 
offer insight into the relevant mechanical defect or failure. It is not clear if the fault is not 
described or not known when the report is submitted. Technically, though, it is the defect 
that is reportable, not the presence of fumes. Even though, there may not be an actual 
defect and, even when there is, maintenance troubleshooting can be “trial and error” [27], 
which can delay defect identification. 

In addition to the SDR regulation for reporting onboard fumes and smoke, U.S. airlines 
must also report each “interruption to a scheduled flight” caused by a known or suspected 
(emphasis added) mechanical difficulty or malfunction, which is not required to be reported 
under the SDR rule (14 CFR § 121.705). However, it is difficult to assess the degree to which 
airlines comply with this rule because there is no central, searchable FAA database; rather, 
these “mechanical interruption reports” are maintained by individual FAA Certificate 
Management Offices, and only for one year. 

Finally, U.S. airlines must also report “occurrences which affect or could affect safety 
of operations”, including onboard fumes/smoke, to the FAA’s “Accident and Incident 
Data System” (AIDS) per Order 8020.11D (2018). The author’s 2014 FOIA request included 
AIDS reports with one or more of the same search terms requested for the SDR data, but 
U.S. airlines only submitted 2.3% as many fume/smoke reports to meet Order 8020.11D 
as compared to the SDR regulation (n = 366 compared to 15,885) during the same time 
period. It would be reasonable to assume that AIDS reports are the more serious events 
during which safety of operations was (or could have been) compromised. To test this 
hypothesis, the author compared the contents in the narrative fields of the AIDS and SDR 
datasets (which were selected based on the same set of search terms during the same time 
frame) and concluded that the nature and seriousness of the reports is comparable. It 
appears that the AIDS regulation is either ignored or misunderstood, and that compliance 
is not enforced. 

In summary, then, if airline staff do not identify a mechanical source of reported fumes, 
or if the fumes occur before take-off, then the event is not reportable. If airline staff do 
identify a mechanical source of fumes, then the airline must report to the SDR database. If 
a mechanical source was suspected and the flight plan was affected, then the airline must 
report to the “mechanical interruption” database. In both of those scenarios, if safety either 
was or could have been compromised, then the airline must also report to the AIDS system. 
However, ground-based fume events can still compromise safety, fumes can occur without 
obvious mechanical defects, and an airline is unlikely to report the same event to two 
reporting systems. As a result, these reports of fumes and smoke are only a subset of the 


Aerospace 2021, 5, 122 


8 of 14 


actual events, and the data that do exist get scattered, all making an accurate assessment 
of the relative and absolute frequencies of different sources of onboard fumes and smoke 
even more difficult. 

In addition to the complying with FAA reporting rules, U.S. airlines must report fume 
and smoke events to the U.S. National Transportation Safety Board (NTSB), if the event 
meets the NTSB definition of either an “accident” or “serious incident” (49 CFR § 830.5, 
49 CFR § 830.2). Some fume events meet the NTSB reporting criteria, such as when a 
crewmember is hospitalized for more than 48 h during the seven days post-flight, sustains 
internal organ damage, or is unable to perform normal flight duties. 

Another source of fume event data for the U.S. fleet is the “Aviation Safety Report- 
ing System” (ASRS) database, hosted by the National Aeronautics and Space Adminis- 
tration (NASA). ASRS catalogues aircraft incident reports for a broad range of unsafe 
conditions from pilots, cabin crew, maintenance workers, air traffic controllers, and dis- 
patchers. Reporting is voluntary and confidential, with the goal of circulating infor- 
mation about safety issues within the industry to prevent reoccurrences. Of 341 ASRS 
“smoke/fire/fumes/odor” events that crewmembers reported in 2016-2017, 35% were 
sourced to the air supply system and 20% were sourced to electrical faults [28]. 


4.2. Underreporting Engine Oil and Hydraulic Fume Events 


Unfortunately, even within the narrowed scope of the reporting rules for fume and 
smoke events, there is evidence of airline underreporting. Some of the underreporting 
is cited for fume events in general. For example, in 2006, the FAA issued a bulletin for 
its maintenance inspectors, acknowledging poor airline compliance with the fume event 
reporting regulations. Specifically, the agency noted that “there have been concerns raised 
about numerous reports of [fume] events on commercial air carrier/operator aircraft. 
During the FAA’s analysis of this [sic] data, it appears as though there are numerous 
[airlines] who may not have reported these events as required by regulation” [21]. 

There is evidence, however, that some underreporting is specific to oil and hydraulic 
fume events, reducing both the number of reports and the relative frequency of those 
events. For example, a review of fume events consistent with oil and hydraulic fluid 
sources at one U.S. airline over 18 months in 2009-2010 found that the airline had only 
reported 38% of 87 fume events documented by crewmembers to the FAA SDR system [11]. 
Part of this discrepancy is explained by the airline not consistently complying with the 
reporting rules, but part of it is explained by regulatory language that does not require 
airlines to report a significant fraction of documented events, as described above. 

An example of underreporting specific to engine oil fumes is a 2013 report that the 
FAA sent to two Congressional committees regarding a law requiring the FAA to conduct 
research and development into sensor and filtration technology, specific to preventing 
exposure to oil fumes in the bleed air on commercial aircraft [29]. The FAA reported to 
Congress that officials had “mined” their fume event databases from 2002-2011, inclusive, 
for reports that included one or more of the search terms “fume, odor, smell, smoke, and 
bleed air”, and that U.S. airlines had only reported 69 oil fume events and no hydraulic 
fluid events during those 10 years [23]. The FAA told Congress that, since the occurrence 
of these events was “extremely low”, it could not justify spending its resources to research 
and develop bleed air sensors and filters. However, this analysis of those same cited reports 
from just one of the FAA databases using the same search terms identified oil as a source of 
fumes in 1353 reports, hydraulic fluid as a source in 173 reports, and bleed air contaminants 
consistent with oil/hydraulic fluid in 1667 reports. This type of underreporting reduces 
the relative frequency of these reports and downplays the need for relevant regulatory 
control measures. 

Finally, underreporting oil and hydraulic fluid fume events is not limited to U.S. 
airlines. In 2011, the European Aviation Safety Agency (EASA) published comments it had 
received in response to the question of whether a rulemaking to address oil fumes in the 
cabin and flight would be justified. EASA noted multiple claims by advocates for pilots, 
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cabin crew, and passengers that European airlines underreport (a word used 53 times in the 
report) fume events and that the safety impacts are downplayed [30]. Likewise, a report 
by the Australian Transport Safety Bureau (ATSB) noted evidence that Australian airlines 
underreport fume events, both to the ATSB and to the aviation regulator [31]. An Australian 
Senate Enquiry into fume events on one particularly problematic aircraft type noted “strong 
evidence of a tendency of pilots to under-report [fume events]”, citing a combination of 
pressure from airline management and lack of training [32]. The International Civil Aviation 
Organization (ICAO) has formally recommended that airlines “encourage flight and cabin 
crew members to report fume events”, referring to crew reports as “essential to assist airline 
maintenance technicians in identifying the root cause of an event and taking corrective 
action to prevent reoccurrence” [14]. Currently, though, fume event reporting to aviation 
regulators is largely limited to maintenance defect reporting. 


4.3, Estimates of the Frequency of Engine Oil and Hydraulic Events 


The answer to the question of how often oil and hydraulic fluid fume events happen 
during commercial airline flights is a subject of considerable debate within the industry. It 
is important to get an accurate answer because, as the FAA’s 2013 report to Congress illus- 
trates [23], the answer will influence the regulatory agenda. The FAA itself acknowledges 
this reality in a 2016 bulletin in which the agency notified airlines “that some data required 
by [the SDR rule] is not being entered properly” and that “[b]ecause of poor data integrity, 
the FAA, manufacturers, and air carriers are unable to accurately detect trends necessary to 
proactively mitigate risk” [22]. 

One of the first estimates of how often oil and hydraulic fluid fume events occur 
was published in a 2002 U.S. National Research Council (NRC) committee report about 
aspects of aircraft air quality associated with potentially negative impacts on crew and 
passenger health. The report referenced fume event data collected at three airlines “over 
several years”, which ranged from 0.9 events per 10,000 flights on Boeing 737 aircraft up to 
39 events per 10,000 flights for the British Aerospace 146 [33]. The NRC committee issued a 
series of recommendations to the FAA, including that the agency “rigorously demonstrate” 
the “adequacy” of its cabin air regulations. In response, the FAA acknowledged that its 
“rulemaking has not kept pace with public expectation and concern about air quality and 
... [nJo present airplane design ... incorporates an air contaminant monitoring system to 
ensure that the air provided to the occupants is free of hazardous contaminants... ” [34]. 
As of this writing, the FAA has still not initiated a rulemaking on this subject. 

In 2010, the FAA published an informational bulletin for airlines, stating that the 
agency “ ... continues to receive over 900 reports a year on [all sources of] smoke or fumes 
in the cabin and or cockpit” [35]. Based on the number of flights in 2010 [36], this translates 
into 0.9 events per 10,000 flights. The FAA noted that, “in fact, it is not unusual to receive 
more than one report during a 24-h period. For instance, on one day in April of 2010, five 
reports of smoke in the cockpit came in from one [airline]. All these incidents prompted 
the flightcrew to declare emergencies and divert to the nearest airport” [35]. 

In 2013, the International Air Transport Association (IATA) reported on its analysis of 
operational safety reports from more than 150 airlines covering almost 31 million flights 
during the period 2008-2012 [37]. Of these, 3444 reports were categorized as “smoke, 
fumes, and odor”, which is equivalent to about one reported event per 10,000 flight 
cycles. Engine-related fumes were the most commonly reported source (23%), followed 
by electrical-related fumes (15%). Maintenance action was the most common operational 
effect (25%). 

In 2015, FAA-funded researchers published their review of engine oil and hydraulic 
fluid fume events identified in FAA, NTSB, and NASA databases from 2007—2012 and sorted 
them by aircraft type. They concluded that U.S. airlines reported fume events on almost 
every aircraft type at an average frequency of two events per 10,000 flights [38], which 
translates into an average of 5.3 events per day based on the number of flights operated 
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during that time [36]. The authors reported “substantial variation” in the frequency of 
reports by aircraft type, with the maximum being eight reported events per 10,000 flights. 


4.4, Recommendations to Improve Data Collection for Onboard Fumes/Smoke 


Beyond the need for better compliance and enforcement of existing regulations, this 
review of SDR data from 2002-2011 provides insight into how fume reporting regulations 
could be improved. 

First, SDR data would be more representative if airlines were required to report the 
presence of onboard fumes/smoke during ground operations, not just inflight. Industry 
sources confirm that the APU (which is the typical bleed air source during ground opera- 
tions, whether at or away from the gate) is a source of oil and hydraulic fluid fumes [25]. 
Reports confirm that APU oil fumes, which contaminate the cabin or flight deck air during 
eround operations, can generate reports of crew ill health [11]. Additionally, if fumes that 
manifest on the ground are not resolved, they can continue inflight [25], which can affect 
both the safety of operations [26] and the health of the crew and passengers [11]. For these 
reasons, ground-based fume events should be documented and prevented. 

Second, reported onboard fumes/smoke for which there is no obvious mechanical 
fault should still be reportable. For example, overfilling an oil or hydraulic fluid reservoir 
is a recognized source of fumes [26,39]. Additionally, some faults (such as oil leakage 
through a worn but not defective engine seal) that introduce fumes to the cabin air may 
only happen during the stress of a high-engine power setting, for example, and may not be 
recreated during routine ground-based troubleshooting [40]. Maintenance staff may not 
identify the relevant mechanical failure even for serious events, such as those involving 
pilot impairment and incapacitation [6]. SDR data would be more representative and useful 
if the presence of fumes/smoke caused by even a suspected defect or failure (defined to 
include spills and leaks) was reportable. 

Third, it would be easier for SDR database users to answer the “how often” and “what 
type” questions regarding onboard fumes/smoke if the FAA’s report categories were more 
tailored to the actual sources of fumes (see Figure 1), or if airlines could choose more 
than one category. For example, a pull-down menu for airlines to select standardized 
responses to the question of “source” would assist with subsequent analyses. The narrative 
field for each report can provide useful details but is too technical to be interpreted by 
software. Currently, when an airline submits a report, they must choose the most rele- 
vant category (“nature condition”), as defined by the FAA. One of the FAA categories is 
“smoke/fumes/odors/sparks”, but only 76% of the 3360 events the author classified as 
“oil”, “hydraulic”, or “bleed-sourced” were categorized as such. U.S. airlines categorized 
another 16% of those bleed air reports as “fluid loss” and assigned the remaining 8% to one 
of 13 additional categories. For the full dataset (n = 12,417, all sources), airlines classified 
84% of the reports as “smoke/fumes/odor/sparks” and 5.4% as fluid loss. In both cases, 
selecting only the “smoke/fumes/odor/sparks” category to search for fume and smoke 
reports seems reasonable but would underestimate the number of relevant reports. 

Lastly, a crewmember reporting system for fume events to supplement the existing 
maintenance reporting systems is needed and must be accompanied by suitable training 
programs to ensure that crewmembers know how to recognize and respond to different 
types of onboard fumes and smoke [14]. This point has been raised since at least 2002, 
starting with an NRC committee, which recommended that the FAA establish a crewmem- 
ber reporting system to collect air quality event data and associated symptom reports [33]. 
The following year, Congress passed a law requiring the FAA “to establish an air quality 
incident reporting system” [41]. Nine years later, Congress again required the FAA to “es- 
tablish a systematic reporting standard for smoke and fume events in aircraft cabins” [42]. 
In 2017, the Cabin Air Safety Act of 2017 was introduced to the floor of the U.S. Senate, 
and the Know Before You Fly Act was introduced to the floor of the U.S. House. Again, 
these bills called on the FAA to establish a system for crewmembers and mechanics to 
report fume events on commercial aircraft, but neither bill was passed. In early 2018, the 
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FAA recommended that airlines “assess current policy and procedures” regarding their 
procedures to differentiate and mitigate smoke and fumes, but the FAA did not mandate 
event reporting [43]. Most recently, in October 2018, Congress mandated that the FAA 
“issue guidance for airline crews and maintenance technicians to report incidents of smoke 
or fumes on board an aircraft operated by a commercial air carrier” [44]. In March 2021, the 
FAA launched a cabin air quality website that includes a link to a voluntary, deidentified 
reporting system operated by NASA and a link to the FAA’s maintenance defect reporting 
system [45]. However, the FAA has still not implemented crewmember reporting systems 
for fume events on commercial aircraft. 


5. Conclusions 


Of the onboard fumes/smoke that U.S. airlines reported to the regulatory agency 
over a 10-year period, electrical-sourced fumes and smoke were the most prevalent (37%), 
reflecting the number and diversity of onboard electrical systems and associated faults 
(Figures | and 3). 

Bleed air issues (i.e., either confirmed or suspected engine oil and hydraulic fluid 
fumes) were the second most prevalent type of reported fumes/smoke (26%), and air- 
lines internally classified 89% of those reports as either air conditioning, engine, or APU 
(Figures 1 and 4). Oil fumes were reported almost eight times more frequently than hy- 
draulic fluid fumes (Figure 1). 

Fans were the third most prevalent source of reported fumes/smoke (8.5%), most 
of them related to the air conditioning systems (Figure 1, Table 2). It would be worth 
investigating how to either improve the design of these fans or to implement an inspection 
program (e.g., measure vibration or electrical current draw) as preventive maintenance. 

It is important to note that, although these data were reported by U.S. airlines, the 
harmonized nature of regulations that govern aircraft design, maintenance, and operating 
practices means that the data likely also reflect the onboard sources of fumes and smoke 
reported in other countries with equivalent aviation systems. 

To more accurately count the number and types of fumes/smoke on U.S.-registered 
aircraft, it would be helpful if the FAA improved its current maintenance reporting reg- 
ulations by requiring airlines to report both ground-based events and those for which a 
mechanical fault is not readily identified. In addition, there is a need for the FAA—and 
its equivalent regulatory bodies around the world—to implement a reporting system for 
crewmembers to gather data on the health and safety implications of exposure to air supply 
system-sourced fumes onboard, as per the recommendation of the ICAO [14]. 

Finally, the proportion of reported fume events that involve engine oil/hydraulic fluid 
described here, in conjunction with a more expansive review of the frequency of reported 
fume events on US airlines [38] and reports that flight safety can be compromised (cited 
above) support the installation and operation of equipment to filter and monitor bleed 
air contaminants in real time when aircraft are in operation. In the U.S., research into 
suitable technologies was most recently called for in legislation passed in 2012 [29] and 
2018 [44]. Progress is slow, but there is progress, nonetheless; in late 2020, an FAA-funded 
research consortium held its first meeting on technologies intended to monitor bleed air 
contaminants and detect events in real time. 
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APU Auxiliary Power Unit 

ASRS _ Aviation Safety Reporting System 
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ATSB Australian Transport Safety Bureau 

FAA United States Federal Aviation Administration 

NASA United States National Aeronautics and Space Administration 
NRC United States National Research Council 

NTSB United States National Transportation Safety Board 
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